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ABSTRACT 


As  the  Marine  Corps  pursues  greater  energy  efficiency  in  expeditionary  operations,  the 
HOMER  micropower  optimization  model  provides  potential  to  serve  as  a  powerful  tool 
for  improving  Marine  Corps  power  planning.  The  HOMER  software  was  developed  for 
the  modeling  and  simulation  of  micropower  systems  over  long  periods  of  time.  Although 
a  deterministic  model,  HOMER  uses  stochastic  input  data,  specifically  solar  irradiance, 
temperature,  and  load  profiles.  HOMER  simulation  fidelity  is  therefore  affected  by  the 
inter-annual  variability  of  these  profiles.  This  research  quantifies  HOMER  robustness  with 
regard  to  solar  irradiance  and  temperature  profile  variability  through  full-factorial  experi¬ 
mental  designs.  The  effect  of  shortening  HOMER  simulation  duration  on  the  variability 
of  HOMER  simulation  outputs  is  also  investigated,  and  though  statistically  significant,  the 
resulting  increase  in  variability  is  not  large  enough  to  preclude  the  use  of  HOMER  for  expe¬ 
ditionary  operations.  This  thesis  also  demonstrates  how  HOMER  can  assist  in  developing 
power  planning  doctrine,  showing  that  the  fuel  consumption  benefits  of  using  multiple  gen¬ 
erators  of  different  sizes  is  no  longer  present  once  a  renewable  energy  asset  is  added  to  the 
micropower  system.  This  analysis  of  HOMER’s  robustness  and  operational  potential  pro¬ 
vides  insight  for  improving  the  Marine  Corps’  use  of  HOMER  for  power  planning  in  an 
expeditionary  environment. 
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CHAPTER  1: 

Thesis  Introduction  and  Discussion  of  Marine  Corps 

Energy 


1.1  Introduction 

This  research  provides  analysis  of  the  Hybrid  Optimization  Model  for  Electric  Renewables’ 
(HOMER)  robustness  and  potential  application  for  Marine  Corps  expeditionary  operations. 
HOMER  is  one  of  the  few  models  in  existence  for  the  modeling  of  micropower  systems 
that  include  both  conventional  and  renewable  energy  resources.  The  model  is  a  time-step 
model  developed  to  analyze  long-term  life-cycle  performance  of  different  micropower  sys¬ 
tems.  This  research  will  investigate  areas  of  concern  in  the  modeling  of  renewable  energy 
power  planning,  with  an  emphasis  on  solar  power.  The  research  will  include  considera¬ 
tion  of  the  impact  of  inter-annual  variability  in  solar  rrradiance  and  temperature  profiles 
on  the  robustness  of  the  deterministic  HOMER  model.  HOMER’s  having  been  developed 
and  used  primarily  for  evaluation  of  micropower  system  long-term  performance  brings  into 
question  whether  the  model  and  simulation  software  is  robust  enough  for  simulating  rela¬ 
tively  short-term  Marine  Corps  expeditionary  operations.  Through  analysis  of  these  areas, 
this  research  will  assist  the  Marine  Corps  in  its  validation  of  the  HOMER  model  as  a  tool 
for  power  planning  in  expeditionary  operations. 

1.1.1  Scope 

The  purpose  of  this  research  is  to  conduct  a  formal  study  of  the  HOMER  model  as  it  is 
tailored  for  the  Marine  Corps,  to  analyze  components  of  the  model,  and  evaluate  the  ro¬ 
bustness  of  the  model.  Variability  over  solar  irradiance  and  temperature  profile  estimates 
are  expected  to  be  influential  on  the  performance  of  HOMER  supporting  expeditionary 
operations.  This  research  assesses  the  usefulness  of  HOMER  for  Marine  logisticians,  engi¬ 
neers,  and  acquisition  professionals  alike  by  evaluating  the  model’s  robustness  with  regard 
to  inter-annual  solar  irradiance  and  temperature  profile  variability  across  multiple  factors. 
These  factors  include  varying  the  duration,  season,  and  location  of  simulated  micropower 
systems,  where  the  term  "micropower  system"  refers  to  "a  system  that  generates  electric- 
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ity...to  serve  a  nearby  load  [1]."  The  faetors  also  include  varying  the  size  of  the  electric 
loads  being  powered  and  the  composition  of  generators  powering  the  loads.  Along  with 
this  analysis  of  model  robustness,  this  thesis  provides  recommendations  on  Marine  Corps 
power  planning  doctrine.  These  additional  recommendations  provide  unique  insights  on 
the  benefits  of  grid-tie  capable  converters  and  the  benefits  of  using  multiple  generators  of 
different  sizes  to  power  a  load. 

The  scope  of  this  research  can  be  more  specifically  understood  by  reviewing  the  research 
questions  for  which  it  provides  answers.  These  research  questions  are: 


1.  Is  the  HOMER  software  robust  enough  for  modeling  and  simulating  micropower 
system  performance  for  Marine  Corps  expeditionary  operations,  given  its  primary 
use  for  the  simulation  of  long-term  micropower  system  performance? 

2.  How  does  inter-annual  variability  of  solar  irradiance  for  a  specified  location  and  time 
frame  impact  the  HOMER  model’s  effectiveness?  How  do  each  of  the  five  factors 
previously  mentioned  influence  the  relationship  between  variability  in  HOMER  esti¬ 
mates  and  inter-annual  solar  irradiance  variability? 

3.  How  does  the  variation  between  National  Renewable  Eaboratory  typical  meteoro¬ 
logical  year  solar  irradiance  profiles  and  measured  annual  solar  irradiance  profiles 
influence  the  fidelity  of  HOMER  simulations? 

4.  How  does  inter-annual  variability  of  the  temperature  for  a  specified  location  and  time 
frame  impact  the  HOMER  model’s  effectiveness?  How  do  each  of  the  five  factors 
influence  the  relationship  between  variability  in  HOMER  estimates  and  inter-annual 
temperature  variability? 

5.  Does  a  micropower  system  composed  of  multiple  generators  of  different  sizes  con¬ 
sume  less  fuel  than  a  micropower  system  composed  of  a  single  peak  load  capable 
generator  when  powering  the  same  load?  Does  this  change  when  a  renewable  energy 
power  system  is  added  to  the  micropower  system? 

6.  Does  a  grid-tie  capable  renewable  energy  power  system  provide  greater  fuel  savings 
than  a  renewable  energy  power  system  that  is  not  grid-tie  capable? 
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1.1.2  Methodology 

As  will  be  discussed  further  in  Chapter  4,  this  research  was  conducted  through  the  exe¬ 
cution  of  two  designs  of  experiments  (DOE),  utilizing  the  NFS  SEED  Center  capabilities 
to  support  analysis.  The  first  DOE  focused  on  analyzing  the  robustness  of  HOMER  with 
regard  to  solar  irradiance  values.  Eor  this  DOE,  solar  irradiance  profiles  from  10  locations 
across  the  United  States  were  acquired  and  included  as  inputs  to  the  HOMER  model  as 
varying  levels  of  the  "location"  factor.  Different  load  profiles  were  built  for  each  of  the 
load  size  factor  levels  and,  through  the  use  of  the  SEED  Center,  HOMER  simulations  were 
conducted  for  the  completion  of  a  full  factorial  analysis. 

The  second  design  of  experiments  focused  on  identifying  HOMER’s  robustness  with  regard 
to  temperature  profile  robustness.  Eor  this  analysis,  inter-annual  variability  of  temperature 
profiles  was  analyzed  similar  to  how  the  solar  irradiance  variability  was  analyzed  in  the  first 
DOE.  Eollowing  this  DOE,  the  results  of  the  two  DOEs  were  combined  to  identify  the  influ¬ 
ence  of  inter-annual  variability  of  solar  irradiance  temperature  profiles  on  the  variability  of 
renewable  energy-derived  fuel  savings  estimates,  estimated  through  HOMER  simulations. 

Eollowing  analysis  of  the  influence  of  inter-annual  variability  of  solar  irradiance  and  tem¬ 
perature  profiles  on  the  HOMER  model,  analysis  was  conducted  on  benefits  of  grid-tie 
capable  converters  and  the  benefits  of  using  multiple  different  sized  generators  to  power  a 
load.  This  additional  analysis  was  conducted  through  the  execution  of  additional  HOMER 
simulations  with  an  added  "generator  composition"  factor  level  and  comparison  of  previ¬ 
ous  DOE  simulation  output.  Now  that  an  introduction  to  this  research  has  been  provided, 
background  information  on  the  Marine  Corps  energy  requirements  and  the  modeling  and 
simulation  challenge  must  be  explained  before  the  methodology,  analysis,  and  findings  of 
this  research  can  be  presented  in  detail. 

1.2  Marine  Corps  Energy  Challenges 

Operation  Enduring  Ereedom  (OEE)  and  Operation  Iraqi  Ereedom  (OIE)  demonstrated  to 
the  Marine  Corps  that  increased  employment  of  technology  on  the  battlefield  is  paired  with 
a  serious  challenge  to  expeditionary  capabilities.  This  challenge  comes  in  the  form  of 
an  over-reliance  on  fuel,  which  shortens  the  tether  between  a  unit’s  mobility  and  its  fuel 
resupply  requirements. 
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At  the  beginning  of  OEF  and  OIF,  Marine  Corps  units  deployed  to  expeditionary  eom- 
bat  zones  employed  Marine  Corps  generators  to  meet  all  base  power  requirements.  These 
power  requirements  ranged  from  lights  and  computers  in  Combat  Operation  Centers  (COC) 
to  camera  systems  and  floodlights  employed  in  support  of  base  security.  Over  the  course  of 
the  Iraq  and  Afghanistan  wars,  the  Marine  Corps  saw  an  increased  energy  consumption  by 
Marine  Corps  units.  In  2011,  the  Marine  Corps  Expeditionary  Energy  Office  (E20)  stated 
that  the  amount  of  fuel  required  for  an  average  Marine  company  to  execute  a  three-day 
patrol  had  increased  from  160  Watts  in  2001  to  1255  Watts  in  2011,  an  increase  of  684  per¬ 
cent  [2].  This  increased  energy  requirement  resulted  in  a  subsequent  increased  reliance  on 
internal  combustion  generators  and  diesel  fuel.  The  increased  energy  consumption  was  due 
to  increases  in  warfighting  capabilities  with  "enhancements  to  command,  control,  commu¬ 
nications,  computers,  and  intelligence  (C4I)  technologies;  hardened  vehicles,  and  weapon 
systems  [3]."  From  2001  to  201 1,  the  Marine  Air  Ground  Task  Force  (MAGTF)  number  of 
radios,  computers,  and  vehicles  increased  by  250,  300,  and  200  percent,  respectively  [3]. 
In  the  extreme  Afghanistan  climates  60  percent  of  Marine  Corps  power  production  was 
dedicated  to  just  keeping  these  C4I  systems  cool  enough  to  operate  [3]. 

To  support  the  energy  requirements  for  these  technology  advances  and  equipment  in¬ 
creases,  the  number  of  internal-combustion  generators  in  the  Marine  Corps  tripled  between 
2001  and  2013,  rising  from  4,000  to  more  than  13,000  [4].  The  amount  of  fuel  consumed 
by  deployed  Marine  Corps  units  also  increased  linearly  as  the  dependency  on  generators 
increased. 

Increased  fuel  consumption  needs  of  a  given  unit  can  be  met  either  by  increasing  the  fuel 
capacity  of  the  unit  or  by  increasing  the  number  of  fuel  resupply  missions  executed  by  sup¬ 
porting  units.  Increasing  capacity  can  decrease  the  unit’s  maneuverability  and  is  often  an 
untenable  solution  due  to  manning  and  equipment  limitations.  That  is,  battlefield  mobility 
requirements  limit  the  fuel  capacity  of  individual  Marine  Corps  unit  positions.  Units  that 
require  a  capability  of  being  mobile  can  only  store  as  much  fuel  as  they  can  carry  with 
them  on  fuel  transport  trucks  when  it  comes  time  to  move.  The  amount  of  fuel  consumed 
by  a  unit,  therefore,  directly  relates  to  the  number  of  logistics  convoys  required  for  support. 
To  give  perspective  to  the  transport  requirements  for  fuel,  one  7-ton  truck  can  carry  800 
gallons  of  fuel  and  requires  two  to  three  Marines  to  man. 
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Increasing  the  number  of  fuel  resupply  missions  required  ean  tether  the  unit  more  to  the 
supporting  units  and  similarly  deerease  maneuverability.  Inereasing  the  number  of  fuel 
resupply  missions  to  support  a  unit  also  inereases  foree  proteetion  risks.  In  the  reeent 
Afghanistan  and  Iraq  wars,  this  foree  proteetion  risk  was  manifested  in  the  improvised 
explosive  deviees  (lED)  employed  by  insurgents  against  eoalition  eonvoys.  In  2012,  the 
Seeretary  of  the  Navy  reported  to  Congress  that  "for  every  50  fuel  eonvoys  in  theater,  there 
is  one  Marine  easualty  [5]."  Reeognizing  the  importanee  of  the  energy  issue,  the  Comman¬ 
dant  of  the  Marine  Corps,  General  James  T.  Conway,  stood  up  the  Expeditionary  Energy 
Offiee  in  2009  to  "analyze,  develop,  and  direet  the  Marine  Corps’  energy  strategy  in  order 
to  optimize  expeditionary  eapabilities  aeross  all  warfighting  funetions  [6]."  Two  years  later 
General  Conway’s  sueeessor.  General  James  E.  Amos,  presented  a  Marine  Corps  strategy 
to  taekle  the  ehallenges  fossil-fuel  dependeney  brings  to  the  Marine  Corps’  expeditionary 
nature. 


1.3  The  Marine  Corps  Expeditionary  Energy  Strategy 

General  Amos  published  the  Marine  Corps  Expeditionary  Energy  Strategy  on  Mareh  21, 
2011.  In  this  energy  strategy,  the  Commandant  of  the  Marine  Corps  set  an  objeetive  of 
making  50  pereent  of  Marine  Corps  bases  “net- zero  energy  consumers”  by  2020  [7].  In 
2011,  fuel  eonsumption  in  Afghanistan  was  estimated  at  8  gallons  per  Marine  per  day  [8]. 
In  the  Expeditionary  Energy  Strategy,  General  Amos  foeused  on  the  deployed  environment 
and  ealled  for  a  50  pereent  reduetion  in  fuel  consumption  per  Marine  down  to  4  gallons  per 
Marine  per  day  [3] .  The  Marine  Corps  reduetion  in  fuel  eonsumption  will  result  not  only 
in  inereased  maneuverability  and  sustainability,  but  also  a  reduetion  in  fuel-supply  related 
easualties. 

To  aeeomplish  this  reduetion  in  fuel  eonsumption,  the  Marine  Corps  first  foeused  on  redue- 
ing  the  fuel  consumption  of  forward  operating  bases  (EOB).  The  Marine  Energy  Assess¬ 
ment  team  identified  that  32  pereent  of  fuel  eonsumed  by  the  Marine  Expeditionary  Brigade 
(MEB)  deployed  to  Afghanistan  in  2009  was  used  for  eleetrie  power  generation  [9]. 

1.3.1  EXFOB  and  Marine  Corps  Renewable  Energy  Technology 

Through  the  efforts  of  E20,  the  Offiee  of  Naval  Researeh  (ONR),  and  their  interaction 
with  industry,  the  Marine  Corps  has  developed  robust,  renewable  teehnologies  to  meet  the 


5 


energy  objeetives  laid  out  in  the  Commandant’s  energy  strategy.  Through  E20,  the  Marine 
Corps  has  investigated  energy  solutions  in  sueh  areas  as  photo-voltaie  (PV)  arrays  for  power 
harvesting,  light  emitting  diodes  (LED)  for  deereased  energy  eonsumption,  and  improved 
battery  and  smart  power  eontroller  teehnologies  for  better  energy  management. 

One  forum  that  has  been  partieularly  useful  in  advaneing  Marine  Corps  renewable  energy 
researeh  and  development  has  been  the  Marine  Corps’  Experimental  Eorward  Operating 
Base  (EXEOB).  The  EXEOB  is  an  annual  event  established  by  the  Commandant  of  the 
Marine  Corps  in  2009  that  "brings  together  stakeholders  from  aeross  the  Marine  Corps  re¬ 
quirements,  aequisition,  and  teehnology  development  eommunities  in  a  dynamie  proeess 
to  quiekly  evaluate  and  deploy  teehnologies  that  reduee  our  need  for  ’liquid  logisties’  to¬ 
day,  and  to  establish  requirements  for  tomorrow  [10]."  In  the  annual  EXEOB  Request  Eor 
Information  (REI)  doeument,  E20  identifies  for  industry  the  partieular  area(s)  of  interest 
for  a  given  year’s  EXEOB.  The  2014  EXEOB  was  announeed  as  foeusing  on  taetieal  en¬ 
ergy  harvesting.  In  2009  the  first  EXEOB  immediately  began  yielding  renewable  energy 
teehnologies  that  were  rapidly  and  sueeessfully  pushed  through  development,  testing,  and 
fielding  to  the  war  fighter.  Two  of  the  E20’s  most  notable  and  widely  deployed  teehnolo¬ 
gies,  the  Ground  Renewable  Expeditionary  Energy  System  (GREENS)  and  Spaee  Portable 
Alternative  Communieations  Energy  System  (SPACES)  systems,  were  part  of  this  first  EX¬ 
EOB. 

1.3.2  Marine  Corps  PV  Systems 

One  of  the  most  sueeessful  renewable  energy  assets  developed  by  the  Marine  Corps  is  the 
GREENS.  The  GREENS  is  a  power  generation  and  eonversion  system  that  allows  Marines 
to  power  systems  with  solar  energy.  Eaeh  GREENS  is  eomposed  of  eight  PV  array  panels, 
four  High  Energy  Lithium  Batteries  (HELB),  and  a  eentral  eontroller.  Eigure  1.1  shows 
how  all  eomponents  of  a  GREENS  are  eombined  to  harvest  solar  power,  reeeive  alternating 
eurrent  (AC)  and  direet  eurrent  (DC)  power,  employ  the  HELB  batteries,  and  to  provide  DC 
output  for  a  load.  While  a  single  GREENS  ean  provide  300  watts  of  power  for  24  hours 
or  1000  watts  of  peak  power,  GREENS  systems  are  modular  and  ean  be  eombined  into  as 
large  as  a  five  GREENS  system,  providing  up  to  1500  watts  of  eontinuous  power  [11].  The 
GREENS’  modular  design  also  allows  the  user  to  employ  fewer  than  the  full  8  solar  panels 
and  fewer  than  the  full  four  HELB  batteries  if  desired,  to  eonfigure  the  GREENS  to  meet 
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Figure  1.1:  GREENS  employed  with  all  PV  arrays  and  HELBs,  from  [11]. 

smaller  load  requirements.  The  first  GREENS  unit  was  tested  in  July  2009  and  fielded  by 

2012  [12], 

Another  solar  power  generation  asset  developed  by  the  Marine  Corps  at  the  same  time  as 
the  GREENS  was  the  SPACES.  The  SPACES  is  mueh  smaller,  lighter,  and  less  powerful 
than  the  GREENS.  Powered  by  two  (12)  volt  folding  portable  solar  panels  a  single  SPACES 
is  designed  to  reeharge  two  BB-2590  lithium  batteries  or  power  battery  operated  deviees 
sueh  as  laptops  or  radios  [13].  Some  basie  speeifieations  for  the  SPACES  and  GREENS 
ean  be  found  in  the  Marine  Corps  Teehnieal  Manual  TM  12359A-OD:  Prineipal  Teehnieal 
Charaeteristies  of  Expeditionary  Power  Systems  Equipment;  more  detailed  speeifieations 
ean  be  requested  from  the  ONR  or  E20  [14]. 

1.3.3  Fielding  of  the  GREENS  and  SPACES  Systems 

The  rapid  development  and  fielding  of  the  GREENS  and  SPACES  demonstrated  the  Marine 
Corps’  drive  to  leverage  eutting  edge  renewable  energy  teehnology  on  the  battlefield.  One 
of  the  near-term  materiel  ineentives  speeified  in  the  United  States  Marine  Corps  (USMC) 
Expeditionary  Energy  Strategy  and  Implementation  Plan’s  Implementation  Planning  Guid- 
anee  is  the  initiative  to  "lead  in  deploying  innovative  energy  solutions  [3]."  This  initia¬ 
tive  represents  the  Marine  Corps  objeetive  of  not  only  developing  energy  teehnologies  but 
rapidly  fielding  them  to  the  warfighter.  The  GREENS  and  SPACES  systems  exemplify 
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Equipment 

Average  Hourly  Power 
Required  (W) 

Peak  Power 
Required  (W) 

Motorola  Battery  Charger 

50 

200 

VRC-110  w/Blue  Force  Tracker 

165 

440 

PRC-150 

57 

375 

Coffee  Pot 

45 

975 

_ 

Figure  1.2:  Examples  of  common  Marine  Corps  equipment  the  GREENS  system  can  power  with 
associated  power  requirements,  after  [11]. 

this  rapid  development  and  fielding  approach.  Fielding  began  in  2012  and  only  four  years 
after  the  Commandant  of  the  Marine  Corps  announced  the  Marine  Corps  Energy  Initia¬ 
tive,  Colonel  Roberto  Charette,  former  director  of  the  E20,  asserted  in  2013  that  both  the 
GREENS  and  SPACES  systems  were  "available  to  all  Marines  in  the  Marine  Corps  [15]." 

1.4  Changing  How  the  Marine  Corps  Decreases  Fuel  Con¬ 
sumption 

There  are  two  general  approaches  for  employing  renewable  energy  technologies  to  reduce 
fuel  consumption.  The  first  technique  includes  decreasing  the  peak  levels  of  a  load  by 
removing  components  of  that  load  and  powering  them  separately  with  renewable  energy 
technology  such  as  PV  arrays,  batteries  and  wind  turbines.  To  ensure  reliable  availability 
of  power,  this  will  usually  also  include  using  a  conventional  generator  as  a  backup  to  supply 
power  in  times  of  poor  weather.  The  Marine  Corps  currently  employs  the  GREENS  in  this 
way.  The  GREENS  controller,  which  limits  output  from  the  PV  arrays  to  IkW,  helps  to 
ensure  that  Marines  employing  the  system  do  not  over  extend  the  system  by  expecting  a 
consistent  provision  of  1.76  kW  from  the  arrays,  which  is  the  peak  amount  of  energy  the 
eight  PV  panels  are  capable  of  producing.  With  the  GREENS  limited  by  the  controller  to 
providing  1  kilowatt  (kW)  of  power,  the  system  is  capable  of  powering  common  Marine 
Corps  equipment  such  as  those  shown  in  Eigure  1.2. 

The  second  approach  for  employing  renewable  energy  to  decrease  fuel  consumption  in¬ 
cludes  employing  the  renewable  energy  technologies  to  decrease  the  operating  time  of 
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an  AC  generator(s)  covering  the  load.  This  eombines  the  generator  and  renewable  en¬ 
ergy  source  into  a  hybrid  system  that  passes  the  load  between  renewable  power  souree 
and  diesel  generator  depending  on  load  size  and  renewable  power  production.  In  many 
situations,  this  method  of  employing  a  renewable  energy  teehnology  can  do  more  to  de- 
erease  fuel  eonsumption  and  allow  for  higher  generator  operating  effieiency  than  the  first 
teehnique  described.  Recognizing  the  efficieney  of  this  method  of  renewable  energy  em¬ 
ployment  approach,  the  Marine  Corps  made  hybrid  systems  the  foeus  of  the  2013  EXFOB. 
This  investigation  of  the  potential  for  hybrids  was  referred  to  as  the  Mobile  Electrie  Hy¬ 
brid  Power  Sources  (MEHPS)  Analysis  of  Alternatives  (AOA)  and  sparked  the  pursuit  of 
MEHPS  technologies  for  the  Marine  Corps. 

1.4.1  EXFOB  2013:  MEHPS  AOA 

Following  the  Marine  Corps’  suceess  in  the  fielding  of  the  GREENS  and  SPACES  system, 
the  Corps  eontinued  to  push  forward  to  exploit  taetieal  and  logistie  opportunities  avail¬ 
able  in  renewable  energy  technology  development.  While  the  early  EXFOBs  succeeded  in 
developing  renewable  teehnologies  for  powering  small  loads  of  up  to  IkW,  with  EXFOB 
2013  the  Marine  Corps  began  its  pursuit  of  hybrid  systems  capable  of  powering  3kW  to 
300kW  loads  [16].  The  U.S.  Army  and  Marine  Corps  teamed  up  to  develop  the  hybrid 
systems  and  divided  the  task  of  developing  the  systems  between  the  two  serviees.  While 
the  Army  has  the  lead  in  development  of  the  10  to  300kW  systems,  referred  to  as  Mo¬ 
bile  Electrie  Microgrids  (MEM),  the  Marine  Corps  has  the  lead  in  the  development  of  the 
"smaller  hybrid  systems  (up  to  lOkW)  [17]."  It  should  be  noted  that  while  this  thesis  refers 
to  stand-alone  systems  composed  of  power  generation  assets  meeting  a  load  as  mieropower 
systems,  the  Marine  Corps  elassifies  these  systems  differently,  depending  on  the  size  of  the 
system’s  peak  load  eapability.  Systems  supporting  an  electrie  load  less  than  lOkW  are  elas- 
sified  within  the  Marine  Corps  as  hybrid  power  systems.  Systems  supporting  an  eleetrie 
load  greater  than  lOkW  are  elassified  as  microgrids.  For  this  thesis,  to  avoid  eonfusion,  all 
hybrid  power  systems  and  mierogrids  are  referred  to  simply  as  micropower  systems. 

These  hybrid  systems,  or  MEHPS,  will  provide  the  next  step  for  deereasing  fossil  fuel  eon- 
sumption  by  forees  in  an  expeditionary  environment.  MEHPS  hybrid  systems  also  have 
potential  for  deereasing  generator  maintenanee  issues  that  result  from  generators  operat¬ 
ing  with  too  low  of  power  loads.  One  such  maintenance  issue,  called  "wet-stacking,"  is 
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seen  when  generators  operate  with  loads  well  below  optimum  levels  and  results  in  "coke 
formation  and  gumming  of  the  engine  [18]."  To  meet  larger  loads  of  3kW  to  lOkW  while 
leveraging  renewable  energy,  the  MEHPS  systems  take  the  approach  of  passing  the  load 
between  the  renewable  power  source  and  diesel  generator  depending  on  the  load  size  and 
renewable  power  available  at  a  given  time.  This  passing  of  the  load  automatically  between 
power  sources  is  accomplished  through  the  use  of  smart  grid  technology.  By  combining 
"smart  controls,  energy  storage,  and  solar  PV  with  traditional  diesel  generators,"  these 
MEHPS  systems  can  significantly  decrease  generator  run  time,  preventing  generator  wet¬ 
stacking  and  decreasing  fuel  consumption  [17].  Eollowing  EXEOB  2013,  reductions  in 
generator  run  time  through  the  use  of  MEHPS  systems  were  seen  as  high  as  80  percent  and 
fuel  consumption  reductions  were  seen  as  high  50  percent  [17]. 

1.5  Marine  Corps  Employment  of  HOMER 

Beyond  the  renewable  technologies  and  hybrid  systems  themselves,  another  renewable  en¬ 
ergy  tool  the  Marine  Corps  has  pursued  as  part  of  its  energy  strategy  is  the  HOMER  model. 
While  several  expeditionary,  renewable  energy  assets  have  been  developed,  the  users  lack 
a  capability  for  predicting  the  fuel  consumption  rates  associated  with  the  renewable  en¬ 
ergy  hybrid  systems.  What  is  now  needed  to  decrease  fuel  consumption  is  a  tool  that  will 
assist  Marine  Corps  power  planners  and  logisticians  in  anticipating  the  fuel  consumption 
savings  to  be  gleaned  from  employing  hybrid  systems  to  support  a  given  load  requirement 
for  a  particular  location  and  time  frame.  Understanding  how  to  optimize  renewable  energy 
power  production  is  necessary  for  Marine  Corps  logisticians  and  engineers  to  optimally 
design  micropower  systems  with  hybrid  power  generation  technologies  and  modular  re¬ 
newable  technologies.  As  a  hybrid  optimization  model  for  electric  renewable  technologies, 
HOMER  has  the  potential  to  fulfill  this  need  for  the  Marine  Corps.  In  order  to  understand 
the  potential  benefits  and  uncertainties  associated  with  employing  the  HOMER  model  in 
support  of  expeditionary  operations,  the  inner  workings  of  the  HOMER  model  must  be 
understood.  Chapter  2  will  provide  insight  into  how  HOMER  models  micropower  systems 
and  the  performance  PV  arrays  as  well  as  the  model’s  origin  and  general  applications. 
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CHAPTER  2: 
HOMER 


2.1  HOMER  and  Its  Potential  Marine  Corps  Application 

HOMER  is  a  deterministic,  time-step  model  that  models  micropower  systems.  It  can  model 
systems  including  both  conventional,  diesel  powered  generators  and  renewable  power 
sources  such  as  photovoltaic  arrays  and  wind  turbines. 

HOMER  has  been  utilized  for  years  by  organizations  and  companies  across  the  globe  to 
assist  them  in  their  cost-effective  employment  of  renewable  energy  sources  and  appro¬ 
priately  sizing  generators  for  the  long-term  powering  of  facilities.  There  is  potential  for 
HOMER  to  serve  as  a  useful  Marine  Corps  Modeling  and  Simulation  tool  for  optimizing 
expeditionary  micropower  systems.  To  be  useful  for  modeling  expeditionary  micropower 
systems,  HOMER  must  provide  system  fuel  consumption  estimates  with  relatively  low 
variability  across  a  wide  range  of  locations,  micropower  system  compositions,  and  load 
profiles  for  short-term  operations. 

2.1.1  HOMER’S  Origin 

The  U.S.  National  Renewable  Energy  Eaboratory  (NREE)  developed  HOMER  in  1993 
as  a  part  of  its  Village  Power  Program.  The  model  was  designed  to  assist  users  in  de¬ 
signing  micropower  systems  and  modeling  their  performance  under  specific  conditions. 
Although  originally  built  only  for  consideration  of  off-the-grid  micropower  systems,  the 
model  was  updated  in  2000  to  accommodate  the  modeling  of  grid-tied  systems  as  well. 
The  HOMER  application  simulates  and  measures  the  performance  of  each  power  genera¬ 
tion  asset,  whether  available  individually  or  in  hybrid  configurations.  In  this  way,  the  model 
is  used  to  identify  the  most  cost-efficient  micropower  system  design  for  meeting  long  term 
electric  load  requirements. 

2.1.2  Discrete  Time  Simulations 

HOMER  is  a  type  of  model  known  as  a  Discrete  Time  Simulation  (DTS)  model,  meaning 
it  employs  a  time-advance  mechanism  (TAM)  driven  by  time-steps.  Understanding  TAMs, 
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and  DTS  models  in  particular,  will  assist  in  understanding  potential  strengths  and  weak¬ 
nesses  of  the  HOMER  model.  TAMs  are  methods  used  in  modeling  and  simulation  for 
representing  time  and  its  effects  in  the  system  being  simulated.  Three  primary  components 
employed  for  the  representation  of  time  advance:  the  states  of  variables  or  entities,  the 
events  that  change  those  states,  and  the  methods  for  timing  those  events.  By  using  these 
components  for  the  representation  of  time,  models  using  TAMs  can  represent  complex,  dy¬ 
namic  systems.  While  DTS  models  use  the  time-step  driven  TAM,  another  type  of  model 
known  as  the  Discrete  Event  Simulation  (DES)  uses  an  "event-driven"  TAM  method  [19]. 

Eor  the  timing  of  events,  DES  models  use  an  "event  list,"  where  future  scheduled  events 
are  populated  on  the  event  list  as  they  are  scheduled  by  state  changes  (events).  In  DTS,  or 
time-step,  simulations,  events  are  not  scheduled  by  using  events  to  populate  future  events 
on  event  lists.  Instead,  states  are  evaluated  on  a  recurring  basis  at  a  specified  time  intervals, 
with  events  depending  on  the  status  of  the  states  at  those  times.  The  "time  step"  is  the 
interval  of  time  between  evaluating  states  and  executing  appropriate  events.  While  both 
DTS  and  DES  models  are  valid  modeling  and  simulation  tools,  modeling  and  simulation 
professionals  must  be  careful  in  determining  which  technique  to  employ  and  in  what  way. 
There  are  potential  problems  with  both  TAM  methods  that  must  be  considered  to  prevent 
or  mitigate  particularly  detrimental  effects. 

An  example  of  a  well-known  issue  seen  with  the  time-step  method  is  the  "bullet  through 
paper"  problem,  an  issue  referred  to  as  "skipping"  in  the  Naval  Postgraduate  School  (NPS) 
thesis  "The  Effect  of  Time  Advance  Mechanism  in  Modeling  and  Simulation  [19]."  As  the 
name  implies,  this  issue  refers  to  the  problems  associated  with  modeling  a  bullet  hitting  a 
piece  of  paper  using  a  time-step  model.  The  "impact"  of  the  bullet  hitting  the  paper  can 
easily  be  modeled  by  a  non-time  step  DES  model  because  the  event  list  assigns  the  event  to 
occur  at  that  specific  time  calculated  for  impact.  When  representing  this  event  through  time 
step,  however,  the  model  will  represent  two  states,  the  state  with  the  bullet  on  one  side  of 
the  paper  and  a  following  state  with  the  bullet  on  the  other  side  of  the  paper,  not  represent 
the  impact.  This  example  may  be  confused  as  applicable  only  to  physics  models,  but  this 
"bullet  through  paper"  problem  can  apply  to  non-physics  modeling  tasks  as  well. 

Another  concern  with  employing  the  time-step  technique  is  that  there  is  no  consideration 
for  the  ordering  of  events  for  those  events  occurring  within  the  same  time  step.  One  way 
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that  time-step  models  are  often  evaluated,  and  a  way  in  whieh  this  partieular  eoneem  eould 
be  addressed,  is  the  analysis  of  the  simulation  time-step  interval  size.  Researehers  vary 
the  size  of  the  time-step  and  analyze  how  those  ehanges  affeet  the  output  of  the  model. 
The  time  step  seleeted  for  a  model  is  an  important  eonsideration  beeause  "too  large  a  time 
step  introduees  errors  in  the  model,  while  too  small  a  time  step  may  eause  extremely  large 
run  times  [20]."  The  HOMER  model  graphieal  user  interface  (GUI)  allows  users  to  change 
the  size  of  the  time-step  in  terms  of  minutes.  As  will  be  seen  in  later  sections  of  this 
thesis,  limitations  on  the  fidelity  of  solar  irradiance  data  and  input  profiles  makes  time 
steps  of  less  than  an  hour  impractical  with  regard  to  reflecting  changes  in  solar  irradiance 
and  temperature  values.  That  being  said,  future  analysis  of  the  HOMER  model  could  be 
useful  in  identifying  how  changes  to  the  time-step  size  affect  simulations’  representation 
of  generator  operating  times  and  the  duration  battery  charging/discharging. 

Statistician  George  Box  once  famously  stated  that  "all  models  are  wrong,  but  some  are 
useful  [21]."  It  is  the  responsibility  of  the  modeling  and  simulation  professionals  to  un¬ 
derstand  how  the  models  they  employ  are  wrong.  This  leads  to  considerations  of  model 
classifications  as  being  either  stochastic  or  deterministic.  This  is  of  particular  importance 
when  considering  the  HOMER  model  because  although  the  HOMER  model  evaluates  solar 
power  on  the  basis  of  stochastic  solar  irradiance  inputs,  the  model  itself  is  deterministic. 
This  is  a  third  area  requiring  greater  research,  and  is  a  primary  focus  of  this  research. 

2.2  How  HOMER  Models  Micropower  Systems 

The  HOMER  engine  is  the  component  of  the  software  that  models  the  performance  of  a 
given  micropower  system  with  a  particular  scenario  including  a  set  of  inputs,  such  as  load, 
temperature,  and  solar  irradiance  profiles.  The  model  engine  is  separate  from  the  optimiza¬ 
tion  software  which  iterates  through  the  different  power  generation  assets  and  identifies  all 
possible  combinations  for  analysis. 

2.2.1  Micropower  Systems  and  Load  Profiles 

In  many  situations  a  “micropower  system”  refers  to  a  power  generator,  such  as  a  diesel 
generator,  that  is  smaller  than  the  average  scale.  In  the  context  of  the  HOMER  model 
and  this  thesis,  the  term  “micropower”  refers  to  the  provision  of  power  to  electric  devices 
using  local  power  generation  assets,  that  is  to  say,  off  the  grid.  A  micropower  system. 
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therefore,  refers  to  the  system  of  power  generation  devices  (solar  arrays,  generators,  wind 
turbines,  etc.),  energy  storage  apparatuses  (batteries),  and  connective  wiring  that  provides 
the  required  energy  for  that  load  in  the  form  of  electricity. 

A  load  profile  is  a  representation  of  the  energy  consumed  by  a  piece  or  pieces  of  equipment 
in  specified  increments  of  time  over  a  specified  amount  of  time.  In  the  HOMER  model,  load 
profiles  of  micropower  systems  are  input  by  the  user  in  one-hour  increments,  matching  the 
one-hour  time  steps  of  the  discrete  event  simulation.  Each  value  in  a  load  profile  is  the  sum 
of  energy,  in  kW  that  is  drawn  from  the  power  sources  (generators,  batteries,  PV  arrays, 
etc.)  during  the  specified  hour. 

Eoad  profiles  can  be  loaded  in  one  of  two  ways  in  HOMER.  Eirst,  the  user  can  specify  a 
generic  one-day  load  profile  for  weekdays  and  weekend  days  of  a  given  month.  In  using  this 
option  of  load  profile  specification,  the  user  can  either  keep  weekday  and  weekend  daily 
load  profiles  the  same,  or  can  account  for  some  difference  in  micropower  system  operation 
over  the  weekend  by  making  the  daily  load  profiles  unique.  When  the  user  specifies  the 
weekday/weekend  daily  load  profiles  for  each  month  the  HOMER  GUI  assigns  that  daily 
load  profile  to  every  weekday/weekend  day  in  the  month.  The  second  option,  and  the 
option  used  to  support  this  research,  is  loading  a  single  Comma  Separated  Value  (CSV)  file 
consisting  of  8760  values.  These  8760  values  provide  the  micropower  system’s  load,  in 
kW,  for  every  hour  of  a  year. 

If  the  HOMER  model  is  run  in  a  retrospective  study  of  a  system’s  performance,  the  actual 
load  profile  of  the  micropower  system  can  be  measured  during  the  experiment  and  then 
loaded  into  the  HOMER  model  in  the  form  of  a  CSV  file.  Often  this  is  not  the  case  though, 
and  HOMER  is  employed  to  analyze  expected  future  micropower  system  performance.  To 
conduct  this  analysis,  the  user  must  provide  the  model  with  a  load  profile  that  represents 
anticipated  future  load  requirements  of  the  micropower  system  in  question.  This  research 
falls  into  the  second  category,  HOMER’s  employment  to  analyze  expectations  of  future 
performance  across  multiple  factors.  Because  of  this,  the  load  profiles  used  in  this  research 
are  not  measured  values  but  are  instead  extrapolations  based  on  historical  data.  To  analyze 
HOMER’s  capabilities  for  to  the  Marine  Corps’  needs,  I  utilized  a  modification  of  the  load 
profile  design  used  in  the  E20’s  MEHPS  AOA  [16].  This  load  profile  design  is  meant  to 
represent  the  load  profile  of  a  typical  Marine  Corps  unit  of  a  specified  peak  power  capacity. 


14 


Figure  2.1:  30kW  modified  seasonal  profile 


Figure  2.2:  30kW  standard  seasonal  profile 


Figure  2.3:  30kW  standard  profile  data  map 


Figure  2.4:  30kW  modified  profile  data  map 


The  profile  used  in  the  MEHPS  AOA,  shown  in  Figure  2.1,  splits  the  year  up  into  four  sea¬ 
sons  and  assigns  a  set  load  profile  for  eaeh  day  of  the  season.  The  load  profile  is  speeified  as 
a  pereentage  of  peak  power.  The  above  visualization  of  a  standard  load  profile,  Figure  2.1, 
was  ereated  using  the  HOMER  GUI  and  represents  the  projeeted  annual  load  profile  of  a 
30kW  peak  power.  This  load  profile,  utilized  by  the  Marine  Corps  for  the  MEHPS  AOA, 
alloeates  the  same  load  profile  values  to  every  day  of  the  season,  with  an  immediate  ehange 
in  the  load  profile  values  upon  reaehing  the  next  season.  These  ehanges  in  load  profiles 
show  how  power  demands  ehange  for  different  seasons  throughout  the  year,  but  they  show 
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the  shift  in  load  profiles  as  oeeurring  instantaneously.  The  winter  and  spring  load  profiles, 
for  example,  suggest  that  the  day  the  season  shifts  from  winter  to  spring  the  first  hour  of  the 
load  profile  will  shift  from  a  load  of  approximately  82  percent  to  a  load  of  5  percent.  This 
is  not  a  realistic  representation  of  how  load  profiles  change  throughout  the  year;  season 
changes  and  their  associated  effects  on  load  profiles  are  not  immediate,  but  gradual. 

For  this  research,  the  MEHPS  load  profile  product  was  modified  to  reflect  the  gradual 
nature  with  which  seasons  change.  This  modified  load  profile  was  created  by  interpolating 
daily  load  profiles  between  the  seasons.  Season  midpoints  were  set  as  the  dates  equidistant 
between  the  summer  and  winter  solstices,  and  the  autumnal  and  vernal  equinoxes  for  the 
respective  seasons  [22].  For  instance,  load  profiles  constructed  throughout  the  month  of 
April,  which  lies  between  the  vernal  equinox  and  summer  solstice,  were  constructed  by 
combining  varying  percentages  of  spring  and  summer  load  profile  levels.  These  four  season 
midpoints  were  set  to  have  load  profiles  equal  to  the  MEHPS  load  profiles.  The  load  profiles 
for  the  rest  of  the  361  days  of  the  year  were  calculated  by  combining  the  load  profiles  for  the 
season  midpoints  before  and  after  the  selected  date.  The  season  midpoints  were  combined 
with  weights  relative  to  the  given  date’s  distance  from  each  season  midpoint  date.  The 
modified  load  profile  is  shown  in  Eigure  2.2.  While  the  two  seasonal  profiles,  Eigure  2.1 
and  Eigure  2.2,  show  very  similar  annual  profiles  by  month,  the  Data  Maps,  Eigure  2.3 
and  Eigure  2.4,  demonstrate  how  blending  the  load  profile  values  provides  a  more  realistic 
representation  of  a  load  profile’s  progression  through  the  seasons. 

2.2.2  Defining  Micropower  Component  Specifications 

Modeling  of  micropower  systems  with  the  HOMER  model  includes  not  only  quantifying 
the  load  profile,  but  also  identifying  the  number  and  characteristics  of  power  generation 
and  storage  equipment  that  comprise  the  system.  These  pieces  of  equipment  can  include 
internal  combustion  generators,  battery  banks,  power  converters,  PV  arrays,  and  other  types 
of  renewable  energy  equipment.  The  technical  specifications  for  each  piece  of  equipment 
are  either  pre-loaded  into  the  HOMER  GUI,  if  it  is  a  common  piece  of  equipment,  or  are 
specified  by  the  user  building  the  scenario  for  simulation.  The  technical  specifications 
allow  the  HOMER  model  to  determine  the  amount  of  power  different  pieces  of  equipment 
can  generate,  store,  or  convert  for  each  time-step  of  the  simulation.  Specifications  for  the 
internal  combustion  generator  include  generator  size  (kW),  minimum  load  ratio  (percent). 
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and  whether  it  is  an  AC  or  DC  generator.  Generator  specifications  also  include  the  fuel  type 
utilized  and  values  necessary  for  demonstrating  the  efficiency  curve  of  the  fuel,  including 
the  fuel  curve  intercept  coefficient  and  slope. 

Modeling  a  PV  array  in  the  HOMER  model  requires  specifying  the  size  of  the  array  (kW), 
the  array  output  current  (AC  or  DC),  and  the  expected  lifetime  of  the  array.  The  HOMER 
model  also  requires  users  to  specify  properties  related  to  the  employment  of  the  PV  array. 
These  include  a  derating  factor  (percent)  associated  with  the  panel,  the  slope  (degrees)  and 
azimuth  with  which  the  PV  array  is  employed.  A  solar  panel’s  derating  value  is  defined  by 
HOMER  as  "a  scaling  factor  meant  to  account  for  effects  of  dust  on  the  panel,  wire  losses, 
elevated  temperature,  or  anything  else  that  would  cause  the  output  of  the  PV  array  to  deviate 
from  that  expected  under  ideal  conditions  [1]."  If  the  user  decides  to  include  the  effects  of 
temperature  on  the  PV  array  in  the  model,  PV  array  technical  specifications  also  include  the 
temperature  coefficient  of  power,  NOCT,  and  efficiency  at  standard  test  conditions  for  the 
array.  The  importance  of  each  of  these  values  and  how  they  assist  the  model  in  forecasting 
PV  energy  production  will  be  discussed  in  the  next  section  as  the  inner  workings  of  PV 
cells  are  explained.  Eurther  information  on  the  modeling  of  micropower  components  in  the 
HOMER  model  and  their  technical  specifications  can  be  found  in  chapter  15  of  the  book 
Integration  of  Alternative  Sources  of  Energy,  titled  "Micropower  System  Modeling  with 
HOMER  [1]." 

2.3  Modeling  Solar  Energy 

2.3.1  Solar  Energy  Modeling  in  the  Solar  Industry 

Anticipating  the  performance  of  solar  energy  systems  requires  an  ability  to  anticipate  the 
solar  resource  available.  Eor  years  the  solar  industry  has  relied  upon  NREE’s  Typical  Mete¬ 
orological  Year  (TMY)  data  set  for  modeling  the  performance  of  solar  power  systems  [23]. 
A  TMY  is  a  statistically  derived  profile  of  a  typical  meteorological  year  for  certain  loca¬ 
tions  across  the  United  States.  NREE’s  TMY  data  consists  of  hourly  measurements  of  solar 
irradiance  and  certain  meteorological  elements  and  has  been  updated  three  times  since  its 
establishment  in  1978  by  Sandia  National  Eaboratories  [24].  The  current  TMY,  TMY3,  was 
created  with  measurements  acquired  from  1976  to  2005.  TMY3  data  is  meant  to  be  used 
for  simulation  of  solar  power  systems  over  "a  longer  period  of  time,  such  as  30  years  [24]." 
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In  fact,  the  first  page  of  the  TMY3  User  Manual  warns  in  italics:  "The  TMY  should  not  be 
used  to  predict  weather  for  a  particular  period  of  time  [24]." 

If  the  user  inputs  solar  irradiance  profiles  in  terms  of  average  monthly  irradiance  values, 
the  HOMER  model  uses  the  Graham-Hollands  algorithm  to  produce  a  synthetic  hourly 
solar  irradiance  profile  based  on  the  average  solar  irradiance  value.  The  Graham-Hollands 
algorithm,  developed  in  1990,  allows  HOMER  to  create  a  synthetic  solar  irradiance  profile 
that  estimates  the  "likely  magnitudes  of  random  fluctuations"  in  solar  irradiance  [25].  This 
algorithm  allows  the  HOMER  model  to  be  used  even  when  the  user  cannot  obtain  a  TMY 
solar  irradiance  profile  or  measured  hourly  solar  irradiance  profile. 

2.3.2  The  Basics  of  PV  Cells 

To  understand  how  PV  arrays  produce  power,  and  the  effects  of  solar  irradiance  and  tem¬ 
perature  on  them,  one  must  first  understand  some  basics  of  semiconductors,  the  process  of 
doping,  and  the  structure  of  PV  cells.  The  following  descriptions  are  abridged  explanations 
providing  only  so  much  insight  as  is  necessary  for  the  purposes  of  this  thesis. 

PV  cells  can  be  made  of  either  specific  elements,  such  as  silicon,  or  compounds,  but  the 
material  must  be  a  semiconductor.  Semiconductors  are  materials  with  lower  electrical  con¬ 
ductivity  than  conductors  and  higher  conductivity  than  insulators.  Conductivity  is  deter¬ 
mined  by  the  material’s  number  of  free  electrons  and  its  number  of  holes,  where  a  hole  is 
simply  "the  absence  of  an  electron  [26]."  A  high  level  of  conductivity  means  having  a  high 
number  of  free  electrons  or  a  high  number  of  holes.  A  pure  semiconductor  has  a  relatively 
low  conductivity  due  to  the  absence  of  naturally  occurring  free  electrons  or  holes. 

Tree  electrons  are  established  in  a  semiconductor,  creating  a  potential  for  current  flow, 
by  applying  energy,  which  breaks  covalent  bonds  between  the  semiconductor’s  atoms  or 
molecules.  The  amount  of  energy  required  to  achieve  this  break  is  known  as  the  band  gap 
energy.  Breaking  the  covalent  bonds  also  creates  a  hole  for  every  electron  that  is  freed. 
While  the  breaking  of  covalent  bonds  must  occur  for  increasing  potential  current  flow, 
semiconductors  must  also  undergo  a  process  called  doping  to  be  used  in  PV  cells. 

Doping  is  the  application  of  an  element  or  compound  to  a  PV  cell’s  semiconductor  that 
increases  either  the  number  of  holes  or  the  number  of  free  electrons  in  the  semiconductor. 


18 


To  do  this,  the  applied  material  must  have  either  greater  or  fewer  electrons  in  its  valence 
level  than  does  the  PV  cell’s  semiconductor.  Semiconductors  doped  with  a  material  having 
fewer  electrons  in  its  valence  level  than  the  semiconductor  have  an  increased  number  of 
holes  and  are  known  as  p-type  semiconductors.  Similarly,  those  semiconductors  doped 
with  materials  having  more  electrons  in  its  valence  level  than  the  semiconductor  have  a 
higher  number  of  free  electrons  and  are  known  as  n-type  semiconductors.  This  distinction 
is  particularly  important  for  understanding  the  construct  and  operation  of  PV  cells. 

Having  identified  the  construct  of  semiconductors,  the  concept  of  how  PV  cells  generate 
electricity  is  relatively  easy  to  explain.  When  sunlight  strikes  a  PV  cell,  power  is  gener¬ 
ated  through  photogeneration.  Brian  Sullivan  explains  how  the  process  of  photogeneration 
works  in  his  NPS  thesis  addressing  temperature  effects  on  PV  cells: 

Photons  from  the  source  of  light  enter  the  semiconductor  material  when  they 
contain  energy  equal  to  or  greater  than  the  band  gap  of  the  material  and  trans¬ 
fer  enough  energy  to  electrons  in  the  valence  band  to  cause  them  to  move  to 
the  conduction  band.  This  transfer  of  energy,  therefore,  creates  electron-hole 
pairs. ..as  electron  pairs  appear,  the  electric  field  of  the  depletion  region  causes 
the  carriers  to  separate  to  opposite  ends  of  the  material... The  flow  of  carriers 
to  opposite  ends  of  the  device  creates  a  current  in  the  material  and  a  voltage 
differential  on  the  two  ends  [26]. 

Greater  detail  on  the  construct  of  PV  cells  and  the  semiconductors  that  they  are  composed 
of  can  be  found  in  the  2010  thesis,  "The  Effect  of  Temperature  on  the  Optimization  of 
Photovoltaic  Cells  Using  Silvaco  ATLAS  Modeling,"  by  Sullivan  [26].  Sullivan’s  thesis 
also  provides  insight  into  the  effects  of  ambient  temperature  levels  on  the  performance  of 
PV  cells,  a  relationship  which  will  now  be  discussed  in  this  thesis  as  well. 

2.3.3  Temperature  Effects  on  PV  Cell  Performance 

Changes  in  a  PV  cell’s  operating  temperature  affect  the  power  production  of  the  PV  cell, 
with  increases  in  PV  cell  temperature  resulting  in  a  net  decrease  in  cell  power  production. 
This  relationship  can  be  explained  in  terms  of  current  and  voltage  changes  relative  to  cell 
temperature,  with  temperature  increases  resulting  in  a  slight  increase  in  current  and  a  much 
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larger  relative  deerease  in  voltage.  In  order  to  understand  the  HOMER  model’s  representa¬ 
tion  of  temperature  effeets,  this  relationship  between  temperature  and  PV  eell  performance 
will  be  explained  in  this  section. 

Between  current  and  voltage  effects,  the  less  significant  effect  of  temperature  change  is  the 
impact  on  current.  Several  different  elements,  or  compounds,  can  be  used  for  constructing 
PV  cells,  and  these  different  materials  have  different  band  gap  energy  levels.  Increases  in 
operating  cell  temperatures  increase  the  thermal  energy  provided  to  the  semiconductor  and 
decrease  the  amount  of  radiant  energy  required  to  reach  the  band  gap  level.  The  addition  of 
this  thermal  energy  due  to  increased  cell  temperature  is  what  increases  short-circuit  current 
flow  in  the  PV  cell,  given  a  constant  solar  irradiance,  "but  only  slightly  [27]."  The  increase 
in  current  is  minimal  when  compared  to  the  effects  temperature  has  on  PV  cell  open-circuit 
voltage  [27]. 

While  increases  in  cell  operating  temperature  slightly  increase  the  cell’s  short-circuit  cur¬ 
rent,  increases  in  cell  temperature  also  increase  the  resistivity  of  components  within  the 
PV  array  and  thereby  decrease  the  voltage  produced  by  the  PV  cell  [26].  The  decrease 
in  PV  cell  voltage  due  to  temperature  increases  is  much  greater  than  the  relative  increase 
in  current.  This  disparity  between  temperature  effects  on  current  and  voltage  results  in  a 
negative  linear  relationship  between  PV  cell  operating  temperature  and  PV  cell  power  pro¬ 
duction.  Figure  2.5  provides  a  graphical  representation  of  the  relationship  between  PV  cell 
temperature,  current  and  voltage  [28]. 

The  relationship  between  temperature  and  PV  cell  power  production  is  quantified  through 
the  specification  of  a  nominal  operating  cell  temperature  (NOCT)  and  a  temperature  coef¬ 
ficient  of  power  for  every  PV  cell.  The  NOCT  specifies  the  "temperature  of  [the  PV  cell] 
at  the  conditions  of  the  nominal  terrestrial  environment  (NTE),"  which  includes  800  VP/m^ 
for  solar  irradiance,  an  ambient  temperature  of  20  degrees  Celsius,  and  an  average  wind 
speed  of  1  m/5  [29].  The  temperature  coefficient  of  power  specifies  "how  strongly  the  PV 
array  power  output  depends  on  the  cell  temperature"  and  is  a  negative  value  due  to  the 
negative  linear  relationship  between  operating  temperature  and  power  production  for  PV 
cells  [30]. 

The  HOMER  model  simulates  the  effects  of  temperature  on  a  PV  cell’s  performance  at 
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Figure  2.5:  Representation  of  temperature  effects  on  current  and  voltage  in  a  PV  cell 


PV  Cell  Power  Production: 

where: 


"  ^  pvf  ?l> 


} py  is  the  rated  capacity  of  the  PV  array,  meaning  its  power  output  under  standard  test  conditions  [kW] 

fpy  is  the  PV  derating  factor  [%] 

Gj  is  the  solar  radiation  incident  on  the  PV  array  in  the  current  time  step  [kW/m^] 

Gy is  the  incident  radiabon  at  standard  test  conditions  [1  kw/m^] 

^p  is  the  temperature  coefficient  of  power  [%/^C] 

is  the  PV  cell  temperature  in  the  current  bme  step  [^C] 

^c57C  temperature  under  standard  test  conditions  [25  "C] 


Figure  2.6:  HOMER  model  equation  for  calculating  PV  cell  power  output 

each  time  step  of  the  simulation.  The  model  calculates  PV  cell  temperature  and  models 
the  effects  of  this  temperature  on  the  cell’s  power  production  for  each  time  step.  The 
model  uses  user-specified  NOCT  and  temperature  coefficient  of  power  specifications  for  the 
given  PV  cell  to  calculate  the  temperature  effects  on  power  production.  The  HOMER  GUI 
identifies  Figure  2.6  as  the  equation  used  by  the  model  to  determine  the  power  production 
for  a  given  PV  cell. 

The  PV  cell  temperature  in  this  equation  includes  the  temperature  effects  of  both  ambient 
temperature  and  the  increased  temperature  derived  from  solar  radiation  striking  the  cell. 
Figure  2.7  shows  the  HOMER  GUI-provided  equation  showing  how  the  HOMER  model 
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PV  Cell  Temperature:  j- =7- +g  I  — 

^  rc^J 

where: 

t  is  the  solar  transmittance  of  any  cover  over  the  PV  array  [%] 
a  is  the  solar  absorptance  of  the  PV  array  [%] 

^7  is  the  solar  radiation  striking  the  PV  array  [kW/m^] 

is  the  electrical  conversion  efficiency  of  the  PV  array  [%] 
is  the  coefficient  of  heat  transfer  to  the  surroundings  [kW/fn^®C] 
T'g  is  the  PV  cell  temperature  [®C] 
is  the  ambient  temperature  [®C] 


Figure  2.7:  HOMER  model  equation  for  calculating  PV  cell  temperature 
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Figure  2.8:  PV  percent  power  production  change  for  a  PV  cell  with  a  temperature  coefFicient  of 
power  of  -0.336 


calculates  the  PV  eell  temperature  based  on  ambient  temperature  and  the  temperature  ef- 
feets  of  solar  irradiance.  When  eonsidering  the  effeets  of  ambient  temperature  changes 
alone,  a  PV  eell’s  power  produetion  percent  change  that  results  from  ambient  temperature 
ehange  ean  be  determined  by  multiplying  the  temperature  delta  by  the  temperature  eoeffi- 
eient  of  power.  An  example  of  how  ambient  temperature  affeets  pereent  power  produetion 
ean  be  seen  in  Figure  2.8.  This  figure  shows  how  a  PV  eell  with  a  temperature  eoefficient 
of  power  of  -0.00336  has  a  PV  power  produetion  deerease  of  .336  pereent  for  every  degree 
Celsius  that  ambient  temperature  inereases  and  an  increase  in  power  of  0.336  percent  for 
every  degree  Celsius  that  ambient  temperature  deereases. 
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2.4  HOMER  Model  Resource  Inputs 

2.4.1  Solar  Resource  Input 

"Solar  resource"  refers  to  the  amount  of  global  horizontal  solar  irradiance  (GHI)  that  the 
specified  location  receives,  per  hour,  over  the  operation’s  time  period.  GHI  is  defined  by 
NREL  as  the  sum  of  Direct  Normal  Irradiance  (DNI)  and  Diffuse  Horizontal  Irradiance 
(DHI)  for  all  practical  purposes.  DNI,  which  will  be  discussed  later,  is  defined  as  "the 
amount  of  solar  radiation  from  the  direction  of  the  sun  [31]."  Similar  to  a  load  profile,  the 
solar  resource  can  be  loaded  by  the  user  in  one  of  two  ways.  First,  the  user  can  provide 
an  average  daily  GHI  value  for  each  month.  The  second  method,  and  the  method  used  in 
this  research,  is  to  load  an  8760- value  CSV  file  containing  a  GHI  value  for  every  hour  of 
the  year.  Also  similar  to  the  load  profile,  the  user  can  either  input  historically  measured 
solar  irradiance  data,  or  estimated  solar  irradiance  data.  A  common  practice  while  using 
HOMER  is  to  use  the  NREL’s  TMY3  data  as  the  input.  When  the  user  cannot  provide 
solar  irradiance  data  for  the  given  location  and  time  frame,  the  HOMER  GUI  obtains  the 
required  solar  irradiance  values  from  either  NASA  or  NREE  databases. 

2.4.2  Temperature  Profile  Input 

In  order  for  the  HOMER  model  to  simulate  the  effects  of  temperature  on  PV  cell  perfor¬ 
mance,  temperature  data  must  be  provided  to  the  model  for  the  location  and  time  frame 
of  the  simulation.  Once  the  user  has  selected  to  include  temperature  effects  in  the  simu¬ 
lation,  temperature  input  data  is  loaded  by  selecting  the  "Temperature"  resource  button  on 
the  HOMER  GUI  [30].  The  temperature  input  data  can  be  provided  in  the  form  of  either 
an  annual  temperature  profile  with  unique  hourly  temperature  values  or  monthly  average 
temperature  values.  If  the  user  enters  monthly  temperature  averages,  the  model  uses  those 
average  values  for  every  time  step  of  the  simulation  for  the  respective  month.  If  the  user 
imports  an  hourly  time  series  data  file,  in  an  8760- value  CSV  file  form,  then  the  HOMER 
model  uses  each  of  the  unique  hourly  values  for  the  associated  hour  of  the  time  step  simu¬ 
lation. 
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CHAPTER  3: 

HOMER  in  Expeditionary  Operations 


3.1  HOMER  Robustness  for  Expeditionary  Operations 

The  Marine  Corps  boasts  an  ability  to  deploy  to  any  clime  or  place.  Due  to  the  unpre¬ 
dictable  nature  of  combat  and  humanitarian  operations,  the  Marine  Corps  does  not  have 
the  luxury  of  being  able  to  gather  data  on  the  solar  irradiance  and  weather  patterns  of  a 
location  for  months  (or  years)  prior  to  deploying  there.  To  be  useful  for  deployment  to 
unanticipated  locations,  HOMER  must  be  robust  enough  to  handle  the  solar  irradiance  and 
temperature  variability  for  those  locations.  Of  particular  concern  when  considering  the  use 
of  HOMER  for  Marine  Corps  power  planning  is  how  HOMER  simulates  micropower  sys¬ 
tem  performance  for  the  short-term  durations  associated  with  Marine  Corps  expeditionary 
operations.  HOMER  was  designed  for  "simulating  the  long-term  operation  of  a  microp¬ 
ower  system  [1]."  While  modifications  will  be  required  to  be  made  to  HOMER  to  simplify 
the  simulation  of  shorter  time  periods,  this  modification  is  expected  to  impact  the  GUI  and 
not  the  model  engine  itself.  This  research  assists  in  determining  whether  the  HOMER  en¬ 
gine  is  compatible  with  the  simulation  of  short  duration  micropower  simulations  or  if  the 
inter-annual  variability  of  environmental  inputs  induce  unacceptable  variability  in  HOMER 
simulation  output. 

Eor  this  research,  it  is  assumed  that  when  employing  HOMER  the  Marine  Corps  will  de¬ 
pend  on  the  historical  solar  irradiance  data  gathered  by  NASA  and  NREE.  This  is  a  standard 
practice  in  the  solar  industry  for  performance  simulation  and  is  a  preferred  option,  when 
possible,  for  using  the  HOMER  model.  A  typical  solar  irradiance  profile  is  utilized  as  the 
solar  resource  input  for  HOMER  because  the  model  was  originally  designed  to  analyze 
the  long-term  performance  of  micropower  systems,  specifically  for  30  years.  This  solar 
irradiance  profile  can  be  derived  either  from  the  NREE  TMY  database  or  built  using  the 
Graham- Hollands  algorithm  and  a  value  for  the  location’s  average  monthly  solar  irradiance. 
Over  the  long  term  of  30  years,  the  effects  of  interannual  solar  irradiance  and  temperature 
variability  are  expected  to  be  negligible.  This  assumption  is  based  on  the  predicted  pro¬ 
files’  errors  above  and  below  the  actual  measured  solar  irradiance  and  temperature  profiles 


25 


of  each  year  averaging  out  over  the  thirty  years.  With  the  Marine  Corps  likely  employing 
the  HOMER  model  for  planning  power  generation  in  operations  as  short  as  15  days,  this 
variability  can  no  longer  be  assumed  to  have  a  negligible  impact  and  must  be  investigated. 

In  2010,  Newell  conducted  an  experiment  where  he  assessed  the  capabilities  of  the 
HOMER  model  in  forecasting  the  power  output  of  a  solar  panel  at  NFS  [32].  In  his  ex¬ 
periment  Newell  utilized  the  Graham-Hollands  algorithm  and  average  monthly  solar  irra- 
diance  measurements  from  the  NASA  Surface  Solar  Energy  data  set  to  build  a  synthetic 
solar  profile  which  served  as  the  solar  resource  input  for  the  HOMER  model.  He  also  uti¬ 
lized  average  monthly  temperature  measurements  as  the  temperature  input  for  running  the 
model.  Newell’s  research  initially  found  the  HOMER  model  to  perform  with  unacceptably 
poor  fidelity,  with  HOMER  estimating  "an  energy  production  level  that  was  over  25  per¬ 
cent  higher  than  the  actual  measured  energy  [32]."  Newell  associated  this  overestimate  of 
power  production  to  three  factors:  a  disparity  between  measured  and  anticipated  surface 
solar  irradiance,  a  disparity  between  measured  and  anticipated  ambient  temperature,  and 
an  underestimate  of  the  derating  value  for  the  solar  panels. 

As  discussed  earlier,  in  order  for  the  HOMER  model  to  simulate  the  performance  of  the 
solar  panels,  the  user  must  load  the  solar  resource  and  temperature  profile  into  the  model. 
The  25  percent  overestimate  in  energy  production,  identified  in  Newell’s  thesis  work,  was 
predicted  with  resource  inputs  based  on  average  monthly  solar  irradiance  and  temperature 
values  and  the  generation  of  a  solar  irradiance  profile  with  the  Graham-Hollands  algorithm. 
To  identify  the  impact  of  discrepancies  between  the  inputted  and  the  real  temperature  and 
solar  irradiance  profiles,  Newell  reran  the  model  with  the  measured  temperature  and  solar 
irradiance  profiles  for  that  period  of  time.  Eigure  3. 1  shows  the  results  of  Newell’s  analysis. 
These  HOMER  runs  showed  that  the  differences  between  real  world  and  originally  inputted 
temperature  and  solar  irradiance  profiles  accounted  for  10  percent  of  the  HOMER  model’s 
25  percent  overestimate  in  PV  power  production  for  Newell’s  experiment.  Six  percent 
of  the  overestimate  in  power  production  was  due  to  differences  between  real  and  original 
temperature  profiles.  Eour  percent  of  the  overestimate  in  PV  power  production  was  due  to 
differences  between  real  and  Graham-Hollands-generated  solar  irradiance  profiles. 

Newell  concluded  that  the  model  performed  satisfactorily  if  the  derating  factor  for  the  array 
was  assumed  to  have  decreased  below  the  standard  80  percent,  based  on  the  weathered  and 
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Figure  3.1:  Effects  of  solar  irradiance  and  temperature  discrepancies  in  Major  Brandon  Newell's 
2010  solar  panel  and  HOMER  model  experiment,  from  [32]. 

unclean  status  of  the  partieular  array  [32].  While  the  effeets  of  solar  irradianee  and  temper¬ 
ature  diserepaneies  were  elearly  presented  in  Newell’s  thesis,  the  question  of  identifying 
the  eorreet  derating  faetor  is  less  justified.  Without  the  support  of  follow-on  experiments 
to  validate  the  derating  faetor  Newell  applied,  the  eonelusion  of  subseribing  the  remaining 
17  pereent  inaeeuraey  to  derating  of  the  solar  panel  is  dangerously  elose  to  assuming  the 
eonelusion.  That  being  said,  this  researeh,  whieh  foeuses  on  the  effeets  of  input  variability, 
assumes  Newell’s  assessment  of  the  HOMER  model’s  aeeuraey  to  be  eorreet.  In  order  for 
this  researeh  to  be  better  grounded  future  researeh  would  need  to  be  eondueted  to  validate 
the  derating  faetor  laid  out  in  Newell’s  assessment  of  the  HOMER  model. 


3.2  Solar  Irradiance  Variability 

Marine  Corps  expeditionary  operations  depend  on  robust,  spaee  and  weight  effieient 
sourees  of  power.  An  important  eomponent  of  expeditionary  operations  is  the  relatively 
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small  size  of  the  units  conducting  fast  paced  and  short  duration  operations.  Logistics  plan¬ 
ners  and  commanders  must  be  able  to  rely  on  their  energy  source,  and  solar  irradiance 
variability  is  one  element  of  solar  irradiance  that  brings  that  reliability  into  question.  Solar 
irradiance  variability  can  be  caused  by  multiple  factors,  including  varying  cloud  cover  and 
pollution. 

Newell’s  2010  thesis  research  demonstrated  that  solar  irradiance  and  temperature  variabil¬ 
ity  from  the  provided  solar  irradiance  and  temperature  model  resources  will  result  in  inac¬ 
curate  power  production  predictions  by  the  HOMER  model.  He  also  demonstrated  that  the 
percentage  of  solar  irradiance  disparity  does  not  account  for  the  percentage  of  inaccuracy 
created  in  predicted  power  output.  While  the  original  user-provided  solar  irradiance  profile 
overestimated  actual  solar  irradiance  by  9.35  percent,  this  did  not  equate  to  the  model  over¬ 
estimating  PV  power  production  by  9.35  percent.  As  Newell  pointed  out,  this  is  "due  to  the 
nature  of  the  comparison,  which  only  compared  the  total  kW  per  m  squared  for  the  month 
and  disregarded  when  the  disparities  occurred  [32]."  The  amount  of  solar  irradiance  vari¬ 
ability  that  occurs  for  a  given  scenario  does  not  necessarily  equate  to  an  equivalent  amount 
PV  power  production  variability,  but  there  is  a  direct  relationship  between  solar  irradiance 
and  PV  production  variability  that  must  be  identified.  In  the  previous  chapter,  it  was  men¬ 
tioned  that  HOMER  uses  the  Graham-Holland  method  to  handle  solar  irradiance  profiles 
provided  in  terms  of  average  monthly  irradiance.  While  the  use  of  Graham-Holland  pro¬ 
vides  for  the  generation  of  "synthetic  sets  of  hourly  solar  irradiation  values,"  utilizing  such 
a  synthetic  solar  irradiance  profile  does  not  account  for  inter-annual  variability  [25].  Even 
the  TMY  solar  irradiance  profiles  provided  by  NREE,  which  may  be  expected  to  exhibit 
higher  fidelity  than  the  Graham-Hollands  derived  profiles,  can  not  diminish  the  effects  of 
real-world  inter-annual  solar  irradiance  variability.  Eor  these  reasons,  an  appreciation  of 
the  relationship  between  solar  irradiance  variability  and  HOMER  model  power  estimate 
variability  must  be  acquired  through  a  DOE  leveraging  the  model  itself. 

3.2.1  NREL  Report  on  Temporal  and  Spatial  Solar  Irradiance  Vari¬ 
ability 

In  a  2010  NREE  report  on  the  variability  of  solar  irradiance,  researchers  analyzed  eight 
years  worth  of  measurements  from  the  National  Solar  Radiation  Database  (NSRDB)  in 
terms  of  temporal  and  spatial  variability  [33].  Temporal  variability  refers  to  the  interannual 
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variability  of  solar  irradiance  measurements  for  set  locations  and  times  of  the  year.  Spatial 
variability  refers  to  the  variability  faced  when  a  solar  irradiance  profile  is  acquired  at  a 
given  measurement  site  and  the  profile  is  attempted  to  be  utilized  for  locations  of  increasing 
distance  from  the  measurement  site.  Both  of  these  types  of  solar  irradiance  variability 
influence  decision  making  on  where  and  when  to  employ  solar  power  systems  in  the  solar 
industry. 

In  this  report,  researchers  Wilcox  and  Gueymard  quantified  temporal  variability  as  the  inter¬ 
annual  variability  of  mean  monthly  DNI  coefficient  of  variation  (CV)  (percent)  across  the 
eight  years,  1998-2005.  When  considered  across  the  full  year,  the  DNI  CV  values  range 
from  0.49  percent  to  15.8  percent.  Figure  3.2  shows  how  temporal  solar  irradiance  vari¬ 
ability  is  affected  substantially  by  both  climate  and  season.  This  graphic  shows  that  the 
winter  seasons  have  substantially  higher  inter-annual  solar  irradiance  variability  than  the 
summer  months.  The  West  Coast  and  northern  regions  of  the  United  States  similarly  dis¬ 
play  higher  degrees  of  inter-annual  solar  irradiance  variability  than  the  other  parts  of  the 
country.  These  findings  are  useful  to  the  Marine  Corps  and  other  institutions  that  employ 
solar  power  assets  as  users  can  see  that  solar  irradiance  profile  forecasts  for  the  winter  and 
west  coast/northern  regions  of  the  United  States  are  more  susceptible  to  profile  variabil¬ 
ity  due  to  temporal  solar  irradiance  variability  than  are  profile  forecasts  for  regions  in  the 
summer  and  midwest  our  southwest  regions  of  the  United  States. 

For  these  and  other  reasons,  temporal  solar  irradiance  variability,  or  inter- annual  solar  irra¬ 
diance  variability,  is  a  primary  focus  of  this  research.  While  Wilcox  and  Gueymard  provide 
valuable  insight  into  the  relationship  between  temporal  solar  irradiance  variability  and  the 
factors  of  climate  and  season,  their  study  was  limited  in  scope  to  eight  years  and  did  not 
directly  identify  the  effects  of  temporal  variability  on  the  HOMER  model.  Spatial  solar 
irradiance  variability  was  not  included  in  this  research’s  experiments.  That  being  said, 
understanding  spatial  variability  provides  valuable  context  for  evaluating  the  options  for 
Marine  Corps  employment  of  NREL  TMY  solar  irradiance  profiles. 

By  quantifying  spatial  solar  irradiance  variability  the  Marine  Corps,  and  other  institutions 
interested  in  forecasting  solar  irradiance,  can  identify  the  number  of  solar  irradiance  mea¬ 
suring  stations  required  to  characterize  and  quantify  a  region’s  solar  irradiance  profiles. 
This  knowledge  can  assist  in  determining  the  degree  of  confidence  with  which  users  can 
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Figure  3.2:  Monthly  temporal  CV  (percent)  of  DNI  across  the  United  States,  1998-2005,  from 
[33]. 


utilize  existing  data  for  modeling  solar  irradianee  and  the  assoeiated  performanee  of  pho- 
tovoltaies.  Wileox  and  Gueymard  found  that  spatial  variability  varies  aeross  the  United 
States  with  elimates  and  is  partieularly  sensitive  to  a  loeation’s  proximity  to  the  west  eoast. 
The  degree  of  spatial  variability  aeross  the  United  States  ean  be  seen  in  Figure  3.3.  In 
Figure  3.3  variability  is  quantified  as  CV  of  solar  irradianee  in  terms  DNI,  for  whieh  solar 
irradianee  is  "always  substantially  more  variable  than  that  in  GHI  [33]."  This  CV  measures 
inter-annual  variability  of  daily  mean  DNI  measures  in  50x50  km  areas  over  eight  years. 
The  results  therefore  "roughly  represent  an  area  within. ..25km  of  a  measurement  site  [33]." 
A  more  detailed  deseription  of  the  methodology  ean  be  seen  in  [33]. 

The  analysis  eondueted  by  Wileox  and  Gueymard  shows  that  spatial  variability  of  solar 
irradianee  data  depends  on  season  and  elimate.  Figure  3.2  shows  that  spatial  variability 
aeross  the  United  States  is  at  its  lowest  in  July  and  August,  with  most  DNI  CV  values  in 
the  2-10  pereent  range,  and  at  its  highest  in  November  and  Deeember,  with  most  DNI  CV 
values  in  the  15-50  pereent  range.  Figure  3.2  also  demonstrates  that  elimate  affects  spatial 
variability  and  temporal  variability  similarly,  with  the  West  Coast  and  northern  regions  of 
the  United  States  displaying  higher  spatial  variability  than  other  regions. 

The  take  away  from  the  Wilcox  and  Gueymard  study  on  spatial  variability  of  DNI  is  that 
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Figure  3.3:  Monthly  spatial  CV  (percent)  of  DNI  across  the  United  States,  1998-2005.  Seasonal 
effects  are  visibly  similar  to  changes  in  mean  temporal  DNI  CV  seen  across  the  seasons  and 
climates  of  the  United  States  in  Figure  3.2,  from  [33]. 

users  can  estimate  spatial  variability  of  solar  irradiance  data  to  determine  the  usefulness 
of  measurement  sites  in  forecasting  solar  irradiance  profiles  of  nearby  locations.  Solar  ir¬ 
radiance  data  for  locations  in  much  of  Southern  California  in  August,  for  example,  can 
expect  locations  within  25km  of  a  measurement  site  to  have  variability  (CV)  of  approxi¬ 
mately  2-5  percent  from  the  DNI  at  the  measurement  site.  For  Washington  in  November  or 
December  locations  within  25km  of  a  measurement  site  can  expect  to  have  variability  of 
approximately  30-50  percent  from  to  DNI  at  the  measurement  site.  This  knowledge  pro¬ 
vides  some  context  to  the  discussion  of  Marine  Corps  employment  of  the  HOMER  model 
and  the  NREL  TMY3  database. 

The  nature  of  Marine  Corps  expeditionary  operations  makes  the  effects  of  temporal  and 
spatial  solar  irradiance  variability  even  more  important  than  in  the  commercial  settings  of 
the  solar  industry.  The  short  term  duration  of  expeditionary  operations  makes  the  Marine 
Corps-employed  PV  systems  more  susceptible  to  the  effects  of  temporal  solar  irradiance 
variability.  Although  the  NREL  research  was  conducted  using  DNI  rather  than  GHI,  which 
is  the  measurement  used  for  this  research,  the  variations  in  variability  are  assumed  to  be 
proportional. 

Spatial  variability  of  solar  irradiance  is  particularly  important  because  the  Marine  Corps 
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units  must  select  where  to  position  themselves  based  on  tactical  concerns  rather  than  where 
solar  irradiance  measurements  were  taken.  While  this  research  focuses  on  the  effects  of 
temporal  solar  irradiance  variability,  the  spatial  variability  findings  of  the  NREL  study 
provide  a  context  for  understanding  the  scope  and  range  of  temporal  solar  irradiance  vari¬ 
ability.  Understanding  spatial  solar  irradiance  variability  will  help  the  Marine  Corps  better 
understand  the  applicability  of  solar  irradiance  measurements  for  regional  analysis  of  the 
solar  resource. 


3.3  Temperature  Variability 

The  previous  chapter  explained  the  effects  of  temperature  on  PV  array  performance.  The 
effects  of  interannual  temperature  variability  can  be  inferred  from  that  discussion  of  tem¬ 
perature  effects,  but  Newell’s  previous  research  provides  a  more  specific  example  of  the 
effects  of  temperature  variability.  Newell’s  research  demonstrated  a  9.35  percent  discrep¬ 
ancy  in  total  solar  irradiance  input  resulting  in  a  4  percent  overestimation  in  the  HOMER 
model’s  prediction  of  power  production.  In  that  experiment,  the  effect  of  temperature  on 
the  model’s  overestimation  of  power  was  even  greater.  The  inclusion  of  the  measured 
temperature  profile  for  the  simulation  demonstrated  that  employing  the  synthetic  TMY 
temperature  profile  accounted  for  a  6  percent  HOMER  model  overestimate  in  PV  power 
production.  Temperature  effects  on  a  PV  cell’s  ability  to  produce  power  were  discussed  in 
Chapter  2,  and  this  experiment  provides  an  example  of  why  the  ability  to  model  the  effect 
of  temperature  on  PV  cells  was  included  the  HOMER  model.  The  experiment  also  exem¬ 
plifies  the  importance  of  understanding  the  effects  of  temperature  variability  on  HOMER 
model  power  predictions. 

3.4  Effects  of  Load  Profile  Variability  on  HOMER 

Uncertainty  in  the  HOMER  model’s  predictions  of  fuel  savings  comes  both  from  the  user’s 
inability  to  predict  exact  profiles  for  solar  irradiance  and  ambient  environmental  tempera¬ 
ture  and  from  the  inability  of  the  user  to  anticipate  their  exact  load  profile.  Marine  Corps 
logisticians  and  engineers  currently  build  unit  power  plans  based  on  peak  load  capacities  of 
the  unit’s  equipment.  Eoad  profiles  are  not  employed  because  it  is  difficult  to  project  load 
profiles  from  one  operation  to  another.  This  concern  over  variability  between  user  specified 
load  profiles  and  reality  was  highlighted  in  a  2007  Institute  of  Electrical  and  Electronics 
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Engineers  (IEEE)  report  assessing  the  HOMER  model’s  performance  with  regard  to  wind 
turbine  power  planning  [34].  In  this  report,  Zdenko  Simic  identified  three  modeling  as¬ 
sumptions  requiring  follow-on  research  to  "verify  sensitivity  of  presented  results  on  some 
modeling  assumptions  [34]."  One  of  those  modeling  assumptions  Simic  cited  as  requiring 
verification  was  the  potential  for  variation  from  "assumed  load  schedule  [34]."  The  lack  of 
requisite  measured  experimental  load  profiles  prevented  this  research  from  analyzing  the 
HOMER  model’s  robustness  with  regard  to  load  profile  variability.  All  analysis  conducted 
during  this  research  employed  synthetic  load  profiles  which  will  be  described  in  the  next 
chapter. 
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CHAPTER  4: 

Methodology 


4.1  Experiment  Design 


To  analyze  the  effects  of  solar  irradiance  variability  and  temperature  variability  on  the  ro¬ 
bustness  of  the  HOMER  model,  a  DOE  was  constructed  that  analyzes  the  HOMER  outputs 
for  different  scenarios  simulated  for  each  year  of  up  to  50  years.  This  was  accomplished 
through  the  synthesis  of  50  years’  worth  of  NREE’s  NSRDB  data.  While  an  analysis  of  the 
effects  of  variability  in  load  profile  input  was  originally  intended,  a  lack  of  historical  data 
that  could  represent  the  accurate  depiction  of  load  profile  variability  prevented  its  inclusion 
in  the  experiments.  The  DOE  instead  focuses  on  quantifying  the  robustness  of  HOMER 
with  regard  to  the  inter-annual  variability  of  solar  irradiance  and  temperature  profiles  as 
resources  for  the  model.  This  DOE  employs  a  variety  of  factors  that  could  reasonably  be 
expected  to  impact  the  solar  power  produced  by  a  PV  array.  By  including  these  factors,  this 
DOE  supports  an  analysis  of  which  scenarios  are  more  affected  by  increased  inter-annual 
fluctuations  in  solar  irradiance  or  temperature.  Understanding  how  different  scenarios  are 
influenced  by  these  variations  will  assist  the  Marine  Corps  in  determining  not  only  the  gen¬ 
eral  robustness  of  the  HOMER  model  but  also  how  best  to  employ  the  model  across  the 
range  of  expeditionary  operations. 

Eactors  included  in  this  DOE  can  be  classified  as  either  micropower  system  design  factors 
or  environment  design  factors.  Micropower  system  factors  include  the  size  of  the  unit  (rep¬ 
resented  by  peak  power  level),  the  number  and  size  of  diesel  generators  in  the  micropower 
system,  and  the  number  of  solar  power  assets  utilized.  Environmental  factors  include  the 
physical  location  of  the  micropower  system,  the  duration  of  the  micropower  system  run 
evaluated,  and  the  season,  represented  by  the  starting  date  of  the  micropower  system  run 
evaluated.  Eigure  4.1  shows  each  of  the  DOE  factors  and  the  levels  across  which  each  is 
varied.  The  sections  below  describe  each  factor’s  levels  and  how  these  levels  were  selected. 


35 


DOE  Factor 

Number  of 
Factor  Levels 

Factor  Description 

Location 

10 

Locations  were  identified  across  the  United  States  based  on  dispersion  and 
completeness  of  data 

Duration 

7 

Duration  effects  are  assessed  by  running  simulations  with  durations  of  15 
days,  and  1,2, 3, 4, 5,  and  6  months 

Season 

24 

Season  is  accounted  for  by  utilizing  scenario  start  times  evenly  distributed 
across  the  year 

Unit  Size 

5 

Unit  size  is  represented  by  varying  the  size  of  the  simuiated  microgrid  peak 
ioad  ievels  between  scenarios 

Generator 

Compositions 

4 

Numbers  and  sizes  of  generators  are  varied  for  scenarios  to  investigate  the 
effects  of  pairing  large  and  small  generators 

Solar/Temperature 
Input  Profiles 

3 

1/3  of  scenarios  input  TMY3  solar  profiles  and  measured  temperature  profiles, 
1/3  input  measured  yearly  solar  profiles  and  no  temperature,  and  1/3  input 
correlated  annual  measured  solar  irradiance  and  temperature  profiles 

_ 

Figure  4.1:  The  6  DOE  factors,  number  of  levels  for  each  factor,  and  a  brief  description  of  each 
factor 

4.1.1  Size  of  the  Unit 

Marine  Corps  units  are  modular,  with  larger  units  being  composed  of  smaller,  standard 
sized  units.  An  artillery  battalion,  for  example,  is  comprised  of  a  headquarters  company 
and  three  artillery  companies.  Each  of  the  artillery  companies  are  in  turn  comprised  of  three 
platoons.  The  Marine  Corps  generally  employs  its  forces  in  these  modular  configurations, 
with  the  size  and  composition  depending  on  the  mission.  In  Afghanistan,  battalion-level 
bases  direct  dispersed  company-level  bases  and  these  companies  often  command  dispersed 
platoon  and  squad-sized  outposts.  Each  of  these  dispersed  bases  and  outposts  have  some 
level  of  energy  consumption  that  has  to  be  met,  and  somewhere  along  the  way  a  leader  made 
the  decision  on  what  energy  equipment  was  provided  to  support  the  unit.  The  HOMER 
model  would  be  employed  to  assist  those  leaders  in  determining  how  to  most  efficiently 
allocate  their  power  equipment  between  their  outposts  and  bases.  Because  the  HOMER 
model  would  be  used  across  this  variety  of  unit  sizes  and  power  requirements,  the  size  of 
the  unit  is  included  as  a  factor  in  this  DOE. 

The  different  unit  sizes  are  represented  by  varying  the  peak  power  requirements.  This 
method  was  employed  in  the  E20’s  MEHPS  AOA  which  was  briefed  to  industry  on  31 
January  2013  [16].  Unit  sizes  utilized  for  the  MEHPS  AOA  and  discussed  in  the  brief  in¬ 
clude  a  platoon  outpost,  a  company  COC,  and  a  battalion  COC.  Replicating  the  load  peak 
power  specifications  used  by  E20,  this  research  will  represent  a  platoon’s  load  profile  as 
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Figure  4.2:  Generator  usage  chart,  based  on  767  historical  records  of  real  world  lOkW,  30kW, 
60kW  and  lOOkW  tactical  quiet  generators.  Probability  refers  to  the  likelihood  of  the  load 
equaling  the  given  percentage  of  peak  power,  from  [16]. 


having  a  peak  power  of  3kW,  a  eompany’s  load  profile  with  5kW  or  lOkW  peak  powers 
values,  and  a  battalion’s  load  profile  with  30kW  or  60kW  peak  power  values.  For  an  iden- 
tifieation  of  the  equipment  used  to  derive  these  peak  power  levels,  the  user  should  refer  to 
the  MEHPS  AOA  brief.  The  load  profile  for  eaeh  unit  size  is  ereated  using  the  modified 
load  profile  pereentages  deseribed  in  the  previous  ehapter. 


4.1.2  Diesel  Generators 

Units  that  employ  a  single  generator  rated  to  their  peak  power  ean  often  be  expeeted  to 
eonsume  more  fuel  than  units  with  the  multiple  smaller  generators  summing  to  the  same 
peak  power  level.  This  is  due  to  how  the  fuel  effieieney  of  internal  eombustion  generators 
deereases  as  the  size  of  the  load  deereases  [35].  Employing  multiple  smaller  generators 
rather  than  a  single  large  eapaeity  generator  often  "allows  the  use  of  one  or  two  units  at 
full  load  rather  than  a  larger  unit  at  redueed  load  [35]."  Eigure  4.2  shows  the  probability  (y- 
axis)  for  lOkW,  30kW,  60kW  and  lOOkW  Taetieal  Quiet  Generators  (TQG)  to  be  utilized 
with  various  loads,  with  loads  speeified  as  a  pereentage  of  generator  peak  power  eapaeity. 
This  figure  shows  that  mieropower  system  loads  rarely  reaeh  the  peak  power  level.  Indeed, 
Eigure  4.2  shows  that  generators  have  a  higher  probability  of  being  run  at  load  levels  less 
than  half  their  peak  load  eapaeity  levels  than  they  have  of  being  run  at  load  levels  above 
half  their  peak  load  eapaeity  levels. 
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Figure  4.3:  Composition  of  diesel  generator  and  GREENS/SPACES  systems  factors  for  DOE 
energy  generation  scenarios.  "Diesel  Generators"  refers  number  of  Marine  Corps  TQGs  included 
in  each  scenario  design. 


Units  with  multiple  smaller  generators  can  run  only  those  generators  needed  to  meet  their 
current  power  requirements  at  any  given  time  and  no  more.  This  was  seen  in  the  MEHPS 
AOA  consideration  of  fuel  consumption  for  three  micropower  systems  with  peak  power 
levels  of  60kW  over  a  120  day  period.  The  system  powered  by  a  single  60kW  generator 
consumed  170  gallons  more  than  the  system  powered  by  two  30kW  generators  and  683 
gallons  more  than  the  system  powered  by  one  30kW  generator  and  two  15kW  generators 
[16].  To  account  for  these  effects  in  the  DOE,  the  number  and  type  of  diesel  generators 
in  each  micropower  system  is  included  as  another  design  factor,  with  each  of  the  units 
having  three  unique  compositions.  The  3kW  platoon  sized  unit  is  the  only  unit  that  only 
includes  two  diesel  generator  composition  options,  because  of  the  small  size  of  the  load. 
The  generator  profiles  can  be  seen  in  Eigure  4.3. 

Including  a  generator  within  the  HOMER  model  requires  an  accurate  representation  of 
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certain  specifications.  The  fuel  curve  intercept  and  slope  are  of  particular  interest  as  they 
dictate  the  speed  at  which  the  generator  will  consume  fuel  at  different  loads.  The  technical 
characteristics  for  each  of  the  3,  5,  10,  15,  30  and  60kW  generators  was  represented  in  the 
HOMER  model  through  the  use  of  fuel  consumption  data  measured  at  Aberdeen  Proving 
Ground  and  provided  by  Harold  Scott  Coombe  of  ONR  on  27  June  2013.  All  generators  are 
Marine  Corps  Advanced  Medium  Mobile  Power  Sources  (AMMPS)  generators,  except  for 
the  3kW,  which  is  a  TQG  generator  due  to  the  current  absence  of  a  3kW  AMMPS  generator 
and  a  need  to  simulate  load  profiles  and  generators  below  5kW. 

4.1.3  Solar  Power  Assets 

As  discussed  earlier,  the  Marine  Corps  has  developed  and  fielded  two  solar  power  genera¬ 
tion  assets;  the  GREENS  and  SPACES  systems.  The  GREENS  and  SPACES  PV  array  and 
component  battery  technical  specifications  were  provided  by  ONR  to  support  this  research. 
Both  the  GREENS  and  SPACES  systems  were  designed  to  be  employed  for  stand-alone 
rather  than  grid-tie  functions.  While  the  GREENS  system  generates  enough  power  that  it 
can  offer  significant  fuel  savings  when  employed  in  a  micropower  system,  the  SPACES  sys¬ 
tem  does  not  provide  enough  energy  to  provide  noticeable  savings  relative  to  the  fuel  con¬ 
sumption  of  micropower  systems  supporting  loads  of  3kW  to  60kW.  The  SPACES  system, 
with  its  provision  of  up  to  320W  of  continuous  power,  is  usually  employed  to  power  com¬ 
munications  equipment  and  charge  batteries  for  smaller  units.  An  additional  key  difference 
between  the  GREENS  and  SPACES  systems,  mentioned  in  the  Marine  Corps  Warfighting 
Eab  guide  for  employing  renewable  technologies,  is  that  "the  GREENS  is  not  as  sensi¬ 
tive  to  optimal  conditions  as  SPACES;  as  it  can  harvest  energy  in  less  ideal  conditions" 
than  SPACES  can  [11].  Early  analysis  of  the  fuel  savings  in  micropower  systems  derived 
from  employing  the  SPACES  system  also  showed  that  its  limited  power  capabilities  made 
it  much  less  useful  in  directly  supporting  a  micropower  system  load.  Eor  these  reasons,  the 
SPACES  was  not  included  in  this  research’s  investigation  of  the  HOMER  model’s  robust¬ 
ness  regarding  inter-annual  variations  in  solar  irradiance  and  temperature. 

Eor  this  DOE,  the  number  of  GREENS  systems  employed  per  unit  size  are  shown  in  Eig- 
ure  4.3.  With  its  capability  of  providing  300W  of  continuous  power,  and  up  to  lOOOW  in 
peak  power,  a  single  GREENS  system  can  be  used  to  significantly  support  the  load  require¬ 
ments  of  platoon  (3kW)  or  company-size  (5kW)  units  [11].  The  GREENS  system’s  mod- 
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ular  design  allows  for  its  employment  in  supporting  larger  company,  company-reinforced 
(10-30kW)  and  battalion-size  (60kW)  load  requirements  as  well.  For  these  larger  units  the 
number  of  GREENS  systems  was  selected  to,  as  much  as  possible,  maintain  the  GREENS 
power  production  to  load  peak  power  ratio  seen  in  the  MEHPS  medium  system  in  the 
MEHPS  AOA  brief  [16].  Additionally,  GREENS  system  is  currently  capable  of  combining 
only  five  GREENS  units  under  one  converter/inverter,  limiting  the  number  of  GREENS 
allocated  to  the  larger  30kW  and  60kW  sized  units  to  only  five. 

Adjusting  the  GREENS  system  size  across  the  range  of  unit  sizes  decreases  concerns  of 
having  the  micropower  system  size  act  as  a  confounding  variable  in  the  analysis  of  solar 
irradiance  variability  effects  relative  to  unit  size.  As  the  Marine  Corps  would  be  using  the 
HOMER  model  for  all  unit  sizes  assessing  such  a  variety  of  unit  micropower  system  sizes 
with  HOMER  is  essential  for  assessing  the  robustness  of  the  HOMER  model. 

4.1.4  GREENS  Converter 

Modeling  the  GREENS  system  included  not  just  the  modeling  of  the  PV  arrays  for  each 
system,  but  also  their  batteries  and  converters.  While  specifications  needed  for  the  model¬ 
ing  of  these  components  were  provided  by  ONR,  questions  regarding  the  simulated  quan¬ 
tity  of  GREENS  systems  and  certain  component  specifications  required  additional  steps. 
In  terms  of  quantities  of  GREENS  systems  simulated,  the  assumption  was  made  for  this 
research  that  the  GREENS  system  is  employed  as  a  complete  individual  unit  or  as  a  combi¬ 
nation  of  complete  individual  units.  Although  the  individual  GREENS  system  is  modular, 
it  is  not  represented  as  such  for  the  purposes  of  this  study.  The  only  modular  capability  that 
is  reflected  is  the  ability  of  multiple  GREENS  systems  to  be  employed  to  support  one  load. 
Eollowing  this  representation  of  GREENS  systems  as  being  employed  as  complete  units, 
the  proportion  of  batteries  and  converters  to  PV  arrays  was  kept  constant  as  the  number  of 
GREENS  systems  was  changed  for  different  scenarios. 

Modeling  the  GREENS  converter  demanded  the  acquisition  of  more  detailed  converter 
specifications.  Specifically,  the  inability  of  the  converter,  when  acting  as  an  inverter,  to 
"operate  simultaneously  with  an  AC  generator"  proved  to  be  an  important  specification 
to  clarify  [30].  Although  the  specifications  provided  by  ONR  dictated  that  the  GREENS 
converter  is  not  currently  capable  of  operating  simultaneously  with  AC  generators,  they 
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suggested  that  sueh  a  eapability  could  soon  be  developed.  Scenarios  simulated  through 
the  use  of  the  HOMER  GUI  showed  that  changes  in  the  grid-tie  capability  of  the  inverter 
had  a  noticeable  impact  on  the  system’s  ability  to  decrease  fuel  consumption  for  certain 
scenarios.  After  consultation  with  Scott  Coombe,  Program  Officer  Code  331,  the  decision 
was  made  to  conduct  this  analysis  with  the  converters  modeled  in  accordance  with  the 
current  specifications,  that  is  to  say  with  the  converters  specified  as  not  capable  of  operating 
simultaneously  with  AC  generators. 

Following  analysis  of  the  HOMER  model’s  robustness  concerning  variability  of  inputs, 
the  opportunity  was  taken  to  compare  the  HOMER  model’s  representation  of  standalone 
and  grid-tie  capable  PV  arrays.  The  simulations  which  included  annual  measured  solar 
irradiance  input  profiles  were  replicated  with  the  only  change  being  that  the  GREENS 
system  converter  was  represented  as  being  capable  of  operating  simultaneously  with  AC 
generators  (grid-tie  capable).  In  one  scenario,  for  example,  if  the  converter  was  simulated 
as  being  grid-tie  capable  the  system  was  able  to  save  2,580  liters  of  fuel.  The  GREENS 
system  simulated  as  not  being  grid-tie  capable  represented  a  savings  of  only  17  liters  of 
fuel  under  the  same  circumstances. 

By  executing  these  additional  grid-tie  system  experiments,  this  research  provides  an  un¬ 
derstanding  of  how  the  HOMER  model  will  benefit  the  Marine  Corps  in  assessing  the 
performance  of  the  current  GREENS  system  as  assessing  a  grid-tie  capable  GREENS  sys¬ 
tem.  This  analysis  provides  the  Marine  Corps  with  an  opportunity  not  only  to  assess  the 
HOMER  model,  but  also  to  assess  the  benefits  of  shifting  to  grid-tie  capable  solar  power 
systems.  As  the  Marine  Corps  continues  its  shift  toward  MEHPS  systems,  it  will  be  helpful 
to  understand  how  analysis  tools  such  as  the  HOMER  model  perform  when  assessing  the 
performance  of  grid-tie  capable  systems.  The  findings  of  this  standalone  versus  grid-tie 
system  analysis  will  be  presented  in  the  analysis  chapter  following  analysis  of  temperature 
and  solar  irradiance  inter-annual  variability  effects  on  the  HOMER  model. 

4.1.5  Location 

The  location  of  the  micropower  system  being  simulated  is  an  important  environmental 
factor  because  it  dictates  the  climate  of  the  environment.  Identifying  how  the  effects  of 
solar  irradiance  variability  associated  with  different  climates  will  help  the  Marine  Corps 
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anticipate  the  robustness  of  HOMER  for  expeditionary  operations  relative  to  the  theater  of 
operations  and  the  climate  of  the  target  location.  While  this  research  focuses  on  the  United 
States  because  of  the  high  quality  data,  future  research  focusing  on  the  impact  of  climate 
can  be  used  to  identify  the  robustness  of  the  HOMER  model  across  regions  of  the  world 
where  the  Marine  Corps  has  a  high  likelihood  of  deploying. 

The  solar  irradiance  variability  visualization  of  the  2010  NREE  solar  irradiance  variability 
report,  Eigure  4.4,  demonstrate  how  interannual  solar  irradiance  variability  is  affected  by 
climate.  The  2012  Eocus  Energy  researchers  who  evaluated  the  characteristics  of  solar  irra¬ 
diance  variability  at  five  locations  across  the  United  States  recommended  that  "the  quantity 
of  sites  be  expanded  to  better  comprehend  the  impact  of  locations  and  climates  upon  so¬ 
lar  variability  [23]."  To  ensure  this  research  includes  the  effects  of  discussed  by  the  Eocus 
Energy  researchers,  this  DOE  includes  10  locations  selected  from  across  the  United  States. 
The  10  cities,  superimposed  on  the  2010  NREE  report’s  irradiance  variability  graphic  in 
Eigure  4.4,  were  selected  based  on  the  completeness  of  their  NSRDB  data  and  their  disper¬ 
sion  in  distance  and  climate  across  the  United  States. 

To  identify  the  effects  of  climate  on  the  inter-annual  variability  in  HOMER  fuel  savings 
predictions  the  10  locations  assessed  in  this  research  were  classified  and  grouped  by  cli¬ 
mate.  The  classification  system  used  was  the  Koppen-Geiger  climate  classification  system. 
The  Koppen-Geiger  system  classifies  climates  by  their  main  climate  which  is  determined 
by  vegetation  group  (arid,  warm  temperate,  snow,  etc.),  precipitation,  and  temperature  [36]. 
A  map  of  the  Koppen-Geiger  climate  classification  for  the  U.S.  is  shown  in  Eigure  4.5. 

Although  there  are  climate  variations  across  almost  all  locations  used  in  this  research,  the 
10  locations  were  able  to  be  grouped  into  3  general  climate  groups:  steppe  climates,  snow- 
fully  humid  climates,  and  warm-temperate  climates.  The  steppe  climate  group  includes 
Albuquerque,  El  Paso,  TX,  and  Boulder.  The  snow-fully  humid  climate  includes  Bismarck, 
ND,  Madison,  WI,  and  Salt  Eake  City.  The  warm-temperate  climate  includes  Seattle,  Ster¬ 
ling,  VA,  Tallahassee,  and  Eugene,  OR. 

4.1.6  Starting  Date 

Starting  dates  are  included  as  an  environmental  factor  in  this  DOE  to  support  the  analysis 
of  seasonal  solar  irradiance  variability  effects  on  the  HOMER  model’s  performance.  The 
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Figure  4.4:  2010  NREL  solar  irradiance  variability  report  visualization  of  1998-2005  solar  irra- 
diance  variability,  adjusted  to  show  the  ten  locations  selected  as  location  factor  levels  for  this 
research,  from  [33]. 


NREL  solar  irradiance  variability  report  that  demonstrated  the  impaet  of  elimate  on  solar 
irradiance  variability  also  demonstrated  the  importance  of  eonsidering  the  season  relative 
to  the  location  [33].  Just  as  10  locations  across  the  U.S.  are  used  in  this  DOE  to  evaluate 
the  effeets  of  spatial  solar  irradianee  variability,  the  starting  dates  of  the  simulations  were 
seleeted  for  this  DOE  for  evaluating  the  effects  of  variability  across  the  seasons  as  well.  To 
ensure  all  seasons  are  eonsidered,  the  levels  for  this  factor  are  the  first  and  fifteenth  day  of 
every  month  of  the  year. 


4.1.7  Duration  of  Operation 

The  final  environmental  faetor,  the  duration  of  the  simulated  operation,  is  eomposed  of 
seven  levels.  The  seven  levels  ineluded  in  this  DOE  are  time  periods  of  either  15  days  or 
one,  two,  three,  four,  five  or  six  months  for  simulated  expeditionary  operations.  Varying  the 
duration  of  the  simulated  operations  (micropower  system  run  times)  in  this  way  represents 
the  wide  range  of  time  frames  of  Marine  Corps  expeditionary  operations.  The  levels  them- 
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Figure  4.5:  Koppen-Geiger  climate  classification  system  applied  to  the  continental  United  States, 
from  [36]. 


selves  were  selected  based  on  the  way  that  E20  defines  expeditionary  time  periods  for 
analysis  of  expeditionary  technologies,  defining  them  as  those  missions  lasting  less  than 
120  days  [37], 

The  longer  time  periods  allow  for  considering  how  HOMER’s  projections  adjust  as  the  op¬ 
eration  moves  between  seasons  of  greater  or  lesser  solar  irradiance  variability.  Shorter  time 
periods  are  expected  to  demonstrate  the  effects  seen  in  locations  containing  seasons  with 
particularly  large  interannual  variability.  To  simplify  the  process  of  data  and  the  manipula¬ 
tion  of  the  HOMER  model  in  support  of  this  DOE,  operations  are  not  included  that  would 
link  data  of  different  calendar  years.  This  means,  for  example,  that  there  are  no  two  month 
operations  simulated  as  occurring  from  1  December  to  1  Eebruary. 
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4.2  Assessing  HOMER  Robustness  with  Regard  to  Solar 
Irradiance  Variation 

One  key  question  asked  by  this  researeh  is  how  eritieal  is  it  that  inter-annual  variation  in  so¬ 
lar  irradianee  is  not  refleeted  in  the  HOMER  model,  whieh  typically  relies  on  TMY  data  to 
provide  the  solar  resource.  It  was  anticipated  that  the  short-term  duration  of  Marine  Corps 
expeditionary  operations,  for  which  the  HOMER  model  would  be  employed  by  the  USMC, 
would  represent  greater  risk  of  variability  in  the  HOMER  fuel  savings  prediction  than  the 
variability  seen  in  the  long-term  scenarios  for  which  HOMER  is  usually  employed.  To 
answer  this  question  the  research  methodology  needed  to  identify  a  measure  of  variation  in 
solar  irradiance  over  the  50  years  of  data  as  well  as  a  measure  of  variation  in  the  associated 
fuel  savings  predicted  by  the  HOMER  model  for  each  of  those  years. 

4.2.1  Acquisition  and  Synthesis  of  Solar  Irradiance  Data 

The  solar  resource  used  in  the  HOMER  model  must  be  loaded  into  the  model  in  the  form 
of  a  list  of  8760  quantifications  of  solar  irradiance  in  kilowatt-hours.  A  source  of  so¬ 
lar  irradiance  data  was  found  in  the  NREE’s  NSRDB.  The  NSRDB  provides  a  variety  of 
meteorological  data  sets  for  locations  across  the  United  States.  The  locations  previously 
specified  in  the  DOE  were  selected  from  the  list  of  solar  irradiance  measurement  sites  in 
the  NSRDB  based  on  the  high  quality  and  completeness  of  their  associated  data. 

All  ten  data  collection  sites  that  were  selected  provided  solar  irradiance  data  from  1961  to 
2010.  The  NSRDB  data  sets  were  acquired  from  the  NSRDB  website.  After  acquiring  the 
500  files,  a  java  script  was  written  to  synthesize  the  .txt  files,  parsing  out  only  the  hourly 
values  for  solar  irradiance  in  terms  of  GHI  and  creating  (500)  8760  value  CSV  files. 

4.2.2  Quantifying  Variability  in  HOMER’s  Predicted  Fuel  Savings 

Quantifying  the  variability  in  HOMER’s  fuel  savings  predictions  first  required  a  quantifica¬ 
tion  of  HOMER’s  fuel  savings  predictions.  To  determine  the  fuel  savings  reaped  from  em¬ 
ploying  the  GREENS,  each  scenario  was  run  twice  for  each  year,  once  with  the  GREENS 
included  in  the  micropower  system  and  once  without  the  GREENS  system.  The  difference 
in  the  micropower  system’s  fuel  consumption  between  these  two  runs  was  calculated  as  the 
fuel  savings  reaped  from  employing  the  specified  number  of  GREENS.  Because  of  the  dif- 
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ferences  in  quantities  of  fuel  eonsumption  between  the  different  mieropower  system  sizes 
(3,  5,  10,  30,  and  60kW)  the  pereent  fuel  saved  (fuel  saved  by  employing  GREENS  di¬ 
vided  by  fuel  eonsumed  without  GREENS)  was  seleeted  as  the  metrie  for  quantifying  and 
eomparing  fuel  savings.  The  analysis  of  variation  in  HOMER  fuel  savings  predietions  was 
therefore  eondueted  using  pereent  fuel  savings  in  this  researeh.  To  identify  the  relationship 
between  inter-annual  variability  in  HOMER’s  predieted  fuel  savings  and  the  inter-annual 
variability  in  solar  irradianee  variability,  a  similar  metrie  had  to  be  seleeted  for  quantify¬ 
ing  variability  for  both.  The  root  mean  square  error  (RMSE)  would  eventually  be  seleeted 
as  the  measure  of  variability,  but  the  first  eonsideration  was  the  CV.  While  eheeking  the 
assumptions  needed  to  employ  the  CV  a  trend  was  found  in  the  solar  irradianee  data  that 
disqualified  the  CV  and  introdueed  a  new  eoneern  for  the  robustness  of  the  HOMER  model 
with  regard  to  the  fidelity  of  its  solar  resouree  input  data.  This  trend  and  its  impaet  on  this 
researeh’s  analysis  of  the  HOMER  model’s  robustness  with  regard  to  inter-annual  solar 
irradianee  will  be  diseussed  over  the  next  5  seetions. 

4.2.3  Inter- Annual  Solar  Irradianee  Trend 

One  well  known  tenet  in  the  modeling  and  simulation  eommunity  is  the  "trash-in-trash- 
out"  understanding  of  modeling.  Avoiding  the  input  of  "trash"  into  a  model  is  important 
for  model  researehers  and  users  alike  to  ensure  there  are  no  eonfounding  faetors  or  trends  in 
the  input  data.  To  eheek  for  any  eonfounding  effeets  in  the  solar  irradianee  input  files,  the 
irradianee  files  for  eaeh  eity  were  analyzed  aeross  the  50  years  for  eaeh  seenario  loeation, 
season  (start  date),  and  duration  aeross  all  50  years.  The  results  of  this  eheek  identified 
a  signifieant  eonfounding  effeet  in  the  solar  irradianee  input  files.  Eor  nine  of  the  ten 
loeations  used  in  this  researeh,  the  majority  of  the  seenarios  demonstrated  a  linear  deeline 
in  surfaee  solar  irradianee  over  the  50-year  period  from  1961  to  2010.  This  trend  was  seen 
in  seenarios  ranging  aeross  all  seasons  and  aeross  all  durations  that  were  investigated.  The 
linear  regression  was  signifieant  in  the  majority  of  oases  for  these  loeations.  The  origin  of 
this  trend  and  the  way  it  is  aeeounted  for  in  this  researeh  will  be  diseussed  over  the  next 
four  seetions. 

4.2.4  Global  Dimming 

After  further  researeh  to  identify  the  souroe  of  this  long  term  trend  of  deolining  solar  ir- 
radianoe  over  the  50  years  of  available  data,  the  trend  was  found  to  be  a  manifestation 
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Station 

1951  -  1983 

1984  -  2010 

W/tm*  decade) 

decade 

W/tm*  decade) 

^</decade 

Hraunschweig 

3.2  ±3.2 

2.8  ±  2.8 

4.3  ±2.4 

3.8  ±2.1 

Hchtclberg 

-3.7  ±  4.6 

-3.3  ±  4.2 

43  ±  3.2 

4.0  ±  2.9 

Hamburg 

-0.6  ±  1.7 

-0.5  ±  1.6 

43  ±  2.6 

4.1  ±2.4 

Hohenpeisscnbcrg 

-4.3  ±2.5 

-3.1±  1.8 

3.6  ±  2.7 

2.6  ±  1.9 

Potsdam 

-l.8±  1.8 

-l.6±  13 

4.7  ±  2.2 

4.1  ±  1.9 

Trier 

-4.1  ±6.0 

-3.3  ±  4.9 

2.7  ±  2.9 

2.2  ±  2.4 

Weihenstcphan 

-3.3  ±3.1 

-2.4  ±  2.3 

4.9  ±  2.7 

3.7  ±  2.0 

Wuer/.burg 

-0.4  ±  3.4 

-0.3  ±2.7 

3.2  ±  2.2 

2.6  ±  1.7 

.Vlean  Anomalies 

-2.1  ±  1.6 

-l.7±  1.3 

4.0  ±  2.0 

3.3  ±  1.6 

Figure  4.6:  Measurements  of  long-term  GHI  that  suggest  a  shift  in  Europe's  solar  irradiance 
trend  from  global  dimming  to  global  brightening,  from  [38]. 


of  a  phenomenon  known  as  "global  dimming  [38]."  Global  dimming  is  deseribed  as  the 
"deeadal  deerease  of  surfaee  solar  radiation  [39]."  The  name  "global  dimming"  was  eoined 
by  Stanhill  and  Cohen  in  their  2001  artiele,  "Global  Dimming:  A  Review  of  the  Evidenee 
for  a  Widespread  and  Signifieant  Reduetion  in  Global  Radiation  [40]."  The  "global"  portion 
of  the  name  does  not  refer  to  the  phenomenon  having  a  universal  effeet  aeross  the  world 
but  instead  refers  to  "the  sum  of  diffuse  and  direet  solar  radiation,"  that  is  to  say  global  ra¬ 
diation  [39] .  The  dimming  refers  to  a  deerease  in  that  global  radiation  whieh  was  measured 
as  oeeurring  at  different  rates  aeross  the  world.  In  the  1980s  and  ’90s,  elimate  modelers 
began  to  distinguish  a  trend  of  deereasing  surfaee  solar  radiation  in  available  surfaee  solar 
radiation  reeords  dating  baek  to  the  1950s.  The  rates  of  the  deerease  differed  by  loeation 
but  were  mostly  seen  as  oeeurring  from  the  1950s  up  to  the  1980s  [38].  Figure  4.6  shows 
some  examples  of  global  dimming  rates,  in  terms  of  Watts  per  meter  squared  per  deeade, 
from  1951  to  1983,  that  were  identified  aeross  Germany  in  the  2014  paper  "Rethinking 
Solar  Resouree  Assessments  in  the  Context  of  Global  Dimming  and  Brightening." 

Figure  4.6  also  identifies  an  apparent  shift  in  the  global  dimming  trend  that  was  identified  as 
having  oeeurred  in  Europe  in  approximately  1985  with  a  shift  from  deereasing  to  inereasing 
solar  irradianee.  This  trend  is  referred  to  as  "global  brightening"  as  researehers  found 
rates  of  inereasing  solar  irradianee  from  the  mid-1980s  into  the  21st  eentury.  The  seeond 
eolumn  of  Figure  4.6  shows  the  reversal  of  rates  to  inereasing  solar  irradianee  from  1984  to 
2010,  showing  global  dimming  reversing  to  global  brightening  for  the  speeified  loeations 
in  Europe. 
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The  decrease  in  surface  solar  radiation  that  characterizes  global  dimming  is  thought  to  be 
caused  by  increases  in  pollution  and  aerosols  into  the  atmosphere  where  the  decrease  in 
solar  irradiance  is  seen.  This  poses  an  interesting  challenge  for  climate  modelers  regard¬ 
ing  international  trends  in  aerosol  usage.  Although  one  nation  may  decrease  its  usage  of 
aerosols,  a  neighboring  nation’s  continued  or  increased  use  of  aerosols  can  have  an  effect 
of  inducing  global  dimming  across  the  region.  The  reversal  in  the  solar  irradiance  trend 
from  global  dimming  to  global  brightening  in  Europe  during  the  1980s-’90s  is  thought  to 
be  largely  attributed  to  Europe’s  significant  reduction  of  aerosol  emissions  [38].  This  rela¬ 
tionship  shows  that  there  is  potential  for  modelers  to  utilize  historical  and  predicted  aerosol 
emissions  by  location  to  forecast  global  dimming  and  brightening. 

4.2.5  Accounting  for  Global  Dimming  in  the  Analysis  Methodology 

The  effects  of  global  dimming  on  the  surface  solar  radiation  are  notable  and  their  effects  on 
the  HOMER  model  warrant  further  discussion  in  this  research.  That  being  said,  the  focus 
of  this  research  is  focused  on  the  HOMER  model’s  robustness  concerning  solar  irradiance 
variation  and  not  the  effects  of  climate  change.  Eor  this  reason,  the  effects  of  global  dim¬ 
ming  were  discounted  from  the  total  inter-annual  solar  radiation  variation.  The  decision 
to  hide  the  variation  imposed  by  global  dimming  was  made  under  assumptions  regarding 
HOMER  employment  and  climate  modeling  of  global  dimming  which  will  be  discussed 
later. 

Recognition  of  global  dimming  as  a  confounding  factor  in  the  solar  irradiance  input  data 
changed  the  way  that  the  HOMER  model’s  robustness  with  regard  to  solar  irradiance  vari¬ 
ability  could  be  assessed.  The  original  intent  was  to  assess  the  HOMER  model’s  robustness 
using  the  CV  in  fuel  savings  across  the  50  years,  as  seen  in  Eigure  4.7,  as  the  measure  of 
performance  to  be  compared  with  the  CV  in  solar  irradiance  in  the  time  frame  across  all 
50  years.  Eor  calculating  the  solar  irradiance  CV,  the  standard  deviation  in  cumulative 
measured  GHI  for  a  given  scenario  across  the  range  of  fifty  years  would  be  represented 
as  sigma.  In  this  equation  n  represents  the  number  of  years  of  complete  data  for  the  given 
scenario,  represents  the  mean  sum  of  hourly  measured  GHI  in  kilowatt-hours  per  square 
meter  (kWh/rn^)  for  the  scenario  across  the  n  available  years,  and  Ei  represents  the  mea¬ 
sured  solar  irradiance  for  the  given  year  i.  This  method  of  quantifying  the  interannual  solar 
irradiance  variability  for  the  different  scenarios  is  an  extension  of  the  method  employed 


48 


^scenario  t  ^ 


1/2 


scenario  ~  ^scenario/ 


i=  1 

scenario  i 


m 


Figure  4.7:  Modification  of  the  Wilcox  and  Gueymard  equation  for  calculating  the  cumulative 
GHI  CV  for  each  scenario  across  the  available  years  of  data,  after  [33]. 

by  Wilcox  and  Gueymard  in  their  analysis  of  solar  irradiance  variability  across  the  United 
States  [33].  The  CV  of  percent  fuel  savings  from  employing  the  GREENS  system  in  a  given 
scenario  would  be  calculated  in  a  similar  manner.  Comparison  of  these  two  CVs,  relative 
to  the  various  DOE  factors,  would  demonstrate  the  effects  of  solar  irradiance  inter-annual 
variations  on  the  HOMER  model’s  fuel  savings  predictions. 

The  CV  would  serve  as  a  suitable  measure  of  variability  for  this  analysis  only  under  the 
assumption  that  there  is  no  trend  in  the  solar  irradiance  data.  The  effects  of  global  dimming, 
seen  in  Eigure  4.8,  required  the  utilization  of  a  measure  of  variability  that  could  account  for 
the  inter-annual  trends  associated  with  global  dimming. 

The  linear  regressions  applied  for  each  location  in  Eigure  4.8  demonstrate  a  need  to  also 
account  for  the  trends  caused  by  global  dimming  relative  to  location.  While  the  majority 
of  locations  and  time  periods  assessed  demonstrated  a  negative  trend  with  a  statistically 
significant  linear  relationship  some  locations,  such  as  Tallahassee,  did  not  demonstrate  such 
a  trend.  Effects  of  global  dimming  on  the  variation  in  inter-annual  surface  solar  radiation 
are  therefore  relative  to  the  location  from  which  the  measurements  were  taken. 

Eigure  4.8  shows  the  cumulative  solar  irradiance,  in  GHI  kilowatt-hours  per  meter  squared, 
measured  at  each  of  the  ten  GHI  measurement  stations  for  the  first  6  months  of  each  year 
from  1961  to  2010.  The  line  superimposed  over  each  plot  shows  the  linear  regression  for 
the  given  location  across  the  50  years  for  which  measurements  are  held  (1961-2010).  As 
these  figures  show,  for  all  locations  except  for  Tallahassee,  the  variation  around  the  mean 
for  each  location  is  much  greater  than  the  variation  around  the  linear  regression  plot. 

Taking  out  the  effects  of  global  dimming  on  inter-annual  solar  radiation  variation  was  ac¬ 
complished  by  quantifying  solar  radiation  variation  as  the  RMSE  of  the  cumulative  solar 
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Figure  4.8:  Linear  regression  plots  show  the  negative  linear  relationship  between  years  since  1961 
and  cumulative  surface  solar  radiation  for  the  first  6  months  of  each  year  per  location.  Top  row: 
Albuquerque,  Bismarck,  Boulder,  El  Paso;  Middle  row:  Eugene,  Madison,  Salt  Lake  City,  Seattle; 
Bottom  row:  Sterling,  Tallahassee 

irradiance  for  each  location,  start  time,  and  duration  specified  in  the  DOE.  RMSE  discounts 
the  effects  of  global  dimming  on  solar  irradiance  variability  by  calculating  the  variation  of 
the  50  surface  solar  radiation  measurements  around  the  scenario’s  line  of  linear  regression 
rather  than  around  the  mean  of  the  50  measurements.  Eater,  when  the  variation  in  each 
micropower  system’s  percent  of  fuel  savings,  the  variation  in  the  percent  fuel  savings  was 
similarly  quantified  as  the  RMSE  of  percent  fuel  savings  across  the  50  years. 

4.2.6  Assumptions  for  Discounting  GHI  Variation  Resulting  from 
Global  Dimming 

Representing  solar  irradiance  and  fuel  savings  variation  in  RMSE  downplays  the  variation 
in  solar  irradiance  due  to  global  dimming.  Had  the  effects  of  global  dimming  not  been  ac¬ 
counted  for  and  removed  from  the  solar  irradiance  variation  through  RMSE,  the  variation 
of  solar  irradiance  over  the  last  50  years  would  appear  to  be  much  greater.  This  varia¬ 
tion  caused  by  global  dimming  was  taken  out  of  the  picture  for  this  research  under  the 
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assumption  that  the  effects  of  global  dimming  will  be  considered  in  the  employment  of  the 
HOMER  model  for  forecasting  PV  performance. 

This  representation  of  solar  irradiance  in  RMSE  also  assumes  that  the  trend  in  global  dim¬ 
ming  will  either  maintain  its  downward  trend  or  that  climate  models  will  succeed  in  fore¬ 
casting  changes  in  global  dimming.  While  the  data  in  Eigure  4.8  show  a  continuing  decline 
in  solar  irradiance  through  2010,  several  reports  highlight  a  reversal  of  the  global  dimming 
trend  in  solar  irradiance,  referred  to  as  "global  brightening." 

4.2.7  Implications  of  Global  Dimming  in  Modeling  of  PV  Arrays 

A  different,  more  significant  implication  of  the  global  dimming  trend  is  a  concern  over 
how  to  employ  TMY3  data  with  the  HOMER  model  to  account  for  a  decreasing  supply  of 
solar  irradiance.  NREE’s  TMY3  profiles  are  currently  employed  in  HOMER  in  a  manner 
that  does  not  allow  users  of  the  HOMER  model  to  account  for  inter-annual  trends  in  so¬ 
lar  irradiance  for  their  location.  This  research  is  primarily  concerned  with  analysis  of  the 
HOMER  model’s  robustness  with  regard  to  expeditionary  operations.  That  being  said,  the 
performance  of  the  model  with  regard  to  long-term  solar  irradiance  trends  is  important  be¬ 
cause  of  its  effect  on  each  solar  power  systems’  performance  over  its  life  cycle.  GREENS 
performance  estimates  based  off  of  ten-year-old  TMY  data,  for  example,  will  overestimate 
the  solar  irradiance  available  to  power  the  system  where  there  is  a  negative  linear  relation¬ 
ship  between  time  and  solar  irradiance.  This  concern  will  be  discussed  later  in  the  thesis, 
after  the  impact  of  interannual  solar  irradiance  variability  is  evaluated. 

4.2.8  Quantifying  Inter-annual  Solar  Irradiance  Variability:  Coeffi¬ 
cient  of  Variation  of  RMSE  (CVRMSE) 

Before  discussing  the  variation  in  HOMER  model’s  power  predictions  due  to  solar  irradi¬ 
ance  variation,  the  variation  seen  in  the  latter  must  first  be  discussed  and  quantified.  The 
units  of  a  RMSE  analysis  are  the  same  as  the  units  of  the  values  being  assessed.  Eor 
variability  in  predicted  fuel  savings  this  means  RMSE  values  are  in  terms  of  the  percent 
difference  of  a  given  prediction  from  the  mean  fuel  savings  predicted  by  the  linear  regres¬ 
sion.  More  specifically,  the  RMSE  value  identifies  the  percent  of  inter-annual  variation  for 
68  percent  of  residuals,  while  95  percent  of  residuals  can  be  expected  to  fall  within  twice 
the  RMSE  value  from  the  mean  fuel  savings. 
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^  Bivariate  Fit  of  Sum  of  Solar  GHI  for  1st  6 
Months  of  Yr  in  {k\Nhlm'^2)  By  Years  Since  1960 


^  Bivariate  Fit  of  Sum  of  Solar  GHI  for  1st  15 
Days  of  Yr  in  (kWh/m*2)  By  Years  Since  1960 


Duration 

Mean  Solar  Irradiance 
(kWh/m''2) 

RMSE 

(kWh/m''2) 

CVRMSE 

(%) 

15  Days 

43.7 

4.65 

0.106 

6  Months 

1013.7 

39.9 

0.0394 

Figure  4.9:  Example  demonstrating  the  effect  of  different  cumulative  solar  irradiance  means  on 
RMSE  representation  of  variability  and  how  CVRMSE  removes  the  effects  of  different  quantities 
of  cumulative  solar  radiation 


The  measure  of  predieted  fuel  savings  was  set  to  be  a  pereentage  in  order  to  remove  the 
effeets  of  varying  levels  of  eumulative  fuel  savings  on  the  fuel  savings  variation.  A  similar 
eoneem  is  seen  in  quantifying  how  the  degrees  of  solar  irradianee  variability  ehange  for 
different  seenarios. 

When  eomparing  RMSE  between  seenarios  of  different  durations,  loeations,  and  seasons, 
the  differenees  in  eumulative  solar  irradianee  make  the  RMSE  less  useful  beeause  its  value 
is  relative  to  the  mean  eumulative  solar  irradianee  of  the  given  seenario.  To  assess  the 
effeets  of  duration,  loeation,  and  start  time  on  solar  irradianee  variability  the  effeet  of  vary¬ 
ing  eumulative  solar  irradianee  must  be  removed.  This  is  aeeomplished  by  employing  the 
CVRMSE  as  the  measure  of  inter-annual  variability  in  solar  irradianee.  CVRMSE  removes 
the  effeets  of  varying  eumulative  solar  irradianee  means  by  dividing  the  RMSE  by  the  solar 
irradianee  for  the  respeetive  seenario.  The  Eigure  4.9  example  demonstrates  the  benefits  of 
this  approaeh  by  eomparing  representations  of  variability  for  the  first  fifteen  days  and  for 
the  first  six  months  of  eaeh  year  in  Albuquerque,  New  Mexieo. 

The  RMSE  of  inter-annual  solar  irradianee  is  expressed  in  units  of  kilowatt-hours  per 
square  meter.  Comparing  the  RMSE  measures  of  variability  for  the  fifteen  day  and  six 
month  duration,  4.65  kWh/rn^  vs  39.9  kWh/rn^,  gives  the  impression  of  six  month  du- 
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Figure  4.10:  The  effect  of  varying  means  on  RMSE  representation  of  variability  and  how  CVRMSE 
removes  these  effects  of  varying  means 


ration  having  a  higher  degree  of  variability.  The  RMSE  does  not  aeeount  for  the  different 
measures  of  mean  solar  irradianee  though.  The  CVRMSE  measures  in  Eigure  4.9  show  that 
the  shorter  duration  seenario  represents  greater  variability  relative  to  the  mean  eumulative 
solar  irradianee  of  the  seenario.  In  faet,  the  shorter  seenario’s  CVRMSE  of  .106  represents 
over  twiee  the  variability  seen  in  the  longer  seenario’s  CVRMSE  of  .0394. 

Expanding  the  seope  of  the  analysis  to  evaluate  the  variability  of  all  seenarios,  Eigure  4.10 
shows  how  CVRMSE  removes  the  effects  of  increasing  cumulative  solar  irradianee  from 
the  measure  of  solar  irradianee  variation  with  regard  to  changing  mean  cumulative  solar 
irradianee.  The  CVRMSE  quantification  of  inter-annual  variability  in  solar  irradianee  seen 
in  Eigure  4.10  removes  the  effects  of  the  increased  cumulative  solar  irradianee  and  allows 
the  user  to  identify  interesting  relationships.  Eor  the  remainder  of  this  thesis,  while  inter¬ 
annual  variability  in  HOMER-predicted  fuel  savings  will  be  quantified  as  RMSE  in  units 
of  percent  fuel  savings,  the  inter-annual  variability  of  measured  solar  irradianee  will  be 
quantified  as  the  CVRMSE  of  measured  cumulative  solar  irradianee. 
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4.3  Assessing  HOMER  Robustness  with  Regard  to  Tem¬ 
perature  Variability 

The  HOMER  model  allows  users  to  include  the  effects  of  ambient  temperature  in  their  sim¬ 
ulations.  Just  as  with  solar  irradiance  the  use  of  typical  meteorological  data  for  temperature 
profiles  provides  additional  uncertainty  associated  with  inter-annual  variability.  In  Chapter 
2  the  effects  of  temperature  on  PV  cell  performance  were  discussed,  and  it  was  identified 
that  temperature  effects  come  not  only  from  ambient  temperature  but  also  from  the  temper¬ 
ature  changes  associated  with  solar  irradiance  striking  the  PV  cell.  The  option  for  including 
ambient  temperature  in  the  model  accounts  for  both  of  these  temperature  influences.  For 
this  thesis,  because  the  solar  irradiance  effects  of  temperature  cannot  be  disassociated  from 
the  power  effects  of  varying  solar  irradiance,  analysis  of  "temperature  variation"  on  the 
HOMER  model  refers  to  the  effect  of  ambient  temperature  variation.  To  accomplish  this 
focus  on  temperature  variation,  solar  irradiance  profiles  are  derived  from  the  TMY3  solar 
irradiance  values  for  the  specific  scenarios  across  all  years.  In  these  temperature  focused 
simulations  inputted  temperature  profiles  are  derived  from  the  NSRDB  temperature  profiles 
for  the  specific  year  being  simulated  for  each  scenario. 


4.3.1  Quantifying  Inter-annual  Temperature  Variability:  RMSE 

With  temperature  means  ranging  above  and  below  zero  degrees  Celsius,  CVRMSE  cannot 
be  employed  for  quantifying  temperature  variation.  This  is  because  CVRMSE  cannot  be 
calculated  when  the  values  that  are  measured  can  be  either  positive  or  negative.  Because  of 
this,  and  the  minimal  linear  relationship  between  scenario  mean  temperature  and  temper¬ 
ature  variability,  temperature  variability  is  quantified  in  this  thesis  in  terms  of  RMSE  with 
units  of  degrees  Celsius. 

Figure  4.11  shows  that  an  inverse  linear  relationship  exists  between  scenarios’  mean  tem¬ 
perature  values  and  temperature  variability  in  RMSE.  Figure  4.12  demonstrates  that  while 
the  linear  fit  has  a  reasonably  strong  fit,  with  an  RSquare  of  0.380,  the  Polynomial  fit  shows 
a  stronger  fit  with  a  RSquare  of  0.534.  This  shows  that  scenarios’  interannual  temperature 
variability  increases  at  a  quadratic  rather  than  linear  rate  as  the  mean  temperature  of  a 
scenario  decreases. 
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Figure  4.11:  Relationship  between  scenarios’  mean  temperature  values  and  variability  in  temper¬ 
ature  in  RMSE 

For  the  purposes  of  this  researeh,  RMSE  is  eonsidered  suffieient  for  quantifying  variability 
in  temperature.  Figure  4.12  shows  the  generally  inverse  relationship  between  temperature 
means  and  temperature  variability  in  RMSE. 

4.3.2  Acquisition  and  Synthesis  of  Temperature  Data 

In  order  to  conduct  an  analysis  of  the  HOMER  model’s  robustness  with  regard  to  tempera¬ 
ture  variability,  a  source  for  hourly  temperature  data  needed  to  be  identified.  To  allow  this 
research  to  also  assess  the  relationship  between  temperature  variability  effects  and  solar 
irradiance  variability  effects,  this  source  needed  to  provide  temperature  measurements  for 
the  same  ten  locations  as  seen  in  the  solar  irradiance  data.  The  NSRDB  contains  the  solar 
irradiance  measurements,  but  it  does  not  contain  measured  temperature  data. 

The  National  Climatic  Data  Center’s  (NCDC)  Climate  Data  Online  system  does  provide 
temperature  data  associated  with  those  same  10  locations.  This  system  provides  Integrated 
Surface  Data  (ISD)  that  includes  temperature  in  the  meteorological  data  section  of  its  mea¬ 
surements.  Additionally,  the  system  identifies  the  locations  where  data  is  collected  by 
USAF  number,  latitude,  and  longitude,  helping  to  ensure  that  the  10  locations  used  for 
temperature  measurements  collection  match  those  used  in  the  solar  irradiance  analysis. 
While  the  NCDC  ISD  provided  the  temperature  profiles  needed  for  this  research,  the  data 
required  significant  synthesis  before  it  could  be  used.  The  two  flaws  in  the  data  were  (1)  a 
lack  of  consistent  "on  the  hour"  measurements  and  (2)  gaps  in  data  throughout  the  years  of 
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Figure  4.12:  Relationship  between  scenarios’  mean  temperature  values  and  variability  in  temper¬ 
ature  in  RMSE  described  with  a  linear  fit  (left)  and  described  with  a  quadratic  fit  (right) 


each  location. 

The  required  synthesis  to  fix  and  properly  identify  the  extent  of  these  flaws  was  accom¬ 
plished  using  two  java  scripts.  The  first  script  accounted  for  missing  temperature  measure¬ 
ments  and  for  those  not  measured  on  the  hour.  Like  with  solar  irradiance  profiles,  tem¬ 
perature  profiles  loaded  into  the  HOMER  model  must  be  provided  with  values  for  every 
hour  on  the  hour  for  the  entirety  of  the  given  time  frame.  Some  of  the  ISD  days  provided 
measurements  for  every  hour  on  the  hour,  but  some  days  of  measurements  provided  more 
or  fewer  values  and  values  at  times  off  the  hours.  To  fix  this,  the  first  script  iterated  through 
every  hour  from  1961  to  2010  for  the  temperature  data  for  each  location,  identifying  the 
temperature  measurement  closest  to  the  given  hour  that  is  within  a  half  hour  of  that  hour. 
This  script  also  flagged  all  hours  lacking  a  temperature  measurement  so  that  the  absence 
of  the  measurement  would  be  considered  when  building  experimental  runs  later.  This  first 
script  provided  an  adjusted  CSV  file  to  was  used  to  facilitate  later  leap  year  and  time  zone 
corrections.  The  second  script  identified  all  missing  data  to  allow  support  the  creation  of  an 
experiment  design  considering  only  the  full  temperature  files.  This  was  done  by  reading  in 
the  adjusted  CSV  file  mentioned  above.  The  script  accounted  for  leap  years  by  bypassing 
all  February  29  values  and,  since  the  ISD  data  is  provided  in  Universal  Time  Code  time. 
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changing  the  time  to  local  time  according  to  the  location. 

Due  to  this  research  including  hourly  data  across  50  years  and  10  locations,  acquiring  per¬ 
fect  data  was  not  possible.  This  analysis  was  therefore  conducted  under  the  assumption  that 
temperature  measurements  obtained  within  30  minutes  of  the  given  hour  have  negligible 
fluctuations  in  value  from  that  hour  mark.  All  dates  and  locations  that  lacked  a  measure¬ 
ment  within  30  minutes  of  any  hour  were  not  included  in  experiments.  A  HOMER  input 
scenario  file  was  written  that  dictated  the  scenario  locations  and  time  frames  to  be  run  as 
well  as  the  particular  years  to  be  included  in  each  scenario’s  simulations.  To  create  this 
input  scenario  file,  a  java  script  was  written  to  compile  the  input  scenario  file  and  only  in¬ 
clude  scenarios  where  at  least  30  years’  worth  of  complete  data  were  available  for  each  time 
frame  and  location.  In  order  to  actually  execute  these  experiments,  an  interface  needed  to 
be  developed  to  provide  access  to  the  HOMER  application  programming  interface  (API). 

4.4  Accessing  the  HOMER  API 

HOMER  EEC  has  continued  to  expand  the  models  application  and  improve  its  GUI  since 
the  software  the  company  was  established  in  1997.  Over  the  course  of  these  improvements 
the  HOMER  model  software  has  been  released  in  multiple  versions.  Version  2.0,  released 
in  2001,  increased  the  model’s  applicability  to  micropower  systems  that  are  connected  to 
an  electrical  grid  [41].  In  the  development  of  version  3,  the  company  detached  the  HOMER 
model  engine  from  the  GUI.  This  means  HOMER  3.0  modeling  and  simulation  engine  was 
left  without  a  graphical  interface  for  accessing  the  API.  Part  of  the  reason  for  this  was 
that  the  HOMER  team  did  not  release  HOMER  3.0  and  instead  progressed  to  developing 
HOMER  4.0.  HOMER  4.0  was  in  beta  development  for  the  duration  of  this  research’s 
experiments. 

In  recent  years,  and  in  coordination  with  ONR  and  E20,  HOMER  EEC  has  sought  to  mod¬ 
ify  the  HOMER  model  to  make  it  meet  specific  Marine  Corps  expeditionary  power  plan¬ 
ning  needs.  This  USMC  version  of  the  HOMER  model  was  developed  by  combining  and 
modifying  capabilities  of  the  HOMER  Explorer  and  HOMER  Professional  versions  [42]. 
Executing  this  research  on  HOMER  4.0  or  the  USMC  version  of  HOMER  was  not  feasible 
due  to  the  status  of  both  versions  in  beta  development.  This  research  was  accordingly  con¬ 
ducted  using  the  HOMER  3.0  engine.  HOMER  3.0  was  sufficient  for  the  purposes  of  this 
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research  since  the  HOMER  engine  itself  was  not  significantly  affected  in  the  modification 
of  the  modeling  and  simulation  software  in  the  development  of  the  USMC  version. 

HOMER  3.0  met  the  requirements  for  analyzing  the  HOMER  engine,  but  a  user  interface 
was  needed  for  accessing  the  API  and  executing  the  DOE.  Running  a  full  factorial  analy¬ 
sis  of  the  HOMER  model  as  specified  requires  the  execution  of  several  million  HOMER 
simulations.  Executing  these  millions  of  HOMER  simulations  would  have  taken  approxi¬ 
mately  771  computer  hours,  not  including  trouble  shooting  time  requirements.  This  would 
not  have  been  feasible  for  my  lone  laptop.  Assistance  was  found  with  the  computer  cluster 
run  by  the  Simulation,  Experiment  and  Efficient  Designs  Center  for  Data  Earming  (SEED 
Center).  The  SEED  Center’s  cluster  was  essential  for  executing  the  DOE  as  it  provided 
for  the  executed  of  these  simulations  in  a  matter  of  hours  rather  than  weeks.  In  order  to 
utilize  the  SEED  Center’s  cluster  a  user  interface  accessible  through  the  command  line  was 
required.  This  java  interface  was  built  over  a  week  in  Boulder,  CO,  with  the  SEED  Center’s 
Steve  Upton  and  HOMER  EEC’s  John  Glassmire. 

During  this  week-long  visit  to  the  HOMER  team  offices  in  Boulder,  Upton  was  instru¬ 
mental  in  the  identification  of  the  HOMER  3.0  API  methods  required  for  manipulating  all 
required  factors.  Upton  then  developed  a  java  interface  that  accesses  the  API  and  provides 
it  with  the  necessary  factor  specifications  and  the  input  files  required  to  execute  the  model 
accordingly.  Greater  detail  regarding  the  construction  and  design  of  the  interface  designed 
and  built  by  Upton  to  support  this  research  can  be  found  in  the  "read  me"  document  he 
wrote  for  the  interface.  This  "read  me"  is  provided  in  appendix  1. 

Eollowing  a  considerable  amount  of  interface  trouble  shooting,  the  final  interface  was  used 
to  execute  the  previously  described  DOE.  The  HOMER  output  for  each  simulation  was 
then  combined  into  large  CSV  files.  These  simulation  output  files,  containing  the  output 
of  hundreds  of  thousands  of  scenarios,  were  used  in  the  R  and  IMP  analysis  discussed 
in  the  previous  sections  of  this  chapter.  SEED  Center  analyst  Mary  McDonald  provided 
considerable  assistance  in  combining  these  files  and  in  calculating  the  several  measures  of 
performance  and  variability  required  for  the  conduct  of  this  research. 
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CHAPTER  5: 
Analysis  and  Findings 


5.1  Analysis  Techniques 

This  analysis  was  conducted  using  JMP  Statistical  Discovery  Software.  IMP  facilitated 
the  analysis  of  inter- annual  variability  in  fuel  savings  predictions,  cumulative  solar  irradi- 
ance  measurements,  and  average  temperature  measurements  across  the  factors  of  location, 
length  of  operation,  start  date,  size  of  the  unit  simulated,  and  the  number  and  size  of  internal 
combustion  generators  included  in  the  micropower  system. 

The  primary  JMP  tools  employed  for  this  analysis  included  partition  trees,  linear  regres¬ 
sions,  and  the  wide  variety  of  JMP  visual  data  representation  capabilities.  For  analysis 
of  solar  irradiance  and  temperature  variability,  partition  trees  were  implemented  first  to 
identify  which  factors  were  most  influential  on  the  inter-annual  variability  of  HOMER  fuel 
savings  predictions.  After  identifying  the  most  influential  factors,  JMP  graphical  displays 
and  linear  regressions  were  employed  to  identify  the  nature  of  the  relationships. 

Where  applicable,  linear  regressions  were  employed  to  quantify  the  linear  relationship  be¬ 
tween  individual  factors  and  the  variation  being  evaluated.  To  ensure  the  thoroughness  of 
this  analysis,  multiple  graphical  depictions  of  the  results  were  constructed  for  each  factor 
and  for  combinations  of  factors  to  facilitate  a  visual  exploration  of  the  data.  This  visual 
exploration  allowed  for  an  identification  of  non-linear  relationships  between  factors  and 
the  measures  of  variability  in  fuel  savings,  solar  irradiance,  and  temperature. 

The  first  section  below  will  present  analysis  of  the  HOMER  model’s  robustness  in  the  way 
of  fuel  savings  estimates’  variability  due  to  inter-annual  solar  irradiance  variability.  This 
analysis  was  accomplished  by  utilizing  the  unique  fifty  solar  irradiance  profiles  for  each 
location  as  described  in  the  previous  chapter,  while  deselecting  the  HOMER  option  of  in¬ 
cluding  temperature  considerations.  The  second  section  compares  fuel  savings  estimates 
derived  using  TMY3  solar  irradiance  profiles  to  those  derived  using  the  NREE  measured 
annual  profiles.  The  third  section  below  presents  analysis  of  the  HOMER  model’s  ro¬ 
bustness  in  the  way  of  fuel  savings  estimates’  variability  due  to  inter-annual  temperature 
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variability.  This  analysis  was  accomplished  by  executing  the  HOMER  model  with  TMY3 
solar  irradiance  files,  selecting  the  HOMER  option  of  including  temperature  considera¬ 
tions,  and  inputting  the  unique  temperature  files  for  each  year  and  location  as  discussed  in 
the  previous  chapter.  The  fourth  section  below  compares  the  variability  in  HOMER  fuel 
savings  estimates  that  is  due  to  inter-annual  solar  irradiance  variability  and  the  variability 
of  HOMER  estimates  that  is  due  to  inter-annual  temperature  variability.  The  fifth  and  fi¬ 
nal  section  below  identifies  how  the  inclusion  of  PV  and  battery  bank  systems  similar  to 
GREENS  in  a  micropower  system  affects  the  benefits  derived  from  using  multiple  different 
sized  generators  over  a  single,  peak  load  capable  generator. 


5.2  Analysis  of  HOMER  Robustness  with  Regard  to  Inter¬ 
annual  Solar  Irradiance  Variability 

The  HOMER  model’s  output  from  each  of  the  simulations  was  compiled  into  single  mea¬ 
sures  of  performance  (MOP)  quantifying  fuel  savings  variability  for  every  scenario  across 
the  range  of  complete  years  of  data  utilized.  Each  MOP  quantifies  the  variation  in  fuel 
savings  derived  from  employing  the  GREENS  system  for  each  scenario. 

As  mentioned  in  the  previous  chapter,  the  RMSE  of  cumulative  percent  fuel  savings  is  the 
MOP  used  to  assess  the  variability  of  HOMER’s  fuel  savings  predictions  relative  to  the 
variability  associated  with  the  effects  of  global  dimming  for  each  location.  Comparing  the 
RMSE  of  each  scenario  allows  this  research  to  analyze  the  HOMER  model’s  robustness 
given  varying  degrees  of  solar  irradiance  variability.  Analyzing  how  the  variation  in  fuel 
savings  predictions  is  affected  by  each  factor  was  conducted  by  first  investigating  the  effects 
of  each  factor  across  all  50  years  of  the  measured  solar  irradiance  profiles  for  each  location. 

The  factors  being  assessed,  as  outlined  in  Chapter  4,  include  the  size  of  the  simulated  unit, 
composition  of  diesel  generators  in  the  system,  physical  location,  duration  of  the  simulated 
operation,  and  the  season(s)  of  the  simulated  operation.  Since  the  MOP  for  this  analysis 
is  the  RMSE  of  fuel  savings  achieved  through  the  use  of  the  GREENS  system  for  each 
scenario,  there  is  only  one  data  point  (RMSE  value)  for  each  set  of  factor  levels. 


60 


^  Partition  for  Fuel(rmse) 


PeakLoadOO.  5) 


PeakLoadO.  10) 


PeakLoad<6 


PeakLoad<30. 5. 3,10) 


All  Rows 


0.045- 

0.04- 

0.035- 

0.03- 

0.025- 

0.02- 


0.01 


1  .11.  jf  Number 

'  RSqiure  RMSE  N  Of  Splits  AJCc 
0.445  0.0066428  5110  2  -36736 


*  All  Rows 

Count  5110  LogWorth  Difference 

Mean  00131491  58462984  0.01281 
StODev  0  0089144 


*PeakLoad(60| 
Count  930 

Mean  0.0026687 

StdOev  0.0042924 

*  PeakLoa(>(30, 5, 3, 10) 

Count  4180  LogWortn  Oifterence 

Mean  0.0154809  274.39793  0.00734 
StODev  0.0079481 

>  Candidates 

,  1  , 

’PeakLoad(30,5) 
Count  2320 

Mean  0.0122144 

StODev  0.0078083 

^PeakLoftO(3.lO) 
Count  1860 

Mean  00195552 
StODev  0  0060024 

^  Candidates 

^  Candidates 

Figure  5.1:  Partition  tree  showing  peak  load  as  the  most  influential  factor  in  determining  RMSE 
of  percent  fuel  savings  with  regard  to  solar  irradiance  variation  effects 


5.2.1  Solar  Irradiance  Variability:  Effects  of  Unit  Size  on  HOMER 
Fuel  Savings  Estimates 

A  partition  tree  was  employed  with  the  RMSE  of  pereent  fuel  saved  as  the  Y  value  and 
with  all  factors  of  the  DOE  as  the  X  variables.  The  factors  included  were  peak  load  of 
the  scenario  load  profile,  duration,  location,  and  operation  start  date.  The  partition  tree  is 
shown  in  Eigure  5.1  with  the  two  most  influential  splits.  The  results  of  this  preliminary 
partition  tree  showed  that  the  peak  load  of  the  micropower  system  has  the  greatest  effect  on 
the  variability  in  HOMER  predictions  of  percent  fuel  savings.  After  identifying  this,  addi¬ 
tional  graphics  were  utilized  to  facilitate  visual  exploration  of  the  nature  of  the  relationship 
between  peak  load  and  variability  in  fuel  savings  predictions. 

Eigure  5.2,  which  includes  RMSE  for  only  the  15-day  and  1-month  scenarios,  shows  that 
the  peak  load  of  the  system  has  a  strong  influence  on  the  variability  of  percent  fuel  savings 
in  RMSE.  The  relationship  between  the  systems  of  different  peak  loads  shows  that  the 
predicted  fuel  savings  of  simulated  3kW  and  lOkW  micropower  systems  are  generally  more 
affected  by  inter-annual  solar  irradiance  variation  than  the  simulated  5kW,  30kW  and  60kW 
systems  are,  with  the  60kW  systems  being  the  least  affected  by  solar  irradiance  variation. 

It  was  expected  that  scenarios  with  different  peak  loads,  with  their  different  generators 
meeting  the  loads  and  different  numbers  of  GREENS  systems  providing  renewable  energy 
to  help  meet  those  loads,  would  have  differing  degrees  of  variability  in  percent  fuel  savings. 
The  reason  for  including  the  different  unit  sizes,  and  representing  them  by  different  peak 
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Figure  5.2:  Graphical  representation  of  the  relationship  between  scenario  peak  load  and  interan¬ 
nual  variability  in  predicted  fuel  savings 


power  loads  and  equipment  eompositions,  was  to  identify  how  other  faetors’  effeets  would 
differ  between  the  different  unit  sizes.  The  faet  that  variability  in  predicted  fuel  savings 
was  at  its  highest  for  the  3kW  and  lOkW  scenarios,  rather  than  increasing  or  decreasing 
relative  to  increasing  or  decreasing  peak  load  levels  of  all  scenarios,  showed  that  some 
characteristic  of  the  different  scenarios  other  than  peak  load  level  was  having  the  greatest 
impact  on  inter- annual  variability  of  fuel  savings  predictions.  After  looking  further  into  this 
relationship  it  was  identified  that  those  scenarios  with  the  greatest  inter-annual  variability 
in  predicted  fuel  savings  were  the  scenarios  with  the  highest  ratio  of  PV  power  production 
capacity  to  peak  load  level. 

Figure  5.3  shows  how  scenarios  with  different  peak  load  levels  differed  in  their  ratio  of 
PV  power  production  capacity  to  peak  load  profile  levels.  Micropower  systems  with  a 
higher  ratio  of  PV  production  capacity  to  peak  load  are  shown  to  be  generally  more  capable 
of  gaining  fuel  savings  from  employing  PV-battery  systems.  Micropower  systems  with 
a  lower  ratio  of  PV  production  capacity  to  peak  load  are  less  capable  of  providing  fuel 
savings  because  they  provide  less  cumulative  PV  power  relative  to  load  requirements.  A 
lower  ratio  also  decreases  the  amount  of  the  load  that  can  be  covered  by  the  PV  arrays 


62 


Peak 

Power 
level  (kW) 

Maximum  Power  Provision  Capacity  for 
PV  Arrays  of  Associated  Scenarios  (kW) 

Ratio  of  Peak  Load  Level  to 

PV  Power  Production 

Capacity 

3 

1.76 

0.59 

5 

1.76 

0.35 

10 

5.28 

0.53 

30 

8.8 

0.29 

60 

8.8 

0.15 

Figure  5.3:  Table  shows  the  maximum  PV  power  available  for  each  scenario  relative  to  the 
scenario's  peak  load  level 

in  these  scenarios  because  of  the  stand-alone  converters  used.  When  using  a  stand-alone 
(not  grid-tie  capable)  converter,  the  PV  array(s)  can  only  provide  power  to  the  load  when 
they  can  replace  the  generator(s)  completely,  this  requires  the  PV  power  available  is  to 
equal  or  exceed  the  load  for  the  given  time  steps.  Because  of  this,  the  systems  with  a 
smaller  ratio  of  PV  power  production  capacity  to  peak  load  have  a  smaller  trade  space  for 
employment  of  PV  power  which  is  where  fuel  savings  variation  can  occur.  The  systems 
with  the  smaller  ratio  are  therefore  more  likely  to  have  scenarios  with  mean  fuel  savings  of 
zero,  and  variability  in  fuel  savings  of  zero  will  be  zero.  This  is  best  seen  with  the  60kW 
scenario,  which  has  the  lowest  ratio  of  PV  power  production  capacity  to  peak  load  and  for 
which  75  percent  of  scenarios  have  a  mean  fuel  savings  predictions  below  3  percent.  For 
all  other  peak  load  scenarios  combined,  75  percent  of  mean  fuel  savings  predictions  are 
greater  than  14  percent. 

For  the  5kW,  30kW  and  60kW  scenarios,  which  all  had  smaller  ratios  of  PV  power  pro¬ 
duction  capacity  to  peak  power  load  than  the  3kW  and  lOkW  scenarios,  a  strong  linear 
relationship  between  mean  predicted  fuel  savings  and  variability  in  predicted  fuel  savings 
is  seen.  Figure  5.4  shows  how  this  linear  relationship  is  present  for  the  5kW,  30kW  and 
60kW  scenarios,  with  RSquare  values  of  .38,  .52  and  .71,  respectively.  Figure  5.4  also 
shows  how  variability  in  predicted  fuel  savings  in  the  3kW  and  lOkW  scenarios  is  much 
less  affected  by  mean  predicted  fuel  savings,  with  RSquare  values  of  .01  and  .02,  respec¬ 
tively.  This  identified  how  differences  in  the  scenario  ratios  of  peak  load  level  to  PV  power 
production  capacity  affect  PV  array  performance  and  fuel  savings  variability.  This  also  pre- 
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Figure  5.4:  Graph  shows  strong  linear  relationship  between  mean  predicted  fuel  savings  and 
variability  in  predicted  fuel  savings  for  5kW,  30kW  and  60kW  scenarios  for  HOMER  prediction 
variability  resulting  from  inter-annual  solar  irradiance  variability 


sented  a  concern  for  how  to  conduct  the  analysis  of  the  other  faetors  in  the  DOE,  ineluding 
duration,  season  and  loeation. 

For  the  analysis  of  the  other  DOE  faetors,  analysis  aeross  the  5kW,  30kW,  and  60kW  see- 
narios  would  suffer  from  the  eonfounding  effeets  of  the  relationship  between  mean  pre¬ 
dieted  fuel  savings  and  variability  of  pereent  fuel  savings.  This  eonfounding  effeet  would 
influenee  the  pereeived  relationship  between  the  assessed  DOE  faetor  and  variability  in 
predieted  fuel  savings.  In  the  end,  it  was  deeided  to  keep  seenarios  of  all  peak  load  levels 
in  the  analysis,  with  the  reeognition  that  the  3kW  and  lOkW  seenarios  would  best  represent 
the  effeets  of  the  subsequent  faetors  to  be  evaluated.  Analysis  of  the  other  faetors  in  the 
DOE  was  eondueted  separately  for  eaeh  of  the  5  different  peak  power  levels.  This  separa¬ 
tion  was  done  to  remove  the  effeets  of  different  peak  load  levels  and  different  PV  power 
to  load  ratios  between  the  seenarios,  minimizing  the  eonfounding  effeets  of  both  peak  load 
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levels  and  mean  predicted  fuel  savings  levels  in  the  analysis  of  other  DOE  factors. 

It  is  important  to  clarify  that  while  the  load  profiles  and  peak  loads  of  these  scenarios  are 
based  on  ONR  and  E20  experiments,  and  while  the  equipment  simulated  for  the  experi¬ 
ments  is  Marine  Corps  equipment,  the  GREENS  system  simulated  in  these  simulations  is 
meant  to  provide  a  representation  of  solar  power  that  would  be  provided  by  other  PV  hybrid 
power  systems.  These  results  therefore  do  not  accurately  represent,  and  are  not  meant  to 
represent,  the  ability  of  Marine  Corps  solar  power  equipment  as  it  currently  exists  to  sup¬ 
port  the  varying  levels  of  load  profiles.  Eor  these  scenarios  it  is  also  worth  reemphasizing 
that  all  variation  is  due  solely  to  inter-annual  differences  in  solar  irradiation,  which  is  the 
only  factor  changed  in  the  model  inputs  across  the  fifty  years  of  data  for  each  design  point. 

5.2.2  Solar  Irradiance  Variability:  Effects  of  the  Duration  of  an  Op¬ 
eration  on  HOMER  Enel  Savings  Estimates 

Having  identified  the  high  degree  of  influence  of  a  micropower  system’s  peak  load  level, 
as  well  as  the  influence  of  the  PV  power  capacity  to  peak  load  ratio,  on  the  variability 
of  predicted  fuel  savings,  partition  tree  analysis  was  used  to  determine  the  second  most 
influential  factor.  This  time  a  partition  tree  was  used  for  the  results  of  each  peak  load 
separately  to  determine  what  factor  was  most  influential  for  scenarios  of  different  peak 
load  levels. 

A  partition  tree,  executed  with  scenarios  of  all  peak  load  levels,  showed  that  the  most 
significant  factor  across  all  scenarios  after  peak  load  level  was  simulation  duration,  with 
an  RSquare  value  of  0.165.  Although  duration  was  found  to  be  the  most  influential  factor 
in  determining  RSME  of  percent  fuel  savings,  the  degree  of  influence  varied  depending 
on  the  peak  load  of  the  micropower  system.  This  can  be  seen  in  Eigure  5.5,  which  shows 
that  duration  explains  more  about  the  RMSE  variability  of  predicted  fuel  savings  when  the 
micropower  system’s  peak  load  is  3kW  or  lOkW.  The  season,  or  start  time,  of  a  simulation 
provides  more  insight  into  the  breakdown  of  RMSE  values  when  the  scenario  peak  load 
is  5kW,  30kW,  or  60kW.  As  discussed  in  the  previous  section,  the  5kW,  30kW,  and  60kW 
scenarios  are  influenced  by  the  confounding  effects  of  mean  predicted  fuel  savings.  The 
confounding  effect  of  mean  predicted  fuel  savings  can  be  seen  through  a  visual  exploration 
of  the  effects  of  simulation  duration  and  start  time,  seen  in  Eigure  5.6  and  Eigure  5.7, 
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Figure  5.5:  Partition  tree  showing  differences  between  peak  load  scenarios  in  what  qualifies  as 
the  second  most  influential  factor  in  determining  RMSE  of  percent  fuel  savings  with  regard  to 
solar  irradiance  variation  effects 

respectively. 

Figure  5.6  and  Figure  5.7  visualize  the  effects  of  duration  and  season  on  mean  HOMER  fuel 
savings  predictions  as  well  as  their  effects  on  the  variability  of  HOMER  derived  percent 
fuel  savings  predictions.  By  visually  assessing  Figure  5.7  across  the  simulation  start  dates 
it  is  seen  that  the  HOMER  variability  associated  with  the  season  factor  is  confounded  by 
the  effects  of  varying  mean  predicted  fuel  savings  across  the  seasons.  The  duration  factor 
shows  no  such  confounding  relationship  between  inter-annual  variability  in  predicted  fuel 
savings  and  variation  in  mean  predicted  fuel  savings. 

Figure  5.7  shows  that  the  season,  or  start  date,  of  simulations  has  a  noticeable  effect  on 
both  the  inter-annual  variability  of  predicted  fuel  savings  and  on  the  mean  predicted  fuel 
savings  of  the  simulations.  The  matching  bi-modal  shape  between  the  HOMER  predictions 
variability  and  HOMER  mean  fuel  savings  graphs  in  Eigure  5.7  suggests  that  the  changes 
in  inter-annual  variability  are  associated  with  the  seasonal  changes  in  mean  predicted  fuel 
savings  levels.  The  relationship  between  the  graphs  can  be  seen  to  be  stronger  for  the  5kW, 
30kW,  and  60kW  scenarios,  which  can  be  ascribed  to  the  previously  discussed  confounding 
effects  of  mean  predicted  fuel  savings.  The  seasonal  variability  in  HOMER’s  predicted 
fuel  savings  for  the  5kW,  30kW,  and  60kW  scenarios  can  therefore  be  ascribed  to  seasonal 
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Mean  Predicted  Fuel  Savings  (%)  &  RMSE  of  Predicted  Fuel  Savings  (Va)  vs.  Duration 

Peak  Load  LeveKkW) 

3  5 


Figure  5.6:  The  relationship  between  simulation  duration  on  mean  predicted  fuel  savings  and  on 
inter-annual  variability  in  predicted  fuel  savings 


variations  in  load  levels  rather  than  seasonal  differenees  in  inter-annual  solar  irradianee 
variability. 

While  Figure  5.6  shows  that  duration  exhibits  a  notieeable  effeet  on  the  variability  of  pre- 
dieted  fuel  savings  between  the  years,  it  does  not  exhibit  a  notieeable  effect  on  the  mean 
predicted  fuel  savings.  This  suggests  that  the  relationship  between  inter- annual  variability 
of  predicted  fuel  savings  and  simulation  duration  is  not  due  to  any  relationship  between 
duration  and  mean  predicted  fuel  savings  levels.  Based  on  this,  and  the  degree  of  influence 
duration  exhibited  on  the  3kW  and  lOkW  scenarios  in  the  partition  trees,  the  duration  of  the 
simulations  is  determined  to  be  a  more  significant  factor  than  seasonality  in  determining 
inter-annual  variability  of  HOMER’s  predicted  fuel  savings. 

As  the  duration  of  a  scenario  has  been  identified  as  the  second  most  influential  factor  on 
the  inter-annual  variability  of  a  scenario’s  HOMER  fuel  savings  predictions,  the  nature  of 
this  factor’s  influence  will  now  be  discussed.  Eigure  5.8  shows  that,  when  considering 
the  effects  of  inter-annual  solar  irradianee  variability  on  projected  fuel  savings  variation, 
an  increase  in  duration  is  associated  with  a  decrease  in  the  variability  of  HOMER  percent 
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RMSE  of  Predicted  Fuel  Savings  (%)  &  Mean  Predicted  Fuel  Savings  (%)  vs. 
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Figure  5.7:  The  relationship  between  simulation  start  time  during  the  year  on  mean  predicted 
fuel  savings  and  on  inter-annual  variability  in  predicted  fuel  savings 


fuel  savings  estimates.  An  inerease  in  the  duration  of  a  seenario  is  also  assoeiated  with  a 
deerease  in  the  range  of  predieted  fuel  savings  variability. 

Considering  the  RMSE  values  for  the  3kW  peak  load  seenarios  helps  to  illustrate  the  vari¬ 
ability  of  HOMER-generated  fuel  savings  estimates  to  be  expeeted  for  operations  of  differ¬ 
ent  durations.  The  inerease  in  the  range  of  the  RMSE  box  plots’  first  and  third  quartiles  as 
durations  of  seenarios  deerease  in  length,  as  seen  in  Eigure  5.8,  illustrates  how  operations  of 
a  longer  duration  have  a  more  predietable  degree  of  inter-annual  variability.  In  the  3kW  see¬ 
nario,  the  RMSE  first  and  third  quartiles  for  pereent  fuel  savings  were  0.00946  and  .0337, 
respeetively  for  15  day  durations;  these  RMSE  quartiles  were  .00565  and  .01886  for  2 
months  durations;  and  these  RMSE  quartiles  were  .00665  and  .01 16  for  6  month  durations. 
As  seen  in  Eigure  5.8,  the  deerease  in  the  degree  of  variation  (RMSE)  and  the  deerease  in 
the  differenee  between  first  and  third  quartiles  of  the  box  plots  is  seen  aeross  seenarios  of 
all  peak  loads  as  seenario  duration  inereases.  This  analysis  supports  the  hypothesis  that 
shortening  the  duration  of  operations  inereases  the  variability  of  HOMER-generated  fuel 
savings  predietions  due  to  inereased  inter-annual  solar  irradianee  variability. 
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Variation  in  Predicted  Fuel  Savings  (RMSE  of  Predited  %  Fuel  Savings)  vs. 
Duration  of  Simulated  Operation  (hrs) 


c^j  Fuel(rmse) 


Duration  of  Simulated  Operation  (hrs) 


Figure  5.8:  RMSE  of  fuel  savings  versus  duration  of  simulated  operation 


The  relationship  between  inter-annual  surfaee  solar  radiation  variability  and  duration  mir¬ 
rors  the  relationship  between  inter-annual  variability  in  HOMER  predicted  fuel  savings 
and  duration,  where  the  inter-annual  variability  in  HOMER  fuel  savings  predictions  is  due 
to  inter-annual  surface  solar  radiation  variability.  These  relationships  can  be  seen  in  Eig- 
ure  5.9.  As  the  duration  of  scenarios  increase  the  variation  in  cumulative  surface  solar 
radiation  decreases.  The  shorter  duration  simulations  show  higher  degrees  of  variability 
in  surface  solar  radiation  just  as  they  represent  higher  degrees  of  variability  in  HOMER- 
generated  fuel  savings  predictions.  Eigure  5.9  shows  this  relationship  between  solar  irra- 
diance  and  duration,  as  well  as  predicted  fuel  savings  and  duration,  in  terms  of  CVRMSE. 
CVRMSE  is  utilized  for  this  comparison  due  to  the  need  to  account  for  increased  solar  irra- 
diance  variation  in  RMSE  that  is  due  to  increased  quantities  of  cumulative  solar  irradiance 
from  longer  durations.  CVRMSE  provides  a  better  representation  of  variability  relative  to 
the  amount  of  cumulative  solar  irradiance  per  scenario.  This  relationship  suggests  that,  by 
identifying  the  inter-annual  variability  of  surface  solar  radiation  associated  with  different 
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Figure  5.9:  Relationship  between  solar  irradiance  variation  (CVRMSE)  and  predicted  fuel  savings 
variation  (CVRMSE)  relative  to  duration 

durations  of  time,  for  a  given  mieropower  system,  loeation,  and  season,  a  user  ean  esti¬ 
mate  the  inter-annual  variability  of  HOMER  fuel  savings  predietions  assoeiated  with  the 
simulated  operation  duration. 


5.2.3  Solar  Irradiance  Variability:  Effects  of  the  Season  of  an  Opera¬ 
tion  on  HOMER  Enel  Savings  Estimates 

The  season  of  simulated  operations  demonstrates  a  minor  influenee  on  the  inter-annual 
variability  of  HOMER  model  fuel  savings  predietions.  As  is  seen  in  Eigure  5.10,  the  season 
of  a  simulated  operation  has  an  inereasing  effeet  on  the  inter-annual  variability  of  HOMER 
model  predieted  fuel  savings  as  the  duration  of  the  simulation  is  shortened.  This  ean  be 
explained  by  noting  that  simulated  operations  of  longer  durations  inereasingly  inelude  more 
than  one  season,  whieh  ean  eause  the  results  to  smooth  out  the  pereeived  effeets  of  the 
different  seasons.  Analysis  of  seasonal  effeets  on  the  inter-annual  variability  of  HOMER 
fuel  savings  predietions  will  therefore  foeus  on  the  fifteen  day  and  one  month  seenarios. 
This  allows  for  a  elearer  identifieation  of  seasonal  effeets  on  inter-annual  variability  of 
HOMER  fuel  savings  predietions. 
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RMSE  of  Predicted  %  Fuel  Savings  vs.  Start  Time  (for  all  scenarios) 
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Figure  5.10:  Relationship  between  season  of  simulation  and  variation  in  HOMER  predictions  of 
fuel  savings 


Figure  5.10  shows  the  degree  to  which  inter- annual  variability  in  predicted  fuel  savings 
fluctuates  with  the  seasons.  These  fluctuations  in  RMSE  of  fuel  savings  predictions  do  not 
relate  to  the  seasonal  fluctuations  in  solar  irradiance  variability.  While  Figure  5.10  shows 
the  seasonal  fluctuations  in  fuel  savings  variability  to  be  bimodal,  Figure  5.11  shows  the 
seasonal  fluctuations  in  solar  irradiance  variability  to  be  unimodal.  Instead,  the  fluctuations 
seen  in  predicted  fuel  savings  variability  mirror  the  seasonal  fluctuations  in  micropower 
system  load  that  is  specified  by  the  user  in  the  model.  The  bimodal  nature  of  inter-annual 
variability  in  HOMER  model  fuel  savings  predictions  reflects  the  bimodal  nature  of  the 
load  profiles  used  as  model  inputs. 

Figure  5.12  shows  the  seasonal  fluctuations  in  the  power  loads  that  were  provided  to  the 
HOMER  model.  For  the  15-day  and  1-month  scenarios  in  particular,  fluctuations  in  the 
inter-annual  variability  of  HOMER-generated  fuel  savings  estimates  mirror  the  seasonal 
fluctuations  in  micropower  system  loads  and  the  subsequent  fluctuations  in  mean  predicted 
fuel  savings  that  come  from  seasonal  changes  in  the  load  profile.  Figure  5.13  highlights  this 
relationship  between  seasons  and  mean  predicted  fuel  savings.  Figure  5.13  also  shows  that 
this  relationship  is  strongest  for  the  5kW,  30kW,  and  60kW  scenarios.  As  discussed  in  the 
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Figure  5.11:  Relationship  between  season  of  simulation  and  variation  in  Surface  Solar  Radiation 


Figure  5.12:  Seasonal  fluctuation  in  power  loads  were  provided  to  the  HOMER  model 

previous  section,  the  5kW,  30kW,  and  60kW  scenarios  are  more  affected  by  the  load  levels 
due  to  their  lower  ratio  of  PV  power  production  capacity  to  peak  load  levels.  This  suggests 
that  seasonal  changes  in  the  variability  of  HOMER-generated  fuel  savings  estimates  are  the 
result  of  seasonal  changes  in  load  profiles  more  than  they  are  the  result  of  seasonal  changes 
in  inter-annual  solar  irradiance  variability. 

The  take  away  from  this  analysis  is  that  seasonal  effects  on  load  profile  fluctuation  have 
a  greater  impact  on  HOMER  model  inter-annual  variability  than  seasonal  effects  have 
on  inter-annual  solar  irradiance  variability.  Seasonal  effects  on  inter-annual  solar  irradi¬ 
ance  variability  demonstrate  minimal  effects  on  the  inter-annual  variability  of  the  HOMER 
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Figure  5.13:  Relationship  between  season  of  simulation  and  variation  in  HOMER  predictions  of 
fuel  savings,  for  15-day  and  1-month  scenario  durations 


model’s  fuel  savings  predietions.  For  the  5kW,  30kW,  and  60kW  seenarios,  this  analy¬ 
sis  demonstrates  a  relationship  between  seasonal  load  profile  ehanges  and  ehanges  in  the 
inter- annual  variability  of  HOMER  model  fuel  savings  predietions.  These  findings  show 
that  seasonal  effeets  of  inter-annual  solar  irradianee  variability  are  negligible  for  the  pur¬ 
poses  of  the  Marine  Corps,  at  least  for  environments  such  as  those  to  be  found  across  the 
continental  United  States.  It  is  worth  noting  that  the  climates  assessed  in  these  simulations, 
which  will  be  discussed  in  the  following  section,  only  allow  for  identification  of  seasonal 
effects  of  inter-annual  solar  irradianee  for  United  States  climates.  These  results  are  not 
representative  of  seasonal  effects  on  inter-annual  solar  irradianee  variability  that  may  be 
seen  in  climates  outside  the  continental  United  States.  With  the  identification  of  how  load 
profiles  affect  variability  in  HOMER  fuel  savings  predictions,  it  is  also  worth  reminding  the 
reader  that  the  load  profiles  used  in  these  simulations  were  derived  from  an  adaption  to  the 
Marine  Corps’  seasonal  load  profile  template.  The  demonstrated  influence  of  load  profile 
seasonal  changes  on  inter-annual  variability  of  HOMER  fuel  savings  predictions  provides 
an  additional  incentive  for  improving  the  fidelity  of  load  profiles  for  HOMER  simulations. 
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Figure  5.14:  Effects  of  climate  on  inter-annual  variability  across  the  seasons,  shown  for  lOkW 
and  15  day  scenarios 


5.2.4  Solar  Irradiance  Variability:  Effects  of  Location  on  HOMER 
Enel  Savings  Estimates 

After  having  identified  the  most  influential  three  factors  of  the  DOE  to  be  the  size  of  the 
load  supported  by  the  micropower  system,  the  duration  of  the  simulation,  and  the  season 
of  the  simulation,  remaining  factors  included  the  scenario  climate  and  the  composition  of 
diesel  generators  within  a  micropower  system.  The  use  of  partition  trees  for  this  analy¬ 
sis  was  hampered  due  to  the  number  of  more  significant  factors  influencing  variability  in 
HOMER  model  fuel  savings  predictions  due  to  solar  irradiance  variability.  Because  of  this, 
visual  exploration  was  used  as  the  method  of  analysis  for  identifying  first  the  effects  of 
climate  and  then  the  effects  of  different  diesel  generator  compositions  on  the  inter-annual 
variability  of  HOMER  model  predictions. 

A  subsection  of  the  total  data  was  selected  for  visual  exploration  of  climate  effects  on 
inter-annual  variability  of  HOMER  fuel  savings  estimates  due  to  inter-annual  solar  irradi¬ 
ance  variability.  Parameters  for  this  subsection  were  selected,  based  on  previous  findings, 
to  best  illustrate  the  effects  of  climate  on  inter-annual  variability  by  using  the  scenarios 
with  the  greatest  inter-annual  variability  due  to  peak  load  and  duration  factor  specifications. 
Eigure  5.14  includes  the  scenarios  with  lOkW  peak  load  levels  and  15  day  durations,  this 
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included  720  data  points,  representing  variability  across  36,000  individual  HOMER  simu¬ 
lations.  With  the  climate  categories  of  "Snow  and  Fully  Humid,"  "Steppe,"  and  "Warm  and 
Temperate"  defined  in  the  previous  chapter.  Figure  5.14  shows  how  the  climate  of  a  loca¬ 
tion  relates  to  the  level  of  inter-annual  variability  in  HOMER  model  fuel  savings  estimates 
across  the  seasons. 

While  the  bi-modal  shape  of  the  plots  shows  the  confounding  effects  of  load  profile  levels 
on  the  measures  of  mean  percent  fuel  savings  estimates  and  inter-annual  fuel  savings  vari¬ 
ability,  the  interesting  information  to  pull  from  these  graphs  is  how  the  climates’  measures 
relate  to  each  other.  Figure  5.14  shows  that  the  steppe  climate  simulations  consistently 
have  mean  percent  fuel  savings  predictions  above  those  of  the  other  two  climate  categories. 
This  matches  some  expectations,  as  steppe  climates  are  generally  characterized  by  greater 
sunlight  and  fewer  clouds  throughout  the  year.  Figure  5.14  also  suggests  that  HOMER 
simulations  for  steppe  climates  have  less  inter-annual  variability  in  fuel  savings  estimates 
for  most  times  throughout  the  year.  It  is  also  interesting  to  note  that  the  inter-annual  fuel 
savings  variability  for  the  steppe  climate  does  not  follow  the  bi-modal  shape  as  much  as 
the  other  two  climates  do.  This  is  particularly  evident  in  the  spring  season,  when  the  steppe 
climate  demonstrates  much  less  inter-annual  variability  in  HOMER  fuel  savings  predic¬ 
tions  compared  to  the  other  climates.  This  is  not  the  case  for  the  winter  season,  when  the 
mean  inter-annual  variability  measures  of  the  steppe  climate  rise  above  the  mean  inter¬ 
annual  variability  measures  of  both  other  climates.  These  findings  are  not  surprising,  but 
they  do  support  the  hypothesis  that  fuel  savings  estimates  of  HOMER  simulations  in  steppe 
climates  include  less  inter-annual  variability  than  do  estimates  for  the  other  two  climates. 

5.2.5  Solar  Irradiance  Variability:  Effects  of  the  Composition  of 

Diesel  Generators  on  HOMER  Enel  Savings  Estimates 

The  final  factor  left  to  assess  in  considering  how  inter-annual  solar  irradiance  variability 
affects  the  inter-annual  variability  of  HOMER  fuel  savings  predictions  is  the  number  and 
size  of  diesel  generators  used  to  power  a  given  load.  Analysis  of  this  generator  composition 
factor  across  each  peak  load  level  showed  that  the  number  and  size  of  diesel  generators 
utilized  for  supporting  a  load  of  a  specified  peak  capacity  had  a  negligible  effect  on  the 
inter-annual  variability  of  the  systems  for  all  peak  load  levels.  To  assist  in  focusing  on 
factors  aside  from  peak  load,  and  given  the  similarity  in  HOMER  results  across  the  multiple 
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Figure  5.15:  Visualization  of  how  inter-annual  variability  in  HOMER  fuel  savings  predictions 
varies  across  the  four  different  lOkW  micropower  system  generator  compositions  for  the  seven 
duration  levels 

peak  load  levels,  analysis  of  the  generator  composition  factor  was  narrowed  down  to  one 
peak  load  level  for  more  specific  analysis.  It  was  identified  earlier  that  the  3kW  and  lOkW 
scenarios  are  preferable  for  continued  analysis,  due  to  their  higher  ratio  of  load  to  solar 
power  generation  capacity.  As  the  lOkW  peak  load  level  is  the  only  load  level,  between 
the  3kW  and  lOkW  scenarios,  that  was  simulated  with  multiple  diesel  generator  size  and 
quantity  configurations,  the  analysis  of  the  lOkW  peak  load  scenarios  is  presented  here. 

To  support  this  analysis  an  additional  generator  composition  was  added  for  the  analysis;  a 
single  lOkW  generator  scenario.  Four  different  generator  compositions  were  therefore  sim¬ 
ulated  for  scenarios  with  a  peak  load  of  lOkW.  These  included  a  scenario  with  one  lOkW 
generator  and  one  5kW  generator,  a  scenario  with  one  lOkW  generator  and  two  3kW  gen¬ 
erators,  a  scenario  with  two  5kW  generators,  and  a  scenario  with  one  lOkW  generator.  In 
Figure  5.15,  the  different  generator  composition  scenarios  with  a  lOkW  peak  load  level 
were  further  separated  by  duration  levels.  This  improves  the  visualization  of  how  gen¬ 
erator  composition  affects  inter-annual  variability  in  HOMER  fuel  savings  predictions.  In 
Figure  5.15,  only  minimal  relationship  is  seen  between  the  generator  composition  scenarios 
and  inter-annual  variability  of  HOMER  fuel  savings  predictions. 

While  generator  composition  does  not  provide  any  interesting  findings  with  regard  to  inter¬ 
annual  variability  of  HOMER  predictions,  this  factor  did  provide  an  opportunity  to  investi¬ 
gate  the  relationship  between  the  generator  composition  of  a  micropower  system  and  fuel 
consumption  for  the  given  micropower  system  and  load  profile.  In  section  4.1.2,  it  was  ex¬ 
plained  that  scenarios  including  the  use  of  multiple  small  generators  have  greater  potential 
for  decreased  fuel  consumption  than  do  those  scenarios  that  include  only  one  large  gen- 
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Figure  5.16:  Fuel  consumed  (liters)  when  micropower  system  includes  GREENS  and  TMY3  solar 
irradiance  profiles  (top)  compared  to  fuel  consumed  (liters)  when  using  no  GREENS  systems 
(bottom) 

erator  supporting  the  same  load.  Comparisons  of  fuel  eonsumption  levels  aeross  the  fore 
mentioned  seenario  generator  compositions  have  supported  this  hypothesis.  Figure  5.16 
shows  how  changing  generator  compositions  for  the  lOkW  scenarios  impacted  fuel  con¬ 
sumption  when  GREENS  systems  were  not  included.  In  this  figure  it  is  seen  that,  when 
GREENS  systems  are  not  considered,  greater  fuel  consumption  is  seen  in  the  scenarios 
including  a  single  large  generator  than  in  any  of  the  other  three  options  simulated.  All 
scenarios  included  cumulative  generator  power  production  capacity  capable  of  providing 
at  least  peak  power  for  the  load  (lOkW). 

Having  identified  the  benefit  of  providing  a  micropower  system  with  more  options  than  just 
a  single  peak  power  capable  generator,  it  is  next  valuable  to  identify  how  inclusion  of  solar 
energy  systems  influences  this  relationship.  Eigure  5.17  shows  that  the  fuel  savings  ben¬ 
efits  of  employing  GREENS,  in  terms  of  percent  fuel  savings  for  the  micropower  system, 
differed  between  the  different  generator  compositions.  This  figure  shows  that  the  scenario 
with  a  single  lOkW  generator  benefited  the  most  from  employing  a  GREENS  system. 

The  bottom  half  of  Eigure  5.16  shows  that,  when  GREENS  systems  are  included,  the  sce¬ 
narios  with  only  one  peak  power  capable  generator  see  the  greatest  decrease  in  fuel  savings 
due  to  inclusion  of  a  PV-battery  system  in  the  micropower  system.  In  fact,  the  6  month  sce¬ 
narios  with  only  the  one  lOkW  generator  go  from  having  a  higher  median  fuel  consumption 
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Figure  5.17:  HOMER  predictions  of  fuel  savings  to  be  gained  from  employing  GREENS  systems, 
with  savings  specified  as  a  percentage  of  fuel  consumed  when  no  GREENS  are  employed.  Graph 
separates  HOMER  predictions  by  duration  of  scenarios  and  by  scenario  generator  compositions 


than  all  other  scenarios,  to  having  lower  median  fuel  consumption  than  two  of  the  other 
three  scenarios  when  GREENS  systems  are  included.  This  is  a  counter-intuitive  result,  as 
one  of  those  two  other  scenarios  has  a  lOkW  generator  in  its  assets  as  well.  Because  of  this, 
it  was  expected  that  the  "OnelOkWGeneratorAndTwoSkWGenerators"  scenario  would  be 
able  to  at  least  match  the  performance  of  the  "OnelOkWGenerator"  scenario.  Instead, 
the  presence  of  the  GREENS  systems  reversed  the  relationship  between  the  two  generator 
composition  scenarios,  with  the  single  peak  power  generator  scenario  becoming  more  fuel 
efficient  than  two  of  the  three  scenarios  including  multiple,  different  sized  generators.  This 
result  shows  that,  when  GREENS  are  present,  multiple  small  generators  can’t  be  assumed 
to  consume  less  fuel  than  a  single  peak  power  generator,  as  is  seen  when  GREENS  are  not 
present  in  the  micropower  system. 

John  Glassmire  of  the  HOMER  Eimited  Eiability  Company  (EEC)  suggested  that  a  pos¬ 
sible  explanation  for  these  results  is  the  way  controllers  are  modeled  to  direct  generator 
interaction  with  battery  banks.  A  potential  reason  the  micropower  systems  don’t  perform 
as  expected  is  that  the  controller  is  operating  without  the  "knowledge  of  the  future"  that  the 
user  has.  He  explained  that  "if  you  had  perfect  knowledge  of  the  future,  you  could  dispatch 
the  generators  in  a  way  that  maximized  the  utility  of  the  battery,  but  that’s  not  realistic  [43] ." 
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Figure  5.18:  Graph  shows  the  difference  between  fuel  savings  predictions  from  employing  HOMER 
with  TMY3  solar  irradiance  profiles  and  mean  fuel  savings  predictions  from  employing  HOMER 
with  the  50  (1961-2010)  annual  measured  NREL  solar  irradiance  profiles  as  solar  resource  inputs. 


Instead,  controllers  operate  without  an  ability  to  anticipate  the  future  load  profiles  or  solar 
energy  that  will  be  available.  This  leaves  the  controllers  unable  to  optimize  the  charging 
and  discharging  of  battery  banks  or  optimize  the  employment  of  different  sized  generators 
that  are  available. 

It  is  unknown  to  the  author  whether  smart  grid  controllers  exist  that  could  overcome  this 
shortfall.  Regardless  of  whether  any  do  exist,  it  is  valuable  to  understand  how  the  addition 
of  GREENS  PV  array  and  battery  systems  can  decrease  or  reverse  the  fuel  savings  gained 
from  employing  multiple  different  sized  generators  instead  of  a  single  peak  power  capable 
generator. 


5.3  HOMER  Fuel  Savings  Predictions  Using  TMY3  Pro¬ 
files  Compared  to  Fuel  Savings  Predictions  Using 
Measured  Solar  Irradiance  Profiles 

In  the  previous  section,  the  effect  of  inter-annual  solar  irradiance  variability  on  the  inter¬ 
annual  variability  of  HOMER  predictions  was  quantified.  This  section  will  identify  how 
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the  fuel  savings  predieted  with  the  use  of  TMY3  solar  irradianee  profiles  relate  to  the  mean 
fuel  savings  predieted  with  the  use  of  the  50  annual  solar  irradianee  profiles  from  1961 
to  2010.  These  50  measured  annual  solar  irradianee  profiles  are  the  NREL  profiles  which 
were  used  to  evaluate  the  effects  of  inter-annual  solar  irradianee  variability  as  discussed 
in  the  previous  section.  This  section  will  also  discuss  how  predictions  achieved  using  the 
TMY3  profile  solar  inputs  relate  to  predictions  achieved  using  annual  measured  profiles 
from  only  the  last  10  years  of  measured  NREL  data  (2001-2010).  This  will  provide  some 
insight  on  how  long-term  trends  in  solar  irradianee,  namely  global  dimming,  impact  the 
use  of  TMY3  profiles  as  solar  resource  inputs  for  HOMER  fuel  saving  predictions. 

Comparison  of  the  fuel  savings  estimates  generated  with  TMY3  solar  profiles  and  the  fuel 
savings  estimates  generated  with  the  1961-2010  annual  measured  solar  profiles  was  con¬ 
ducted,  and  is  shown  in  Eigure  5.18.  Eor  each  individual  scenario,  Eigure  5.18  shows  a  plot 
of  the  mean  fuel  savings  predicted  using  HOMER  across  the  50  annual  solar  irradianee 
profiles  and  a  plot  of  the  fuel  savings  predicted  using  HOMER  with  TMY3  solar  irradianee 
profiles  for  each  individual  scenario.  The  difference  between  estimates  from  TMY3  inputs 
and  estimates  from  the  50  annual  measured  profiles  was  calculated  for  each  individual  sce¬ 
nario  by  subtracting  the  latter’s  percent  fuel  savings  prediction  from  the  TMY3-generated 
percent  fuel  savings  estimate.  The  visualization  of  the  differences  in  estimates  is  catego¬ 
rized  both  by  scenario  peak  load  and  by  scenario  duration.  As  can  be  seen  in  Eigure  5.18, 
the  HOMER  fuel  savings  estimated  using  TMY3  solar  irradianee  profiles  had  a  median  of 
variation  of  less  than  half  a  percent  above  the  mean  HOMER  fuel  savings  estimated  using 
the  50  measured  annual  profiles.  This  shows  that  there  is  a  difference  of  less  than  a  half 
of  a  percent  between  TMY3-generated  fuel  savings  predictions  and  the  mean  fuel  savings 
predictions  generated  derived  from  using  the  50  measured  annual  solar  profiles. 

The  second  analysis  mentioned  at  the  beginning  of  this  section,  the  identification  of  global 
dimming  effects  on  the  use  of  TMY3  profiles  in  HOMER,  was  conducted  in  a  manner 
similar  to  this  first  analysis.  The  difference  between  TMY3 -generated  fuel  savings  esti¬ 
mates  and  estimates  derived  from  using  the  10  (2001-2010)  annual  profiles  was  calculated 
for  each  individual  scenario  by  subtracting  the  latter  mean  percent  fuel  savings  prediction 
from  the  TMY3-generated  percent  fuel  savings  estimate.  The  results  of  this  analysis,  again 
categorized  by  scenario  peak  load  and  by  scenario  duration,  can  be  seen  in  Eigure  5.19. 
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Figure  5.19:  Graph  shows  the  difference  between  fuel  savings  predictions  from  employing  HOMER 
with  TMY3  solar  irradiance  profiles  and  mean  fuel  savings  predictions  from  employing  HOMER 
with  the  10  annual  measured  NREL  solar  irradiance  profiles  for  the  years  2001-2010  as  solar 
resource  inputs. 


Figure  5.19  shows  that  the  fuel  savings  estimated  with  the  HOMER  model  using  the  most 
reeent  10  years  of  annual  solar  irradianee  profiles  were  lower  than  the  TMY3-generated 
fuel  savings  estimates  for  the  different  seenarios. 

Figure  5.20  shows  the  differenee  between  TMY3 -generated  fuel  savings  estimates  and  fuel 
savings  estimates  generated  using  the  50  years  of  solar  irradianee  profiles.  Figure  5.20 
also  shows  the  differenee  between  TMY3 -generated  fuel  savings  estimates  and  the  fuel 
savings  estimates  generated  using  the  most  reeent  10  years  of  solar  irradianee  profiles. 
In  Figure  5.20  it  ean  be  seen  that  the  median  differenee  between  HOMER  fuel  savings 
predietions  derived  from  using  2000-2010  annual  solar  irradianee  profiles  and  HOMER 
fuel  savings  predietions  derived  from  using  TMY3  profiles  is  1.36%  and  1.35%  for  the 
3kW  and  lOkW  seenarios,  respeetively.  The  median  differenee  between  HOMER  fuel 
savings  predietions  derived  from  using  TMY3  solar  profiles  and  the  HOMER  fuel  savings 
predietions  derived  from  using  all  50  years  of  solar  irradianee  profiles  is  0.37%  and  0.33% 
for  the  3kW  and  lOkW  seenarios,  respeetively.  This  suggests  that  HOMER  fuel  savings 
predietions  made  with  TMY3  solar  irradianee  input  profiles  ean  be  expeeted  to  overestimate 
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Figure  5.20:  Graph  shows  the  difference  between  fuel  savings  predictions  derived  from  employing 
HOMER  with  TMY3  solar  irradiance  profiles  and  mean  fuel  savings  predictions  from  employing 
HOMER  with  the  10  annual  measured  NREL  solar  irradiance  profiles  for  the  years  2001-2010 
as  solar  resource  inputs.  Graph  also  shows  the  difference  between  fuel  savings  predictions  from 
employing  HOMER  with  TMY3  solar  irradiance  profiles  and  mean  fuel  savings  predictions  from 
employing  HOMER  with  the  50  (1961-2010)  annual  measured  NREL  solar  irradiance  profiles  as 
solar  resource  inputs. 


fuel  savings  to  be  derived  from  using  PV-battery  systems  by  roughly  one  pereentage  point 
due  to  the  effeets  of  global  dimming  alone. 


5.4  Analysis  of  HOMER  Robustness  With  Regard  to  Inter- 
Annual  Temperature  Variability 

Analysis  of  how  inter-annual  temperature  variability  affeets  the  robustness  of  the  HOMER 
model’s  fuel  savings  predietions  was  aeeomplished  in  a  similar  manner  to  the  analysis 
demonstrated  above,  whieh  investigated  solar  irradianee  variability  effeets.  RMSE  is  the 
MOP  for  inter-annual  variability  of  pereent  fuel  savings  predieted  by  the  HOMER  to  be 
aehieved  by  employing  GREENS  systems.  This  is  the  same  MOP  for  HOMER  robust¬ 
ness  as  was  used  in  the  analysis  of  inter-annual  solar  irradianee  variability  effeets  on  the 
model.  While  CVRMSE  was  utilized  as  the  measure  of  inter-annual  solar  irradianee  vari¬ 
ability,  RMSE  is  utilized  as  the  measure  of  inter-annual  temperature  variability.  By  using 
an  average  temperature  rather  than  eumulative  temperature  measure  in  this  analysis,  the 
effeets  of  varying  the  duration  of  simulations  is  already  aeeounted  for.  This  allows  RMSE 
to  be  employed  rather  than  CVRMSE  for  identifying  and  eomparing  the  inter-annual  tem- 
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perature  variability  of  the  scenarios.  The  following  subsections  iterate  through  the  factors 
of  the  DOE,  identifying  how  inter-annual  temperature  variability  affects  the  robustness  of 
HOMER  fuel  savings  predictions  and  identifying  how  each  factor  of  the  DOE  affects  that 
relationship.  Eor  analysis  of  temperature  effects,  generator  composition  was  not  included. 
The  additional  experiment  with  scenarios  generator  compositions  consisting  of  one  lOkW 
generator  alone  was  only  executed  with  variations  in  the  solar  irradiance  input  and  not  tem¬ 
perature  input.  Without  this  scenario,  there  is  little  that  can  be  gleaned  from  the  analysis 
of  generator  composition  effects  with  regard  to  inter-annual  temperature  variability.  The 
ordering  of  the  factors  considered  (unit  size,  duration,  season,  and  climate)  is  based  on 
partition  tree  identification  of  decreasing  influence  on  inter-annual  variability  of  HOMER 
fuel  savings  predictions. 

5.4.1  Temperature  Variability:  Effects  of  Unit  Size  on  Variability  in 
HOMER  Enel  Savings  Estimates 

The  first  partition  tree  executed  in  this  analysis  found  that  the  peak  load  of  the  scenario  con¬ 
sidered  was  the  factor  with  the  greatest  influence  on  the  relationship  between  inter-annual 
temperature  variability  and  inter-annual  variability  in  HOMER  fuel  savings  predictions. 
Eigure  5.21  shows  the  first  split  of  this  initial  partition  tree.  This  split  shows  that  24.5% 
of  inter-annual  variability  in  HOMER  fuel  savings  predictions  resulting  from  inter-annual 
temperature  variability  can  be  explained  by  identifying  whether  the  peak  load  of  the  sce¬ 
nario  is  above  or  below  30kW.  It  should  be  recalled  that  the  peak  load  factor  was  similarly 
identified  as  the  most  influential  factor  on  the  effects  of  inter-annual  solar  irradiance  vari¬ 
ability  on  HOMER  fuel  savings  prediction  variability. 

As  with  the  analysis  of  solar  irradiance,  the  peak  load  factor  was  further  evaluated  through 
linear  regression  analysis  and  visual  exploration  of  inter-annual  HOMER  fuel  savings  pre¬ 
diction  variability  across  the  five  peak  load  levels.  Box  plots  were  utilized  to  identify  the 
quantity  and  range  of  inter-annual  variability  in  the  HOMER  fuel  savings  predictions  for 
each  of  the  five  load  profile  peak  load  levels.  Eigure  5.22  shows  that  the  quantity  of  inter¬ 
annual  variability  in  HOMER  fuel  savings  predictions  is  far  less  when  due  to  temperature 
variability  than  it  was  when  due  to  inter-annual  solar  irradiance  variability,  seen  in  Eig¬ 
ure  5.2.  The  median  inter-annual  variability  of  HOMER  predictions  for  the  3kW  scenarios, 
for  example,  is  approximately  1.7%  when  resulting  from  inter- annual  solar  irradiance  vari- 
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Figure  5.21:  Partition  tree  with  one  split  identifying  scenario  peak  load  level  as  the  most  influential 
factor  for  inter-annual  HOMER  fuel  savings  prediction  variability  caused  by  inter-annual  mean 
temperature  variability 


ability  and  only  approximately  0.13%  when  it  is  resulting  from  inter- annual  temperature 
variability. 

Although  the  inter- annual  variability  in  HOMER  fuel  savings  predietions  resulting  from 
temperature  variability  is  elose  to  a  magnitude  of  ten  less  than  the  variability  resulting 
from  solar  irradianee  variability,  the  influenee  of  seenario  load  profile  peak  load  levels 
on  variability  of  HOMER  predietions  is  similar.  Eigure  5.22  shows  how  the  relationship 
between  inter-annual  variability  of  HOMER  fuel  savings  predietions  and  the  different  peak 
load  levels  of  the  seenarios  is  similar  to  what  was  seen  when  assessing  the  effeets  of  solar 
irradianee  inter- annual  variability.  Here  again  the  3kW  and  lOkW  seenarios  demonstrate  a 
higher  degree  of  inter-annual  variability  in  HOMER  fuel  savings  predietions  than  the  5kW, 
30kW,  or  60kW  seenarios.  Eigure  5.23  helps  to  explain  this  relationship  between  peak 
load  levels  and  inter- annual  variability  in  HOMER  fuel  savings  predietions.  It  shows  that 
the  ratio  of  peak  load  level  to  PV  power  produetion  eapaeity  onee  again  appears  to  serve 
as  a  eonfounding  faetor  on  inter-annual  variability  in  fuel  savings  predietions.  The  table 
showing  the  peak  load  level  to  PV  power  produetion  eapaeity  ratios  ean  be  reviewed  in 
Eigure  5.3. 

Eigure  5.23  visualizes  the  eonfounding  effeet  of  these  ratios  with  the  use  of  linear  regres¬ 
sion  lines,  as  was  done  in  the  analysis  of  solar  irradianee  variability  effeets.  Similar  to  what 
was  seen  in  this  previous  solar  irradianee  analysis,  when  eonsidering  the  effeets  of  inter- 


84 


Figure  5.22:  Graphical  representation  of  the  impact  of  scenario  peak  load  on  variability  in  HOMER 
percent  fuel  savings  predictions  (RMSE) 

annual  temperature  variability  the  linear  relationship  between  mean  predieted  fuel  savings 
and  variability  in  predieted  fuel  savings  is  more  prevalent  in  the  5kW,  30kW,  and  60kW 
seenarios.  Figure  5.23  shows  the  linear  regressions  for  the  five  peak  load  levels,  with  the 
5kW,  30kW  and  60kW  seenarios  having  RSquare  values  of  0.401,  0.569  and  0.775,  respee- 
tively.  The  3kW  and  lOkW  seenarios,  as  seen  in  Figure  5.23,  demonstrate  a  weaker  linear 
relationship  between  mean  predieted  fuel  savings  and  inter-annual  variability  in  HOMER 
fuel  savings  predietions,  with  RSquare  values  of  0.112  and  0.199,  respeetively.  The  dif- 
ferenee  seen  in  HOMER  predietion  variability  between  the  two  groups  of  peak  load  levels 
(3kW  and  lOkW  versus  5kW,  30kW,  and  60kW)  is  assumed  here  to  be  attributed  to  the 
different  ratios  of  peak  load  to  PV  energy  produetion  eapabilities  between  the  two  groups. 

Eor  eonsideration  of  the  remaining  faetors  in  the  analysis  of  variability  in  HOMER  fuel 
savings  estimates  due  to  inter-annual  temperature  variability  all  five  peak  load  levels  are 
again  ineluded  in  the  analysis.  Also,  the  3kW  and  lOkW  peak  load  seenarios  are  again 
eonsidered  to  be  more  representative  of  the  effeets  of  temperature  inter-annual  variabil¬ 
ity  on  inter- annual  variability  of  HOMER  fuel  savings  predietions.  This  is  beeause  of  the 
higher  ratio  of  PV  power  produetion  eapaeity  to  system  peak  load  levels  in  the  3kW  and 
lOkW  seenarios  and  the  subsequently  weaker  linear  relationship  between  HOMER’s  mean 
pereent  fuel  savings  predietions  and  inter-annual  variability  in  HOMER  pereent  fuel  sav- 
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Figure  5.23:  Graph  shows  strong  linear  relationship  between  mean  predicted  fuel  savings  and 
variability  in  predicted  fuel  savings  for  5kW,  30kW  and  60kW  scenarios  for  HOMER  prediction 
variability  resulting  from  inter-annual  scenario  temperature  variability 

ings  predictions. 

5.4.2  Temperature  Variability:  Effects  of  Duration  on  Variability  in 
HOMER  Euel  Savings  Estimates 

The  factor  with  the  second  greatest  influence  on  HOMER  fuel  savings  predictions  inter¬ 
annual  variability  due  to  inter-annual  temperature  variability  was  identified  through  the  use 
of  a  partition  tree.  This  partition  tree  did  not  include  peak  load  as  a  factor  and  instead 
considered  the  three  factors  location,  start  time  (season),  and  duration  as  x  variables.  The 
RMSE  variability  of  HOMER  percent  fuel  savings  predictions  was  set  as  the  y  variable 
for  the  partition  tree.  This  partition  tree,  seen  in  Eigure  5.24,  shows  that  the  duration  of 
the  scenario  is  the  second  most  influential  factor  concerning  how  inter-annual  temperature 
variability  affects  inter-annual  variability  of  HOMER  fuel  savings  estimates.  The  partition 
tree  shows  that  13.4%  of  inter-annual  variability  in  HOMER  fuel  savings  predictions  can  be 
explained  by  identifying  whether  a  given  scenario’s  duration  is  greater  than  or  less  than  two 
months.  This  is  of  course  given  that  the  variability  in  HOMER  fuel  savings  predictions  are 
due  to  inter-annual  temperature  variability,  with  the  solar  irradiance  resource  kept  constant 
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Figure  5.24:  Partition  tree  showing  duration  of  the  scenario  to  be  the  second  most  influential 
factor  on  the  inter-annual  variability  of  HOMER  fuel  savings  predictions  when  that  variability  is 
due  to  inter-annual  temperature  variability 


for  all  years  using  the  NREL  TMY3  solar  irradiance  profiles. 

To  better  quantify  the  effects  of  duration  as  a  factor  across  each  of  the  five  peak  load  levels, 
five  partition  trees  were  reviewed,  one  for  each  peak  load  level.  In  Figure  5.25  the  first  split 
for  these  partition  trees  shows  that  duration  proved  to  be  the  second  most  influential  factor 
for  four  of  the  five  peak  load  levels.  The  duration  factor’s  influence  is  greatest  for  the  3kW 
and  lOkW  scenarios,  with  RSquare  values  of  0.443  and  0.354,  respectively.  That  being 
said,  even  the  5kW  and  30kW  scenarios  show  duration  to  be  their  second  most  influential 
factor  with  RSquare  values  of  0.23 1  and  0. 1 69,  respectively.  Only  the  60kW  scenarios  show 
scenario  start  time  to  be  the  second  most  influential  factor,  with  an  RSquare  of  0.051,  for 
explaining  inter-annual  variability  of  HOMER  fuel  savings  predictions  when  variability  is 
due  to  inter-annual  temperature  variability.  Even  for  the  60kW  scenario,  however,  duration 
can  be  seen  as  an  influential  factor  on  the  fourth  split  of  the  partition  tree,  when  this  factor 
increases  the  RSquare  value  from  0.261  to  0.501.  It  is  notable  that  for  the  all  of  the  other 
four  peak  load  level  scenarios  the  important  split  occurs  at  the  same  level  for  the  duration 
factor;  the  two  month  mark.  As  can  be  seen  in  Figure  5.25,  there  is  a  noticeable  increase 
in  both  the  average  variability  and  range  of  variability  in  fuel  savings  predictions  when 
duration  values  are  less  than  two  months,  or  1440  hours. 

Now  that  duration  has  been  identified  as  the  second  most  influential  factor  for  explaining 
the  relationship  between  inter-annual  temperature  variability  and  inter-annual  variability  in 
HOMER  fuel  savings  predictions,  a  visual  exploration  of  the  relationship  will  support  a 
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Figure  5.25:  The  first  split  in  partition  trees  associated  with  each  peak  load  level,  identifying 
duration  as  the  second  most  influential  factor  for  inter-annual  variability  in  HOMER  fuel  savings 
predictions  when  variability  is  due  to  inter-annual  temperature  variability 


qualification  of  the  relationship  between  HOMER  variability  and  seenario  duration.  Fig¬ 
ure  5.26  supports  sueh  a  visual  exploration  of  the  relationship  between  duration  and  inter¬ 
annual  variability  of  HOMER  model  fuel  savings  predietions  when  this  variability  is  due 
to  inter-annual  temperature  variability.  This  visualization  shows  that  inter-annual  variabil¬ 
ity  of  HOMER  fuel  savings  predietions  inereases  at  a  rate  that  inereases  substantially  as 
duration  of  the  seenario  deereases  below  three  months.  The  relationship  between  duration 
and  variability  predictions  is  the  greatest  for  the  scenarios  with  3kW  and  lOkW  peak  load 
levels;  this  was  antieipated,  given  the  output  of  the  previously  diseussed  partition  trees. 

For  the  3kW  and  lOkW  peak  load  levels,  approximately  75%  of  fifteen  day  seenarios  show 
greater  inter-annual  variability  than  50%  of  one  month  seenarios.  For  these  same  peak 
load  levels,  50%  of  one  month  seenarios  show  greater  inter- annual  variability  than  75% 
of  all  seenarios  with  larger  duration  levels.  While  the  relationship  is  more  prominent  for 
seenarios  with  3kW  and  lOkW  peak  load  levels,  the  relationship  seen  between  seenario 
duration  and  inter-annual  variability  in  HOMER  fuel  savings  predietions  is  seen  aeross  all 
peak  load  levels. 
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Figure  5.26:  Visualization  of  the  relationship  between  duration  and  inter-annual  variability  of 
HOMER  model  fuel  savings  predictions  when  variability  is  due  to  inter-annual  temperature  vari¬ 
ability 
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Figure  5.27:  The  first  split  in  the  partition  tree  identifying  start  time  as  the  third  most  influential 
factor  for  inter-annual  variability  in  HOMER  fuel  savings  predictions  when  variability  is  due  to 
inter-annual  temperature  variability 


5.4.3  Temperature  Variability:  Seasonal  Effects  on  Variability  in 
HOMER  Enel  Savings  Estimates 

When  duration  and  peak  load  are  not  ineluded  as  faetors  in  the  partition  tree  analysis,  the 
start  time  of  seenarios  is  seen  to  be  the  third  most  influential  faetor  for  understanding  inter¬ 
annual  variability  in  HOMER  fuel  savings  predietions.  Figure  5.27  shows  this  partition 
tree,  with  the  first  split,  about  the  start  time  of  September  1st,  yielding  an  RSquare  value  of 
0.045. 
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Visual  exploration  of  this  factor’s  influence  on  inter-annual  variability  of  HOMER  fuel 
savings  predictions  showed  that  the  relationship  between  scenario  start  time  and  the  vari¬ 
ability  in  inter-annual  HOMER  fuel  savings  predictions  was  related  to  the  relationship  be¬ 
tween  scenario  start  time  and  mean  fuel  savings  predictions.  This  relationship  is  similar  to 
that  seen  when  assessing  the  relationship  between  scenario  start  time  and  HOMER  predic¬ 
tion  variability  due  to  solar  irradiance  variability.  Once  again,  the  confounding  effects  of 
the  seasonal  fluctuations  in  mean  load  profiles  appear  to  explain  more  of  the  inter-annual 
variability  in  HOMER  fuel  savings  predictions  than  does  the  inter-annual  temperature  vari¬ 
ability  associated  with  the  seasons.  While  the  relationship  between  scenario  start  time  and 
scenario  mean  temperature  is  unimodal,  the  graphical  representation  of  start-time  relation¬ 
ship  with  HOMER  prediction  variability  is  bi-modal,  with  the  latter  matching  the  bi-modal 
relationship  between  mean  fuel  savings  and  scenario  start-time.  Through  this  analysis, 
the  relationship  between  scenario  season  and  the  inter-annual  variability  of  HOMER  fuel 
savings  predictions  is  considered  to  be  minimal  for  the  purposes  of  this  research.  An  im¬ 
portant  caveat,  as  stated  with  the  solar  irradiance  analysis,  is  that  seasonal  effects  seen  in 
climates  outside  the  continental  United  States  could  potentially  have  greater  influence  on 
inter- annual  variability  of  HOMER  simulations. 


5.4.4  Temperature  Variability:  Climate  Effects  on  Variability  in  HOMER 
Enel  Savings  Estimates 

The  next  factor  considered  is  the  climate  of  the  scenarios.  Eigure  5.28  shows  that  climate 
has  little  influence  on  the  inter-annual  variability  of  HOMER  fuel  savings  predictions  when 
that  variability  is  due  to  inter-annual  temperature  variability.  While  the  steppe  climate 
provides  for  greater  mean  predicted  fuel  savings  percentages  than  the  other  two  climates 
throughout  the  year,  there  is  no  one  climate  that  shows  a  consistently  higher  or  lower  inter¬ 
annual  variability  in  fuel  savings  predictions  relative  to  the  other  climates.  This  suggests 
that,  when  considering  the  effects  of  inter-annual  temperature  variability,  the  scenario’s 
climate  demonstrates  minimal  influence  on  the  inter-annual  variability  of  HOMER  fuel 
savings  estimates. 


90 


Figure  5.28:  Effects  of  climate  on  inter-annual  variability  of  HOMER  fuel  savings  predictions  due 
to  temperature  variability,  shown  across  the  seasons  for  lOkW  and  15-day  scenarios 


5.5  Comparing  the  Effects  of  Solar  Irradiance  and  Tem¬ 
perature  Inter-annual  Variability  on  the  Robustness 
of  HOMER  Predictions 

When  the  effeets  of  inter-annual  temperature  variability  are  eonsidered,  variability  of 
HOMER  fuel  savings  predietions  for  individual  peak  load  and  duration  faetor  levels  is 
less  than  the  variability  for  equivalent  seenarios  assoeiated  with  only  inter-annual  solar  ir- 
radianee  variability.  While  the  relationship  between  duration  and  inter-annual  HOMER 
fuel  savings  predietions  is  signifieant  for  the  temperature  variability  seenarios,  the  quantity 
of  inter-annual  variability  is  smaller  than  it  is  when  due  to  solar  irradianee  variability.  Eig- 
ure  5.29  shows  how  solar  irradianee  inter-annual  variability  has  a  mueh  greater  influenee 
on  the  inter-annual  variability  of  HOMER  fuel  savings  predietions  than  does  temperature 
variability.  Eor  example,  the  HOMER  fuel  savings  predietions  for  3kW,  15-day  seenarios 
have  a  median  inter-annual  variability  of  approximately  0.275%  when  the  variability  is  due 
to  inter- annual  temperature  variability.  The  3kW,  15-day  seenarios  have  a  median  inter¬ 
annual  variability  of  approximately  2.2%  when  the  variability  is  due  to  inter-annual  solar 
irradianee  variability. 
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Figure  5.29:  Visualization  of  how  inter-annual  variability  (RMSE)  of  percent  fuel  savings  differ 
between  HOMER  simulations  that  utilize  measured  temperature  profiles  and  measured  annual 
solar  irradiance  profiles  ("R"  graph)  and  HOMER  simulations  that  utilize  measured  temperature 
profiles  and  TMY3  solar  irradiance  profiles  ("TMY"  graph) 


Interestingly,  when  both  solar  and  temperature  inter-annual  variability  are  considered  in 
the  HOMER  simulations,  the  inter-annual  variability  of  HOMER  fuel  savings  predictions 
is  actually  less  than  the  variability  resulting  from  solar  irradiance  variability  alone.  This  de¬ 
crease  in  inter- annual  HOMER  prediction  variability  is  seen  in  Eigure  5.30.  The  decrease 
in  variability  can  be  explained  by  understanding  the  relationship  between  the  correlated 
solar  irradiance  and  temperature  input  files.  Eor  a  given  location  solar  irradiance  and  tem¬ 
perature  values  tend  to  fluctuate  with  the  seasons  in  an  inverse  relationship  relative  to  each 
other.  In  the  summer,  for  example,  temperature  values  tend  to  be  at  their  highest  while  solar 
irradiance  values  tend  to  be  at  their  lowest  values.  In  the  winter  season,  temperature  values 
tend  to  be  at  their  lowest  while  solar  irradiance  values  tend  to  be  at  their  highest.  This 
inverse  relationship  between  temperature  and  solar  irradiance  fluctuations  across  the  sea¬ 
sons  explains  why  the  combination  of  solar  irradiance  and  temperature  data  in  a  HOMER 
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Figure  5.30:  Visualization  of  how  inter-annual  variability  (RMSE)  of  percent  fuel  savings  differ 
between  HOMER  simulations  that  utilize  measured  temperature  profiles  and  measured  annual 
solar  irradiance  profiles  ("Temp  Only"),  HOMER  simulations  that  utilize  measured  temperature 
profiles  and  TMY3  solar  irradiance  profiles  ("TempAndSolar"),  and  HOMER  simulations  that 
utilize  measured  annual  solar  irradiance  files  and  don't  consider  temperature  effects  ("SolarOnly") 


simulation  would  decrease  inter- annual  variability  of  the  predictions.  Increases  in  power 
generation  that  would  be  associated  with  increases  in  solar  irradiance  beyond  the  mean 
level  would  be  partially  dampened  by  the  associated  decreases  in  ambient  temperature  for 
the  same  period  of  time.  This  could  explain  why  including  both  solar  irradiance  input  pro¬ 
files  and  correlated  temperature  input  profiles  in  HOMER  simulations  results  in  a  decrease 
in  the  inter-annual  variability  of  HOMER  simulations. 


5.6  Using  HOMER  to  Assess  Grid-tie  versus  Standalone 
Converters  for  PV  Systems 

All  of  the  experiments  conducted  for  this  research  included  Marine  Corps  AMMPS  gen¬ 
erators  and  a  slightly  modified  representation  of  the  Marine  Corps  GREENS  renewable 
power  system.  The  modification  of  the  GREENS  specifications  was  the  representation  of 
a  single  GREENS  as  being  able  to  produce  1.76kW  of  power  with  its  1.76kW  PV  arrays, 
rather  than  limiting  the  power  output  to  IkW  as  is  seen  with  the  real  GREENS  system. 
The  GREENS  specifications  were  modeled  in  this  way  in  the  HOMER  model  because  the 
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Mean  Fuel  Saved  By  Using  PV  Arrays  and  Batteries  With  Existing  Converter 
(Liters)  vs.  Duration  of  Simulation 
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Figure  5.31:  Mean  fuel  savings  derived  from  employing  the  GREENS  PV  arrays  and  batteries 
with  the  existing  standalone  converter  and  without  the  controller  decreasing  power  provided  by 
the  PV  arrays 


model  was  unable  to  represent  the  GREENS  eontroller’s  limitation  on  power  produetion 
without  deereasing  the  modeled  size  of  the  PV  arrays  themselves.  It  is  important  to  note 
that  the  HOMER  EEC,  being  made  aware  of  this  differenee  between  the  existing  model  and 
the  GREENS  capabilities,  will  be  able  to  correct  the  HOMER  representation  of  GREENS 
in  future  versions.  Although  this  representation  of  GREENS  facilitated  the  analysis  of 
HOMER  robustness  concerning  variability  of  inputs,  it  also  facilitated  an  assessment  of 
how  PV  arrays  and  batteries  equivalent  in  size  to  the  GREENS  components  could  support 
Marine  Corps  micropower  systems  if  employed  in  a  way  different  from  the  current  method. 
This  analysis  has  demonstrated  the  notable  fuel  savings  that  can  be  derived  from  employing 
PV  arrays  and  battery  banks  in  direct  support  of  loads  in  conjunction  with  other  microp¬ 
ower  system  assets.  The  mean  fuel  savings  experimentally  calculated  across  the  scenarios 
are  shown  in  Eigure  5.31,  the  mean  values  are  organized  by  quantity  of  GREENS  and  by 
the  duration  of  the  simulation.  This  research  went  a  step  further  in  the  analysis  of  Marine 
Corps  renewable  energy  power  planning  by  assessing  how  the  implementation  of  grid-tie 
converters  would  affect  the  fuel  savings  to  be  derived  from  employing  renewable  power 
sources. 
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Figure  5.32:  Comparison  of  mean  fuel  savings  derived  from  employing  the  GREENS  PV  arrays 
and  batteries  with  a  grid-tie  capable  versus  using  a  standalone  converter 


Following  analysis  of  inter-annual  variability  in  temperature  and  solar  irradianee  inputs 
for  existing  eonverter  type,  this  researeh  shifted  to  an  assessment  of  HOMER  predietions 
of  grid-tie  eapable  PV  systems  for  the  same  seenarios.  This  analysis  of  eonverter  effeets 
on  the  fuel  savings  to  be  derived  from  GREENS  systems  was  eondueted  only  with  mea¬ 
sured  annual  solar  irradianee  profiles  as  the  solar  resouree  input  and  with  the  effeets  of 
temperature  not  ineluded  in  the  model.  The  only  ehange  made  to  the  PV  array  systems 
was  the  speeifieation  of  the  inverters  as  being  grid-tie  eapable.  This  speeifieation  ehange 
was  exeeuted  to  identify  the  potential  fuel  savings  of  the  GREENS  PV  arrays  if  the  eon¬ 
verter  was  ehanged  to  be  a  grid-tie  eapable  eonverter.  This  modifieation  to  the  GREENS 
speeifieations  is  in  addition  to  the  previously  explained  modifieation  of  not  limiting  the  PV 
array  power  produetion  from  1.76kW  to  IkW  per  individual  GREENS  system.  Figure  5.32 
shows  the  mean  fuel  savings  predieted  through  the  use  of  HOMER  for  the  seenarios  when 
the  GREENS  eonverter  is  made  to  be  grid-tie  eapable  eompared  to  the  mean  fuel  savings 
predieted  for  seenarios  utilizing  the  existing  stand-alone  eonverter  speeifieation. 

Eigure  5.32  demonstrates  an  interesting  relationship  between  the  peak  load  size  and  the 
benefit  of  employing  a  grid-tie  eapable  eonverter.  While  having  a  grid-tie  eapability  showed 
to  be  of  little  value  for  the  smaller  mieropower  systems  (even  slightly  detrimental  at  times) 
the  use  of  grid-tie  eonverters  resulted  in  signifieant  fuel  savings  over  the  use  of  standalone 


95 


converters  when  the  peak  load  level  is  higher. 


As  the  Marine  Corps  pushes  forward  in  its  development  and  fielding  of  renewable  energy 
technologies,  this  analysis  of  grid-tie  eapable  PV  array  performanee  provides  valuable  in¬ 
sight  for  critieally  assessing  the  Marine  Corps’  approach  to  renewable  power  employment. 
It  ean  help  answer  the  question  of  whether  stand-alone  systems  are  still  the  best  way  to  go, 
or  if  fuel  savings  derived  from  grid-tie  capable  PV  systems  would  be  worth  the  investment 
of  time  and  resources.  These  results  show  that  the  pursuit  of  grid-tie  eapabilities  for  PV- 
battery  bank  systems  would  likely  only  inerease  fuel  savings  for  mieropower  systems  with 
peak  load  levels  of  30kW  or  greater. 
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CHAPTER  6: 

Conclusions,  Recommendations,  and  Future  Research 


This  research  included  analysis  of  the  robustness  of  the  HOMER  model  with  regard  to  solar 
irradiance  and  temperature  variability  and  quantified  how  several  factors  relate  to  that  ro¬ 
bustness.  This  research  also  included  an  analysis  of  global  dimming  and  its  effects  on  solar 
power  generation,  the  benefits  of  employing  a  PV  array  and  battery  bank  system  similar  to 
GREENS  in  direct  support  of  a  micropower  system,  and  the  potential  fuel  saving  benefits 
of  making  such  a  GREENS-like  system  grid-tie  capable.  In  this  chapter,  these  findings  will 
be  briefly  discussed  and  recommendations  will  be  provided  for  improving  Marine  Corps 
use  of  the  HOMER  model,  PV  arrays  and  battery  banks,  and  power  planning  in  general. 
Recommendations  for  future  work  are  also  provided  to  illustrate  additional  research  that  is 
needed  to  further  support  these  findings  and  to  investigate  additional  questions  highlighted 
by  this  research. 

The  variability  in  HOMER-generated  fuel  savings  estimates,  from  including  PV-battery 
bank  systems,  associated  with  the  inter- annual  variability  of  solar  irradiance  and  temper¬ 
ature  profiles  was  found  to  be  statistically  significant,  but  relatively  minor.  The  highest 
degree  of  variability  seen  in  HOMER  fuel  savings  estimates  was  seen  with  the  15-day, 
lOkW  scenarios.  When  considering  solar  irradiance  variability,  these  simulations,  with 
mean  percent  fuel  savings  estimates  of  33.8%,  demonstrated  a  mean  variability  in  per¬ 
cent  fuel  savings  estimates  of  2.2%,  with  a  95%  confidence  interval  of  1.2  to  3.4%.  This 
means  that  for  a  HOMER  estimate  of  33%  GREENS -derived  fuel  savings,  a  consideration 
of  inter-annual  solar  irradiance  variability  results  in  anticipating  the  fuel  savings  estimate 
to  lie  between  30.8  and  35.2%.  A  more  conservative  consideration  of  variability,  using 
the  3.4%  variability  level  seen  as  the  highest  value  within  the  95%  confidence  interval, 
would  result  in  a  fuel  savings  estimate  between  29.6  and  36.4%.  As  was  discussed  ear¬ 
lier,  scenarios  with  only  slightly  longer  durations  and  larger  load  to  PV  power  capacity 
ratios  demonstrated  far  less  variability.  This  shows  that  there  is  promise  for  Marine  Corps 
implementation  of  the  HOMER  model  in  support  of  expeditionary  power  planning. 

This  analysis  of  the  variability  associated  with  the  HOMER  model’s  stochastic  inputs  has 
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quantified  how  changes  to  certain  simulation  factors,  such  as  duration,  influence  the  confi¬ 
dence  that  users  can  hold  when  using  HOMER  fuel  consumption  estimates.  The  following 
sections  provide  a  summary  of  the  analysis  as  well  as  some  additional  insights  gleaned  from 
this  research.  The  findings  and  recommendations  below  provide  information  to  improve 
Marine  Corps  power  planning  with  the  HOMER  model  and  renewable  energy  technology. 

6.1  Evaluation  of  HOMER  Model  Robustness  with  Re¬ 
gard  to  Inter-Annual  Variability  in  Solar  Irradiance 
and  Temperature  Profiles 

This  research  quantified  the  variability  of  HOMER  fuel  savings  predictions  that  is  due  to 
inter- annual  solar  irradiance  and  temperature  variability  across  multiple  factors.  It  was  first 
identified  that,  of  the  five  factors  assessed,  the  peak  load  level  (including  ratio  of  peak  load 
level  to  PV  power  production  capacity)  and  the  duration  of  a  simulated  operation  are  the 
two  most  useful  factors  for  explaining  the  variability  of  HOMER  fuel  savings  predictions. 
Eigure  6.1  and  Eigure  6.2  quantify  the  degrees  of  variability  seen  across  the  levels  of  these 
two  factors  for  the  50  years  of  measured  data.  Both  Eigure  6.1  and  Eigure  6.2  illustrate 
how  decreases  in  duration  and  increases  in  the  peak  load  to  PV  production  capacity  ratio 
increase  the  variability  of  a  scenario’s  HOMER  fuel  savings  predictions.  This  can  prove 
useful  for  power  planners,  as  the  tables  quantify  how  changes  in  operation  parameters 
influence  the  variability  of  HOMER  fuel  savings  estimates. 

Comparing  Eigure  6.1  and  Eigure  6.2  shows  that  inter-annual  solar  irradiance  variability 
has  a  greater  impact  on  the  variability  of  HOMER-generated  fuel  savings  predictions  than 
does  the  inter-annual  variability  of  temperature  input  profiles.  Variability  in  HOMER  fuel 
savings  estimates  that  is  due  to  inter-annual  solar  irradiance  variability  is  roughly  a  power 
of  10  greater  than  variability  due  to  inter- annual  temperature  variability. 

The  other  three  factors  included  in  the  DOE,  simulation  start  time,  location/climate,  and 
composition  of  diesel  generators,  demonstrated  weaker  relationships  with  the  variability  of 
HOMER  fuel  savings  predictions.  The  relationships  between  these  factors  and  the  variabil¬ 
ity  in  HOMER  fuel  savings  estimates  can  be  reviewed  in  sections  5.2.3,  5.2.4,  and  5.3.5  for 
solar  irradiance  variability  and  5.4.3  and  5.4.4  for  temperature  variability.  These  relation- 
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Variability  of  Fuel  Savings  Predicted  in  RMSE  from  the  Mean  Percent  Fuel  Savings  Prediction 
(when  variability  is  due  to  inter-annual  solar  irradiance  variability  across  50  years) 

Duration  of  Simulated  Operations 

15  Days 

1  Month 

2  Months 

3  Months 

4  Months 

5  Months 

6  Months 

Ratio  of  Scenario  Peak  Load  Level  to  PV  Power 
Production  Capacity  (Peak  Load  Level) 

0.33 

{3kW) 

97.5%  Quantile 

3.10% 

2.31% 

1.78% 

1.51% 

1.35% 

1.19% 

1.11% 

50%  Quantile 

2.17% 

1.61% 

1.29% 

1.10% 

1.00% 

0.92% 

0.89% 

2.5%  Quantile 

1.17% 

0.90% 

0.71% 

0.66% 

0.66% 

0.68% 

0.65% 

0.30 

(lOkW) 

97.5%  Quantile 

3.41% 

2.61% 

1.96% 

1.61% 

1.41% 

1.27% 

1.17% 

50%  Quantile 

2.21% 

1.69% 

1.31% 

1.14% 

1.04% 

0.95% 

0.90% 

2.5%  Quantile 

1.20% 

0.95% 

0.82% 

0.77% 

0.73% 

0.73% 

0.70% 

0.20 

(SkW) 

97.5%  Quantile 

2.91% 

2.20% 

1.57% 

1.28% 

1.06% 

0.86% 

0.78% 

50%  Quantile 

1.59% 

1.23% 

0.96% 

0.79% 

0.72% 

0.68% 

0.63% 

2.5%  Quantile 

0.00% 

0.09% 

0.21% 

0.37% 

0.44% 

0.47% 

0.46% 

0.16 

{30kW) 

97,5%  Quantile 

2.87% 

2.12% 

1.47% 

1.13% 

0.89% 

0.68% 

0.60% 

50%  Quantile 

1.34% 

0.98% 

0.75% 

0.61% 

0.54% 

0.51% 

0.48% 

2.5%  Quantile 

0.00% 

0.19% 

0.32% 

0.32% 

0.18 

(eokw) 

97.5%  Quantile 

1.62% 

1.08% 

0.59% 

0.36% 

0.25% 

0.18% 

0.15% 

50%  Quantile 

0.00% 

0.00% 

0.13% 

0.16% 

0.15% 

0.13% 

0.11% 

2.5%  Quantile 

0% 

0% 

0% 

0% 

0.01% 

0.06% 

0.08% 

Figure  6.1:  Quantifying  the  variability  in  HOMER  Fuel  Savings  predictions  associated  with  dif¬ 
ferent  duration  and  peak  load  levels,  when  the  variability  is  due  to  inter-annual  solar  irradiance 
variability 


ships  provide  less  value  than  the  peak  load  level  and  duration  factors  for  quantifying  the 
variability  of  HOMER  fuel  savings  estimates. 

6.1.1  Recommendation:  Use  Quantification  of  Variability  for  TMY3- 
Generated  Vice  Graham-Hollands-Generated  Solar  Irradiance 
Profiles 

It  is  recommended  that  the  Marine  Corps  use  this  quantification  of  HOMER  output  vari¬ 
ability  when  using  the  NREE  TMY3  solar  irradiance  profiles  and  subsequent  versions  of 
TMY  profiles  only.  While  this  thesis  quantifies  the  variability  of  fuel  savings  estimates 
due  to  inter-annual  solar  irradiance  and  temperature  variability,  these  measures  can  only  be 
marginally  useful  when  using  Graham-Hollands-generated  solar  irradiance  profiles  without 
additional  assessment  of  the  relationship  between  TMY  profiles  and  Graham-Hollands- 
generated  profiles.  This  is  because  the  Graham-Hollands-generated  profiles  are  not  built 
using  historical,  measured  solar  irradiance  profiles  like  TMY3  profiles  are,  making  it  diffi¬ 
cult  to  identify  the  variability  of  Graham-Hollands-generated  profiles  relative  to  the  typical 
profiles  seen  with  NREE’s  TMY3  solar  irradiance  profiles.  Additionally,  based  on  the  lim¬ 
ited  evaluation  of  HOMER  output  variability  due  to  Graham-Hollands  algorithms,  TMY3 
solar  irradiance  profiles  are  suggested  to  be  preferable  to  the  Graham-Hollands  generated 
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Variability  of  Fuel  Savings  Predicted  in  RMSE  from  the  Mean  Percent  Fuel  Savings  Prediction 

(when  variability  is  due  to  inter-annual  temperature  profile  variability  across  50  years) 
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Figure  6.2:  Quantifying  the  variability  in  HOMER  Fuel  Savings  predictions  associated  with  differ¬ 
ent  duration  and  peak  load  levels,  when  the  variability  is  due  to  inter-annual  temperature  profile 
variability 


solar  irradiance  profiles  when  employing  the  HOMER  model.  This  research  has  quantified 
the  inter-annual  variability  of  HOMER  results  due  to  measured  inter-annual  solar  irradi¬ 
ance  variability  and  has  quantified  the  difference  between  fuel  savings  predictions  from 
using  TMY3  solar  irradiance  and  the  measured  annual  profiles  across  50-year  and  10-year 
periods.  The  author  knows  of  no  such  analysis  of  variability  in  HOMER  fuels  savings  pre¬ 
dictions,  and  global  dimming  effects,  associated  with  Graham-Hollands  synthetic  profiles 
or  of  any  comparison  of  Graham-Hollands  generated  fuel  savings  predictions  and  mean 
fuel  savings  predictions  for  multiple  measured  annual  profiles. 


6.1.2  Recommendation:  Develop  and  Employ  Typical  Meteorological 
Year  Temperature  Profiles 

It  is  recommended  that  a  library  of  synthetic  temperature  profiles  be  developed  or  acquired 
for  locations  for  which  the  HOMER  model  is  likely  to  be  used  by  the  Marine  Corps.  In 
section  5.5  it  is  shown  that  the  variability  of  HOMER  simulations  associated  with  inter¬ 
annual  temperature  profile  variability  is  a  power  of  ten  less  than  the  variability  associated 
with  solar  irradiance  profile  variability.  It  is  assumed  in  this  research,  from  the  example 
provided  by  Newell’s  thesis  and  the  explanation  of  NOCT  coefficient  effects  in  section 
2.3.3,  that  inclusion  of  temperature  averages  or  profiles  in  the  HOMER  model  is  important 
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for  maintaining  the  fidelity  of  the  simulations  at  an  acceptable  level.  This  suggests  that 
while  a  synthetic  temperature  profile  would  increase  simulation  fidelity  it  would  not  notably 
increase  the  variability  in  HOMER  fuel  savings  estimates.  In  fact,  inclusion  of  correlated 
temperature  profiles  could  even  decrease  variability  in  HOMER-generated  estimates. 

Section  5.5  showed  that  utilization  of  correlated  annual  solar  irradiance  and  temperature 
profiles  for  a  given  scenario  provides  for  a  slight  decrease  in  the  variability  of  HOMER 
simulation  output.  That  being  said,  it  is  unknown  if  this  effect  will  be  seen  when  TMY 
solar  and  temperature  input  files  are  utilized  rather  than  correlated  solar  irradiance  and 
temperature  input  files  for  individual  years.  It  may  also  be  useful  for  additional  research  to 
identify  whether  average  temperature  profiles  would  be  sufficient  for  providing  increased 
fidelity  with  minimal  variability,  similar  to  the  results  seen  from  employing  measured  an¬ 
nual  temperature  profiles. 


6.2  HOMER  Robustness  with  Regard  to  Global  Dimming 

This  research  used  RMSE  measures  of  variability  to  discount  the  effects  of  global  dimming 
on  the  robustness  of  HOMER  with  regard  to  inter-annual  solar  irradiance  variability.  The 
benefit  of  excluding  the  effects  of  global  dimming  in  this  analysis  of  inter-annual  solar 
irradiance  and  HOMER  fuel  savings  estimate  variability  is  that  this  analysis  shows  the 
robustness  of  the  HOMER  model  if  the  user  considers  global  dimming  in  their  employment 
of  the  model.  Had  this  research  ignored  the  effects  of  global  dimming  and  instead  measured 
inter-annual  variability  in  terms  of  CV,  the  HOMER  simulations  would  likely  have  shown  to 
be  much  less  robust  than  they  were  shown  to  be  when  global  dimming  effects  are  accounted 
for.  This  is  because  the  more  recent  measured  annual  solar  irradiance  profiles  demonstrate 
significantly  less  cumulative  surface  solar  irradiance  and  GREENS-derived  fuel  savings 
than  do  HOMER  simulations  using  measured  solar  irradiance  profiles  from  several  decades 
ago. 

The  analysis  shown  in  section  5.3  quantified  the  effects  of  global  dimming  on  HOMER 
fuel  savings  estimates.  HOMER  simulations  conducted  using  TMY3  solar  irradiance  input 
profiles  demonstrated  a  mean  over-estimate  of  less  than  half  a  percent  in  fuel  savings  when 
compared  to  mean  results  of  simulations  utilizing  annual  profiles  from  1961-2010.  The 
over-estimate  in  fuel  savings  estimates  increased  by  approximately  1  percent  for  the  same 
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scenarios  when  TMY3  simulations  were  compared  to  mean  results  of  simulations  using 
annual  profiles  from  only  2000-2010.  This  shows  that  a  failure  to  consider  the  effects  of 
global  dimming  when  conducting  HOMER  simulations  with  TMY3  profiles  can  result  in 
over-estimates  of  fuel  savings. 

6.2.1  Recommendation:  Include  Global  Dimming/Brightening  Ef¬ 
fects  in  HOMER  Modeling  or  Output  Analysis 

It  is  recommended  that  Marine  Corps  use  recent  TMY  solar  irradiance  profiles  as  resource 
inputs  for  the  HOMER  model  along  with  an  estimate  of  global  dimming  and  brightening 
effects.  It  was  shown  in  this  research  that  comparisons  can  be  made  between  simulations 
using  the  full  range  of  annual  measured  profiles  and  simulations  using  only  the  most  recent 
annual  measured  profiles  to  identify  the  degree  to  which  global  dimming  is  not  reflected 
in  TMY3  profiles.  Global  dimming  effects  could  be  replicated  by  modifying  the  solar 
irradiance  profiles,  or  the  user  could  simply  be  informed  of  the  effects  of  global  dimming 
as  a  consideration  for  their  decision  making. 

The  influence  of  the  global  dimming  and  brightening  phenomenon  is  inherently  unclear 
because  it  is  unclear  for  any  given  location  and  time  whether  it  is  dimming  or  brightening 
that  is  occurring.  The  nature  of  global  dimming/brightening  influence  is  discussed  here 
with  regard  to  the  locations  included  in  this  study.  If  global  brightening  continues  to  occur 
in  North  America  over  the  next  few  years,  then  the  difference  between  TMY3-generated 
simulations  and  real-world  results  will  be  diminished  as  the  degree  of  HOMER  fuel  savings 
over-estimations  will  decrease.  On  the  other  hand,  if  global  dimming  occurs  in  North 
America  over  the  next  several  years,  the  disparity  between  TMY3-generated  simulations 
and  real-world  results  will  increase  as  TMY3  profiles  provide  increasingly  larger  over¬ 
estimates  of  available  solar  irradiance  compared  to  reality. 

6.3  Using  HOMER  to  Weigh  the  Benefits  of  Different 
Generator  Compositions  for  a  Micropower  System 

This  research  not  only  discussed  the  robustness  of  the  HOMER  model  but  also  demon¬ 
strated  ways  in  which  the  model  can  be  used  to  improve  power  planning  doctrine.  Chapter 
4  discussed  the  fuel  savings  and  maintenance  benefits  of  employing  multiple  generators  of 
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different  sizes  for  a  mieropower  system  rather  than  utilizing  only  a  single  generator  that 
ean  provide  the  peak  power  required  by  the  load.  The  first  group  of  HOMER  simulations 
diseussed  in  seetion  5.2.5,  those  where  the  mieropower  systems  did  not  include  GREENS 
systems,  supported  the  claim  that  increased  fuel  savings  can  be  achieve  by  employing  mul¬ 
tiple  different  sized  generators.  The  second  group  of  HOMER  simulations  discussed  in  that 
section,  those  where  mieropower  systems  did  include  GREENS  systems,  contradicted  the 
initial  claim  that  scenarios  with  only  a  single,  peak-load  generator  would  consume  more 
fuel  than  the  other  scenarios. 

These  surprising  findings  suggest  that  while  using  multiple  generators  results  in  substan¬ 
tial  fuel  savings  for  a  non-hybrid  mieropower  system,  those  benefits  disappear  once  a  PV 
array  and  battery  bank  system  are  added  to  the  mieropower  system.  John  Glassmire  from 
HOMER  EEC,  suggested  that  this  increased  relative  performance  of  the  single  peak- load 
generator  may  be  attributed  to  the  way  controllers  operate  in  selecting  when  to  charge  and 
when  to  discharge  batteries  based  on  the  controllers’  lack  of  knowledge  about  the  future. 

These  findings  can  prove  useful  for  decreasing  generator  maintenance  requirements  in  ad¬ 
dition  to  decreasing  fuel  consumption.  While  these  results  do  not  specifically  address  the 
question  of  wet  stacking  due  to  generators  operating  at  sub-par  loads,  it  can  be  inferred  by 
the  decreased  fuel  consumption  of  the  single  generator  scenarios  that  the  use  of  GREENS 
does  decrease  the  amount  of  time  that  a  single-peak  load  generator  is  operating  at  sub-par 
power  production  levels.  This  decrease  in  generator  time  spent  operating  at  sub-par  power 
production  levels  would  relate  to  decreased  wet-stacking  and  associated  maintenance  is¬ 
sues. 

6.3.1  Recommendation:  Use  Multiple  Different  Sized  Generators 
When  Not  Employing  GREENS  in  a  Mieropower  System 

It  is  recommended  that  mieropower  system  configurations  including  multiple  different 
sized  generators  be  employed  when  a  GREENS  or  other  such  battery  bank  and  solar  ar¬ 
ray  system  is  not  included  in  supporting  the  load.  It  is  also  recommended  that  the  benefit  of 
having  multiple  generators  when  a  GREENS  system,  or  similar  system,  be  reviewed  with 
regard  to  the  cost  in  transportation  space,  heavy  lift  engineering  equipment,  and  mainte¬ 
nance  requirements  associated  with  the  additional  generator(s).  While  this  research  did  not 
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definitively  show  that  the  inelusion  of  a  PV-battery  bank  system  makes  a  single  peak-power 
generator  preferable  to  multiple  generators,  it  does  provide  a  new  insight  for  weighing  the 
benefits  of  having  multiple  generators  relative  to  the  eosts. 


6.4  Fuel  Savings  Associated  with  Existing  versus  Grid-Tie 
Capable  Converters  for  the  GREENS  System 

Seetion  5.6  showed  mixed  results  from  employing  grid-tie  eapable  eonverters  with  the 
GREENS  system.  The  benefits  of  the  grid-tie  eonverters  relative  to  existing  eonverters 
seem  to  be  inversely  related  to  the  ratio  of  peak-load  level  to  PV  power  produetion  capaeity 
for  the  seenario.  Eittle  differenee  was  seen  between  the  fuel  savings  estimates  of  existing 
and  grid-tie  eapable  eonverters  for  the  3kW,  5kW,  and  lOkW  seenarios,  with  load  to  PV 
power  ratios  of  0.59,  0.35,  and  0.53,  respeetively.  Eor  the  30kW  and  60kW  seenarios,  how¬ 
ever,  the  fuel  savings  estimates  derived  the  simulated  use  of  a  grid-tie  eapable  eonverter 
were  higher  than  the  fuel  savings  estimates  derived  from  utilizing  the  existing  eonverter. 
These  seenarios  had  peak  load  to  PV  power  produetion  capaeity  ratios  of  0.29,  and  0.15, 
respectively.  Eor  the  30kW  scenarios,  50%  of  grid-tie  derived  fuel  savings  estimates  were 
higher  than  75%  of  those  derived  from  simulating  a  stand-alone  converter.  Eor  the  60kW 
scenarios,  75%  of  grid-tie  derived  fuel  savings  estimates  were  higher  than  75%  of  fuel 
savings  estimates  derived  from  simulating  a  stand-alone  converter. 

One  possible  explanation  for  this  relationship  is  that,  when  the  stand-alone  converter  is 
utilized,  the  solar  array  and  battery  bank  must  be  able  to  provide  all  of  the  power  required 
by  the  load  at  a  given  time  in  order  to  provide  any  power  to  the  load  at  that  time.  When 
the  load  size  is  increased  and  the  size  of  the  GREENS  is  not  increased  at  an  equivalent  rate 
(i.e.,  where  there  is  a  decrease  in  the  ratio  of  peak  load  to  PV  power  production  capacity) 
the  amount  of  time  that  the  GREENS  can  provide  power  to  the  load  is  decreased.  When  a 
grid-tie  converter  is  used,  on  the  other  hand,  the  GREENS  can  provide  power  to  the  system 
regardless  of  whether  it  can  cover  the  entire  load  by  itself.  This  means  that  at  certain  times 
the  GREENS  can  provide  power  simultaneously  with  a  small  generator  to  cover  a  load  that 
would  otherwise  require  the  use  of  a  larger  generator.  This  results  in  an  increased  utilization 
of  GREENS  in  grid-tie  converter  scenarios  and  a  related  increase  in  GREENS-derived  fuel 
savings. 
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6.4.1  Recommendation:  Develop  or  Acquire  Grid-Tie  Converters 
Only  for  Hybrid  Micropower  Systems  with  Peak  Power  Levels 
Greater  than  30kW 

The  question  of  whether  or  not  to  acquire  rugged,  grid-tie  capable  converters  does  not  have 
a  simple  answer.  While  giving  users  the  option  of  operating  individual  GREENS  in  either 
a  grid-tie  or  stand-alone  method  would  provide  for  increased  fuel  savings,  it  is  not  clear 
whether  that  would  be  worth  the  cost  of  developing  or  acquiring  a  grid-tie  capable  con¬ 
verter  for  all  hybrid  micropower  systems.  While  grid-tie  converters  exist  commercially,  it 
would  be  more  difficult  to  find  or  develop  grid-tie  converters  that  could  meet  the  ruggedness 
requirements  for  use  in  an  expeditionary  environment.  If  the  Marine  Corps  intends  to  use 
GREENS-like  systems  to  maximize  the  fuel  savings  of  larger  micropower  systems,  those 
with  peak  power  levels  of  30kW  or  higher,  then  it  is  more  likely  that  a  grid-tie  converter 
would  be  worth  the  costs.  If  the  Marine  Corps  intends  to  only  use  GREENS-like  systems  to 
support  smaller  micropower  systems,  as  seen  in  the  3kW,  5kW  and  lOkW  scenarios,  then 
it  is  more  likely  that  investing  in  grid-tie  converters  would  not  be  worth  the  costs. 


6.5  Future  Work 

6.5.1  Future  Work:  Analysis  of  the  Inter- Annual  Variability  of  Graham- 

Hollands  Generated  Solar  Irradiance  Profiles 

The  Graham-Hollands  algorithm  was  utilized  by  Newell  in  his  2010  thesis,  and  the  re¬ 
sulting  HOMER  fuel  savings  estimate  demonstrated  an  unacceptable  degree  of  error  when 
compared  to  a  HOMER  simulation  using  a  measured  solar  irradiance  profile.  While  it 
would  be  nice  to  always  be  able  to  use  TMY3  solar  irradiance  profiles,  the  TMY  profile 
requirement  of  historical  data  collection  precludes  the  constant  availability  of  TMY  pro¬ 
files  for  solar  resource  inputs.  Eor  this  reason,  it  is  necessary  to  quantify  the  variability  of 
Graham-Hollands  solar  irradiance  profiles. 

As  was  discussed  earlier,  this  thesis  identified  the  variability  of  HOMER  fuel  savings  es¬ 
timates  associated  with  the  inter-annual  variability  of  measured  solar  irradiance  profiles. 
While  this  provides  a  quantification  of  variability  in  HOMER  fuel  savings  estimates  gen¬ 
erated  using  TMY3  profiles,  these  measures  of  variability  cannot  be  used  when  Graham- 
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Hollands-generated  solar  irradiance  profiles  are  used  as  the  HOMER  solar  resource.  Anal¬ 
ysis  of  the  variability  between  Graham-Hollands-generated  solar  irradiance  profiles  and 
TMY3  solar  irradiance  profiles  would  allow  this  thesis’  quantification  of  inter- annual  vari¬ 
ability  due  to  solar  irradiance  profile  variability  to  be  appropriately  applied  for  simulations 
using  Graham-Hollands-generated  solar  irradiance  profiles.  This  could  be  accomplished 
with  an  experiment  that  quantifies  the  variability  of  fuel  savings  estimates  associated  with 
employing  Graham-Hollands-generated  solar  irradiance  profiles  and  compares  Graham- 
Hollands  generated  fuel  savings  predictions  and  fuel  savings  predictions  for  measured  an¬ 
nual  solar  irradiance  profiles. 

6.5.2  Future  Work:  Additional  Research  Into  the  Effects  of  Global 
Dimming  and  Brightening 

It  is  recommended  that  additional  research  into  global  dimming  and  brightening  be  con¬ 
ducted  to  increase  the  Marine  Corps’  understanding  of  the  phenomenon’s  implications  for 
solar  power  assets.  An  improved  understanding  of  the  causes  of  global  dimming  and  bright¬ 
ening  could  potentially  assist  in  improved  forecasting  of  solar  irradiance  profiles  for  future 
years  and  provide  associated  improvements  in  the  fidelity  of  HOMER  simulations. 

6.5.3  Future  Work:  Quantifying  the  Maintenance  Benefit  of  Different 
Generator  Composition  Scenarios  and  Investigating  Controller 
Improvements 

Euture  research  would  be  required  to  quantify  how  much  the  generator  compositions  of 
different  scenarios  relate  to  frequency  of  generators  operating  at  low  loads  relative  to  their 
peak  capacity.  This  would  allow  both  those  using  HOMER  operationally  and  those  using 
HOMER  to  improve  power  planning  doctrine  to  identify  which  composition  of  generators 
could  be  expected  to  result  in  the  least  amount  of  wet  stacking  and  maintenance  hours 
required  for  each  generator. 

It  is  also  recommended  that  future  research  be  conducted  into  the  way  in  which  controllers 
dictate  the  charging  and  discharging  of  batteries  for  a  hybrid  power  system.  If  advances 
have  been  made  in  controllers,  improving  the  controller  decision  making  process  for  charg¬ 
ing  and  discharging  of  batteries,  such  controller  improvements  could  potentially  provide  for 
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increased  efficiency  in  Marine  Corps  expeditionary  power  systems. 


6.5.4  Future  Work:  Continued  Validation  of  the  HOMER  Model 

This  research  was  conducted  under  the  assumption  that  Newell’s  assessment  of  the 
HOMER  model  as  providing  sufficient  fidelity  was  correct.  While  nothing  has  been  seen 
to  counter  this  assumption,  additional  experimentation  is  required  to  more  definitively  as¬ 
sess  the  model’s  fidelity.  More  specifically,  Newell’s  assessment  of  the  HOMER  model  as 
being  sufficiently  accurate  is  dependent  on  ascribing  the  last  17  percent  of  variability  seen 
in  the  initial  experiment  to  an  incorrect  derating  factor.  A  follow-on  experiment  was  not 
executed  to  validate  this  assumption  and  support  the  assessment  of  the  HOMER  model’s 
fidelity.  It  is  recommended  that  future  research  be  conducted  to  validate  the  assumption 
that  the  HOMER  model’s  over-estimate  in  fuel  savings  was  due  to  an  incorrect  solar  panel 
derating  factor  for  the  NFS  solar  panels  utilized  in  Newell’s  experiment. 

6.5.5  Future  Work:  Marine  Corps  Approach  to  Solar  Power  Asset 
Evaluation 

The  Marine  Corps  currently  evaluates  the  capabilities  of  solar  power  assets  by  assessing 
their  projected  performance  in  extreme  environments  rather  than  assessing  their  projected 
performance  across  the  range  of  most  likely  environments.  Marine  Corps  analysis  of  power 
generation  assets  includes  an  assessment  of  performance  in  extreme  heat  and  extreme  cold. 
While  evaluating  equipment  for  these  extreme  conditions  will  ensure  the  Marine  Corps 
has  equipment  that  is  capable  of  operating  in  any  clime  and  place,  it  does  not  help  in  the 
pursuit  of  optimal  system  performance  and  power  management  across  the  range  of  likely 
environments. 

An  analysis  of  equipment  should  include  a  balance  between  assessing  equipment  reliabil¬ 
ity  and  performance.  In  both  of  these  areas  one  train  of  thought  argues  that  quality  spread 
is  a  better  mode  of  assessment  than  assessing  performance  in  extreme  conditions.  To  do 
this,  one  could  identify  the  locations  where  Marine  Corps  deployments  are  most  likely  to 
occur  and  weight  the  performance  of  equipment  in  each  environment  according  to  the  per¬ 
ceived  likelihood  of  deployment  to  that  region  or  environment.  In  this  regard,  the  weighted 
performance  refers  to  both  measured  reliability,  in  terms  of  operating  time  and  mainte¬ 
nance  requirements,  and  energy  production  performance  of  the  equipment.  Methods  for 
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weighting  the  likelihood  of  Marine  Corps  deployment  can  include  values  such  as  national 
instability  rankings,  perceived  economic  and  political  challengers  to  U.S.  hegemony  in  re¬ 
gions  across  the  world,  and/or  the  inclusion  of  regions  falling  within  the  arc  of  instability. 
The  arc  of  instability  is  a  term  describing  a  region  of  the  world  "stretching  from  the  Andean 
region  of  Latin  America  across  Sub-Saharan  Africa,  the  Middle  East  and  the  Caucasus,  and 
through  the  northern  parts  of  South  Asia,"  which  is  deemed  to  contain  countries  particularly 
prone  to  political  instability  [44].  The  prevalence  of  "migration,  demographic  decline,  and 
economic  stress"  throughout  the  countries  within  the  arc  of  instability  increases  the  likeli¬ 
hood  of  conflict  in  the  region  [45].  This  thesis  has  shown  how  historical  solar  irradiance 
and  temperature  data  can  be  used  to  assess  equipment  performance  across  many  different 
environments  and  seasons.  By  basing  equipment  analysis  on  factors  such  as  these,  rather 
than  only  the  most  extreme  scenarios,  the  Marine  Corps  could  better  work  toward  maxi¬ 
mizing  its  expeditionary  energy  efficiency  and  effectiveness  across  the  range  of  operating 
environments. 

6.5.6  Future  Work:  HOMER  Robustness  Concerning  Load  Profile 
Variability 

The  need  for  an  analysis  of  how  variability  in  load  profiles  affects  HOMER  fuel  savings 
estimates  was  presented  in  Chapter  3.  While  the  effects  of  variability  in  load  profiles  was 
originally  intended  to  be  assessed  in  this  research,  a  lack  of  sufficient  measurements  of  load 
profiles  over  time  prevented  the  execution  of  that  analysis.  While  it  would  be  preferable  for 
future  analysis  of  this  factor  to  be  based  on  real  world  measures  of  load  profile  variability, 
a  DOE  that  produces  synthetic  load  profiles  with  specified  degrees  of  variability  could  also 
provide  valuable  insight  into  variability  in  the  load  profile  model  resource. 

6.5.7  Future  Work:  Spatial  Solar  Irradiance  and  Temperature  Vari¬ 
ability 

A  HOMER  user’s  decision  to  employ  either  a  TMY3  solar  irradiance  profile  or  a  profile  de¬ 
rived  with  the  Graham-Hollands  algorithm  is  a  cost-benefit  decision.  Eor  the  TMY3  option, 
greater  fidelity  of  the  profile  is  acquired  at  the  expense  of  limitations  on  locations  for  which 
the  method  can  be  employed.  TMY3  solar  irradiance  profiles  can  only  be  utilized  near  the 
sites  that  have  been  established  to  collect  measurements  and  NREE  has  synthesized  the 
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data  to  compile  a  TMY3  profile.  For  the  Graham- Hollands  option,  the  inereased  flexibility 
in  loeation  applieability  eomes  with  a  eost  in  the  fidelity  of  the  solar  irradianee  profile.  A 
elarifieation  of  this  eost-benefit  relationship  is  neeessary  to  identify  to  best  solar  resouree 
option  for  the  Marine  Corps  in  different  situations.  Sueh  a  elarifieation  eould  be  aehieved 
through  a  quantifieation  of  spatial  variability  effeets  on  solar  irradianee  and  temperature 
profiles. 

While  an  introduetion  to  spatial  solar  irradianee  variability  was  provided  in  seetion  3.2, 
greater  analysis  is  required  to  make  this  understanding  of  spatial  variability  useful.  Quan¬ 
tifieation  of  how  distanee  from  the  TMY3  measurement  loeation  affeets  the  variability  of 
TMY3  profiles  and  HOMER  fuel  savings  estimates  is  needed  to  help  determine  the  num¬ 
ber  and  spread  of  TMY3  measurement  sites  neeessary  to  make  the  Marine  Corps’  use  of 
TMY  solar  irradianee  profiles  for  the  HOMER  model  more  feasible  for  expeditionary  see- 
narios.  Similar  researeh  for  improving  the  Marine  Corps’  understanding  of  temperature 
spatial  variability  would  be  benefieial  for  the  same  reasons. 
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APPENDIX:  ReadMe  for  the  Interface  Built  to  Access 

HOMER  3.0  API 


README  for  managing  and  submitting  HOMER  experiments  on  reaper,  SEED  Center’s 
eluster 

Stephen  C.  Upton 

SEED  Center  for  Data  Earming 

12/3/2013 

modified:  5/13/2014 

Overview  Running  HOMER  on  our  eluster  required  development  of  several  software  arti- 
faets,  as  well  as  using  the  HOMER  3  software  eomponents.  The  HOMER  3  eomponents 
ineluded  a  HOMER  exeeutable  and  several  java  libraries.  The  other  software  artifaets 
described  below  were  written  to  configure  a  set  of  HOMER  "experiments"  and  gener¬ 
ate  supporting  files  in  order  to  run  HOMER  on  our  cluster.  Our  cluster  uses  the  open 
source  software  application  called  "HTCondor"  (http://research.cs.wisc.edu/htcondor/)  for 
distributing  the  runs  across  the  machines  in  the  cluster. 

The  other  two  pieces  of  software  that  were  written  are  called  HOMERRunner  and  HOME- 
REarmer.  HOMERRunner,  written  in  java,  uses  the  HOMER  3  APIs,  and  takes  a  single  set 
of  input  in  the  form  of:  (1)  a  HOMER  3  scenario;  (2)  a  city  name;  (3)  a  specific  year,  e.g., 
1961;  (4)  a  start  time;  (5)  an  end  time,  where  the  times  specify  a  duration  during  the  year; 
and  (6)  a  solar  irradiance  data  file  and  a  temperature  file,  both  associated  with  an  individual 
city  and  a  specified  year.  HOMERRunner  then  takes  those  as  input,  modifies  the  scenario 
input  file,  using  the  HOMER  3  API,  sets  the  latitude  and  longitude  of  the  city,  the  azimuth 
of  the  PV  component  based  on  city  data,  as  well  as  adjusting  the  solar  and  temperature 
input  to  conform  to  the  time  period  specified,  and  calls  the  HOMER  executable.  Once  the 
HOMER  executable  has  completed,  output  is  directed  to  a  CSV  (comma-separated  value) 
output  file,  with  the  file  name  constructed  from  the  input  settings. 

HOMEREarmer,  also  written  in  java,  is  a  wrapper  around  HOMERRunner  and  basically 
runs  HOMERRunner  over  multiple  inputs,  collected  together  in  what  we  call  an  "experi- 
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ment".  A  collection  of  experiments  is  considered  a  "run".  An  experiment  is  just  a  set  of 
input  settings  specified  in  a  CSV  input  file.  An  experiment  file,  or  a  subset  of  an  entire 
run,  is  considered  a  single  job.  A  single  job  is  then  run  on  a  single  processor  on  a  machine 
as  scheduled  by  the  HTCondor  software.  All  the  output  is  then  gathered  in  a  centralized 
location  and  post-processed. 

HOMER  Input:  1 .  DOE  file  an  experiment  file  that  lists  scenario  file,  solar  file,  temperature 
file,  start  time,  and  end  time,  in  that  order,  in  a  .csv  file  (the  city  and  year  are  extracted  either 
from  the  solar  file  or  the  temperature  file,  depending  on  the  type  of  experiment)  2.  scenario 
file  a  .hm3  HOMER  formatted  input  file  3.  solar  file  a  plain  text  file  with  8760  lines,  one 
per  hour,  representing  solar  irradiance  at  that  time  4.  temp  file  a  plain  text  file  with  8760 
lines,  one  per  hour,  representing  temperature  at  that  time 

Other  software  artifacts  (in  /Users/condor/models/HOMER  on  the  reaper  cluster): 

1.  clean_DOE_file.scala  usage:  scala  <pathtoscript>/clean_DOE_file.scala  <full  path  to 
DOE  file  name>  purpose:  removes  paths  from  DOE  file  for  scenario  file,  solar  file,  and 
temperature  file  entries 

2.  split_experiments. scala  usage:  scala  <pathtoscript>/split_experiments. scala  <full  path 
to  experiment  file>  <name  of  your  experiments  dir>  <source  to  get  year  from;either  "solar" 
or  "temp">  [optional  splitBy  -default  "scenario:city:year"]  purpose:  splits  the  DOE  file  into 
a  set  of  more  manageable  job  groups;  default  is  to  split  by  scenario,  city,  and  year 

3.  generate_submit_files. scala  usage:  scala  <pathtoscript>/generate_submit_files. scala 
<full  path  to  study  directory>  <name  of  your  experiments  dir,  relative  to  study  directory> 
<full  path  to  your  scenario  dir>  <full  path  to  your  data  dir>  <name  of  your  submit  dir,  rela¬ 
tive  to  study  directory>  [optional:  full  path  to  your  template  file;  default  is  studyDir/submit- 
template.dat]  purpose:  generates  the  condor  submit  files  for  a  HOMER  experiment  using 
the  submit-template.dat  file 

4.  submit-template.dat  file  used  as  template  to  generate  condor  submit  files  for  running  a 
HOMER  experiment 

5.  generate_ppsubmit_files. scala  usage:  scala  <pathtoscript>/generate_ppsubmit_files. scala 
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<full  path  to  study  directory>  [optional:  full  path  to  your  post-processor  submit  template 
file;  default  is  <homerdir>/ppsubmit-template.dat] 

6.  ppsubmit-template.dat  file  used  as  template  to  generate  eondor  submit  files  for  post- 
proeess  a  HOMER  experiment 

7.  eoneat.output.by.eity.R  R  seript  to  eoneatenate  the  output  from  a  set  of  experiments,  by 
eity 

8.  eombine.fuel.saved.R  ((ran  from  the  parent  direetory  of  the  PV/NoPV  direetories,  e.g., 
Input)  usage:  Rseript  <pathtoseript/eombine.fuel.saved.R  <use.for.year  (solar  I  temp)>  <so- 
lar  eategory  (R  I  T)>  <temp  eategory  (R  I  G)>  <inverter  type  (Existing  I  GridTie)>  purpose: 
coneatenates  the  output  from  an  experiment  group  that  has  PV  and  NoPV  direetories,  with 
output  assoeiated  with  running  a  seenario  that  has  PV  eomponents  and  one  that  has  NoPV 
eomponents,  but  is  the  same  seenario  otherwise.  Computes  the  fuel  saved,  i.e.,  the  differ- 
enee  in  fuel  usage  with  PV  eomponents  and  one  without  PV  eomponents.  This  is  caleulated 
using  the  same  seenario,  eity,  year,  start  time,  and  duration  settings. 

9.  eheek.missing. output. R  (ran  from  study  direetory)  usage:  Rseript  <pathtoseript/eheek. 
missing.output.R  <full  path  to  study  direotory>  <name  of  your  experiments  dir,  relative 
to  study  direotory>  <name  of  your  Output  dir,  relative  to  study  direotory>  <souree  to  get 
year  from;either  "solar"  or  "temp">  purpose:  eheeks  to  see  whieh  seenario/eity/year/start 
time/duration  experiment  output  is  missing  so  that  the  user  ean  re-run  if  there  was  an  issue. 

Proeedure:  1.  first  run  elean_DOE_file.seala  for  eaeh  DOE  file  to  remove  paths,  if  needed. 
Current  eonfiguration  assumes  that  the  seenario  files,  solar  files,  and  temperature  files  will 
be  loeated  in  the  exeeute  direetory  of  the  eondor  job.  If  desired,  ehanges  ean  be  made  to  the 
submit-template.dat  default  and  the  runHomerOnNode.bat  file  to  use  zipped  up  files;  just 
ensure  the  paths  in  the  DOE  files  refleet  what  those  paths  will  be 

2.  run  split_experiments.seala  to  divide  up  the  DOE  file  into  more  manageable  "ehunks"; 
default  is  to  group  by  seenario,  eity,  year,  i.e.,  in  the  generated  experiment  file  will  only  be 
the  jobs  for  a  speeifie  seenario  in  a  speeifie  eity  for  a  speeifie  year,  e.g..  Experiment l.hmS, 
Boulder,  1961 
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3.  run  generate_submit_files  to  generate  the  condor  submit  file  for  each  of  the  experiment 
files  generated  in  step  2. 

4.  go  into  the  submit  directory  and  submit  all  the  condor  jobs,  i.e.,  for  f  in  ‘ls‘  ;  do  con- 
dor_submit  $f ;  done 

Procedure  to  run  post-processor:  1.  run  generate_ppsubmit_files.scala  in  the  study  direc¬ 
tory  2.  go  into  the  ppsubmit  directory  and  submit  all  the  condor  jobs,  i.e.,  for  f  in  ‘ls‘  ;  do 
condor_submit  $f ;  done 

Procedure  to  compute  fuel  saved  1 .  run  combine.fuel.saved.R  (ran  from  the  parent  directory 
of  the  PV/NoPV  directories,  e.g..  Input) 

Procedure  for  running  the  regressions: 

1.  open  up  R  in  the  combined  directory  (where  your  combined  file  for  the  experiment  is 
located. 

2.  source("  /models/HOMER/fit.regressions.R") 

3.  library (plyr) 

4.  rr  <-  read.csv("nameoffile  with  all  data"),  e.g.,  read.csv("Existing.withPercent.withPV.csv") 

5.  rg  <-  ddply(rr,.(InverterType,scenario,city,StartTime,duration,solar,temp),  function(x) 
fit.model(x,"percent.fuel. saved", "percent.fuel. saved")) 

6.  rg.pv  <-  ddply(gg,.(InverterType,scenario,city,StartTime,duration,solar,temp),  func- 
tion(x)  fit.model(x,"PV.production","PV. production")) 

7.  hh<-join(rg,rg.pv,by=c("InverterType", "scenario", "city",  "StartTime", "duration", "solar", 
"temp")) 

8.  write.csv(hh,"Existing.withPercent.withPV.regressions.csv", 
quote=EAESE,row.names=EAESE) 
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